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An Introduction to Atmospheric Thermodynamics 


This new edition is a self-contained, concise but rigorous book 
introducing the reader to the basics of the subject. It has been 
brought completely up to date and reorganized to improve the 
quality and flow of the material. 

The introductory chapters provide definitions and useful 
mathematical and physical notes to help readers understand the 
basics. The book then describes the topics relevant to atmospheric 
processes, including the properties of moist air and atmospheric 
stability. It concludes with a brief introduction to the problem of 
weather forecasting and the relevance of thermodynamics. Each 
chapter contains worked examples to complement the theory, as 
well as a set of student exercises. Solutions to these are available 
to instructors on a password protected website at 
www. Cambridge. org/9780521696289. 

The author has taught atmospheric thermodynamics at 
undergraduate level for over 20 years and is a highly respected 
researcher in his field. This book provides an ideal text for short 
undergraduate courses taken as part of an atmospheric science, 
meteorology, physics, or natural science program. 

Anastasios A. Tsonis is a professor in the Department of 
Mathematical Sciences at the University of Wisconsin, Milwaukee. 
His main research interests include nonlinear dynamical systems 
and their application in climate, climate variability, predictability, 
and nonlinear time series analysis. He is a member of the 
American Geophysical Union and the European Geosciences 
Union. 



We have explained that the causes of the elements are four: 
the hot, the cold, the dry, and the moist. In every case, 
heat and cold determine, conjoin, and change things. Thus, 
hot and cold we describe as active, for combining is a sort 
of activity. Things dry and moist, on the other hand, are the 
subjects of that determination. In virtue of their being acted 
upon, they are thus passive. 

Aristotle, Meteorology , Book IV 
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PREFACE 


This book is intended for a semester undergraduate course in 
atmospheric thermodynamics. Writing it has been in my mind 
for a while. The main reason for wanting to write a book like 
this was that, simply, no such text in atmospheric thermodynam- 
ics exists. Do not get me wrong here. Excellent books treating 
the subject do exist and I have been positively influenced and 
guided by them in writing this one. However, in the past, atmo- 
spheric thermodynamics was either treated at graduate level or 
at undergraduate level in a partial way (using part of a general 
book in atmospheric physics) or too fully (thus making it diffi- 
cult to fit it into a semester course). Starting from this point, my 
idea was to write a self-contained, short, but rigorous book that 
provides the basics in atmospheric thermodynamics and prepares 
undergraduates for the next level. Since atmospheric thermody- 
namics is established material, the originality of this book lies in 
its concise style and, I hope, in the effectiveness with which the 
material is presented. The first two chapters provide basic defini- 
tions and some useful mathematical and physical notes that we 
employ throughout the book. The next three chapters deal with 
more or less classical thermodynamical issues such as basic gas 
laws and the first and second laws of thermodynamics. In Chapter 
6 we introduce the thermodynamics of water, and in Chapter 7 
we discuss in detail the properties of moist air and its role in 
atmospheric processes. In Chapter 8 we discuss atmospheric sta- 
bility, and in Chapter 9 we introduce thermodynamic diagrams 
as tools to visualize thermodynamic processes in the atmosphere 
and to forecast storm development. Chapter 10 serves as an epi- 
logue and briefly discusses how thermodynamics blends into the 
weather prediction problem. At the end of each chapter solved 
examples are supplied. These examples were chosen to comple- 
ment the theory and provide some direction for the unsolved 
problems. 


IX 
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Finally, I would like to extend my sincere thanks to Ms Gail 
Boviall for typing this book and to Ms Donna Genzmer for draft- 
ing the figures. 


Anastasios A. Tsonis 
Milwaukee 



CHAPTER ONE 


Basic definitions 


• Thermodynamics is defined as the study of equilibrium states of a 
system which has been subjected to some energy transformation. 
More specifically, thermodynamics is concerned with transforma- 
tions of heat into mechanical work and of mechanical work into 
heat. 

• A system is a specific sample of matter. In the atmosphere a par- 
cel of air is a system. A system is called open when it exchanges 
matter and energy with its surroundings (Figure 1.1). In the 
atmosphere all systems are more or less open. A closed system is 
a system that does not exchange matter with its surroundings. 
In this case, the system is always composed of the same point- 
masses (a point-mass refers to a very small object, for example 
a molecule). Obviously, the mathematical treatment of closed 
systems is not as involved as the one for open systems, which 
are extremely hard to handle. Because of that, in atmospheric 
thermodynamics, we assume that most systems are closed. This 
assumption is justified when the interactions associated with 
open systems can be neglected. This is approximately true in the 
following cases, (a) The system is large enough to ignore mixing 
with its surroundings at the boundaries. For example, a large 
cumulonimbus cloud may be considered as a closed system but a 
small cumulus may not. (b) The system is part of a larger homo- 
geneous system. In this case mixing does not significantly change 
its composition. A system is called isolated when it exchanges 
neither matter nor energy with its surroundings. 

• The state of a system (in classical mechanics) is completely speci- 
fied at a given time if the position and velocity of each point-mass 
is known. Thus, in a three-dimensional world, for a system of N 
point-masses, 6 N variables need to be known at any time. When 

I 
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Figure l.l 

In an open system mass 
and energy can be 
exchanged with its 
environment. A system is 
defined as closed when it 
exchanges energy but not 
matter with its 
environment, and as 
isolated if it exchanges 
neither mass nor energy. 



N is very large (like in any parcel of air) this dynamical defini- 
tion of state is not practical. As such, in thermodynamics we are 
dealing with the average properties of the system. 

If the system is a homogeneous fluid consisting of just one com- 
ponent, then its thermodynamic state can be defined by its 
geometry, by its temperature, T, and pressure, p. The geometry 
of a system is defined by its volume, V, and its shape. How- 
ever, most thermodynamic properties do not depend on shape. 
As such, volume is the only variable needed to characterize geom- 
etry. Since p, V, and T determine the state of the system, they 
must be connected. Their functional relationship f(p,V,T) = 0 
is called the equation of state. Accordingly, any one of these vari- 
ables can be expressed as a function of the other two. It follows 
that the state of a one-component homogeneous system can be 
completely defined by any two of the three state variables. This 
provides an easy way to visualize the evolution of such a sys- 
tem by simply plotting V against p in a rectangular coordinate 
system. In such a system, states of equal temperature define an 
isotherm. Any other thermodynamic variables that depend on the 
state defined by the two independent state variables are called 
state functions. State functions are thus dependent variables and 
state variables are independent variables; the two do not differ 
in other respects. That is why in the literature there is hardly 
any distinction between state variables and state functions. State 
variables and state functions have the property that their changes 
depend only on the initial and final states, not on the particular 
way by which the change happened. If the system is composed of 
a homogeneous mixture of several components, then in order to 
define the state of the system we need, in addition to p, V, T, the 
concentrations of the different components. If the system is non- 
homogeneous, we must divide it into a number of homogeneous 
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parts. In this case, p, V, and T of a given homogeneous part are 
connected via an equation of state. 

For a closed system, the chemical composition and its mass 
describe the system itself. Its volume, pressure, and tempera- 
ture describe the state of system. Properties of the system are 
referred to as extensive if they depend on the size of the sys- 
tem and as intensive if they are independent of the size of the 
system. An extensive variable can be converted into an intensive 
one by dividing by the mass. In the literature it is common to 
use capital letters to describe quantities that depend on mass 
(work, W, entropy, S ) and lower case letters to describe inten- 
sive variables (specific work, w, specific heat, q). The mass, m, 
and temperature, T, will be exceptions to this rule. 

An equilibrium state is defined as a state in which the system’s 
properties, so long as the external conditions (surroundings) 
remain unchanged, do not change in time. For example, a gas 
enclosed in a container of constant volume is in equilibrium if 
its pressure is constant throughout and its temperature is equal 
to that of the surroundings. An equilibrium state can be stable, 
unstable, or metastable. It is stable when small variations about 
the equilibrium state do not take the system away from the equi- 
librium state, and it is unstable if they do. An equilibrium state 
is called metastable if the system is stable with respect to small 
variations in certain properties and unstable with respect to small 
changes in other properties. 

A transformation takes a system from an initial state i to a final 
state /. In a (p. V ) diagram such a transformation will be rep- 
resented by a curve / connecting i and /. We will denote this 
as i — T -$ /. A transformation can be reversible or irreversible. 
Formally, a reversible transformation is one in which the succes- 
sive states (those between i and /) differ by infinitesimals from 
equilibrium states. Accordingly, a reversible transformation can 
only connect those i and f states which are equilibrium states. 
It follows that a reversible process is one which can be reversed 
anywhere along its path in such a way that both the system 
and its surroundings return to their initial states. In practice a 
reversible transformation is realized only when the external con- 
ditions change very slowly so that the system has time to adjust 
to the new conditions. For example, assume that our system is 
a gas enclosed in a container with a movable piston. As long as 
the piston moves from i to / very slowly the system adjusts and 
all intermediate states are equilibrium states. If the piston does 
not move slowly, then currents will be created in the expanding 
gas and the intermediate states will not be equilibrium states. 
From this example, it follows that turbulent mixing in the atmo- 
sphere is a source of irreversibility. If a system goes from i to 
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/ reversibly, then it could go from f to i again reversibly if the 
same steps were followed backwards. If the same steps cannot 
be followed exactly, then this transformation is represented by 

another curve /' in the (p, V) diagram (i.e. / — -A i) and may or 
may not be reversible. In other words the system may return to 
its initial state but the surroundings may not. Any transforma- 
tion i — > / — > i is called a cyclic transformation. Given the 
discussion above we can have cyclic tansformations which are 
reversible or irreversible (Figure 1.2). A transformation i — -» / 
is called isothermal if I is an isotherm, isochoric if I is a constant 
volume line, isobaric if / is a constant pressure line, and adiabatic 
if during the transformation the system does not exchange heat 
with its surroundings (environment). Note and keep in mind for 
later that adiabatic transformations are not isothermal. 

• Energy is something that can be defined formally (we have to 
wait a bit for this), but its concept is not easily understood by 
defining it. We all feel we understand what is meant by energy, 
but if we verbally attempt to explain what energy is we will get 
upset with ourselves. At this point, let us just recall that for 
a point-mass with a mass m p moving with speed v in a uni- 
form gravitational field g , Newton’s second law takes the form 


Figure 1.2 

(a) A reversible cyclic 
process; (b) an 
irreversible cyclic process. 


P 



f 


1= reversible 
/' = reversible 



(b) 


/= reversible 
/' = irreversible 
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d(K + P)/dt = 0 where K = m p v 2 /2 , P = m p gz, t is the 
time, and z is the height. K is called the kinetic energy and P is 
called the potential energy. The total energy of the point-mass 
E = K + P is, therefore, conserved. If we consider a system of 
N interacting point-masses that may be subjected to external 
forces (other than gravity), then the total energy of the system 
is the sum of the kinetic energy about the centre of gravity of 
all point-masses (internal kinetic energy), the kinetic energy of 
the centre of gravity, the potential energy due to interactions 
between the point-masses (internal potential energy), and the 
potential energy due to external forces. The sum of the internal 
kinetic and internal potential energy is called the internal energy 
of the system, U. A system is called conservative if dU/dt = 0 
and dissipative otherwise. For systems considered here we are 
interested in the internal energies only. 




CHAPTER TWO 


Some useful mathematical 
and physical topics 


2.1 Exact differentials 


If z is a function of the variables x and y, then by definition 

dz= {t*)f x+ {%), iv (2 ' 1) 

where dz is an exact differential. Now let us assume that a quan- 
tity Sz can be expressed according to the following differential 
relationship 

Sz = Mdx + Ndy (2-2) 

where x and y are independent variables and M and N are functions 
of x and y. If we integrate equation (2.2) we have that 


/ 


Sz = 


Mdx + 



Since M and N are functions of x and y, the above integration 
cannot be done unless a functional relationship f(x, y) = 0 between 
x and y is chosen. This relationship defines a path in the (a;, y) 
domain along which the integration will be performed. This is called 
a line integral and its result depends entirely on the prescribed path 
in the (a;, y) domain. If it is that 


M = 




then equation (2.2) becomes 


(2.3) 


The right-hand side of the above equation is the exact or total 
differential dz. In this case Sz is an exact differential. If we 
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now integrate Sz from some initial state i to a final state / we 
obtain 

/ <52=/ dz = z(x f ,y f ) - z(Xi,yi). (2.4) 

J i J i 

Clearly, if <52 is an exact differential its net change along a path 
i — > f depends only on points i and / and not on the particular 
path from i to /. We say that in this case 2 is a point function. 
All three state variables are exact differentials (i.e. Sp = dp , ST = 
dT , SV = dV). It follows that all quantities that are a function of 
any two state variables will be exact differentials. 

If the final and initial states coincide (i.e. we go back to the initial 
state via a cyclic process), then from equation (2.4) we have that 


Sz = 0. 


(2.5) 


The above alternative condition indicates that <52 is an exact differ- 
ential if its integral along any closed path is zero. At this point we 
have to clarify a point so that we do not get confused later. When we 
deal with pure mathematical functions our ability to evaluate f Sz 
does not depend on the direction of the closed path, i.e. whether 
we go from i back to i via i ——t f — T —> i or via i — f — ^-A i 
(see Figure 1.2). When we deal with natural systems we must view 
the condition j> Sz = 0 in relation to reversible and irreversible 
processes. If, somehow, it is possible to go from i back to i via 
i — A / —A i but impossible via i -—t f —A i (for example, 
when /' is an irreversible transformation), then computation of <52 
depends on the direction and as such it is not unique. Therefore, 
for physical systems, the condition § Sz = 0 when <52 is an exact 
differential applies only to reversible processes. 

Note that since 

d dz d 2 z d 2 z d dz 

dy dx dydx dxdy dx dy 


it follows that an equivalent condition for <52 to be an exact 
differential is that 


d M _ dN 
dy dx ' 


( 2 . 6 ) 


Equations (2.3)-(2.6) are equivalent conditions that define 2 as a 
point function and subsequently <52 as an exact differential. If a 
thermodynamic quantity is not a point function or an exact differ- 
ential then its change along a path depends on the path. Moreover, 
its change along a closed path is not zero. Such quantities are path 
functions. For a path function the thermodynamic processes must 
be specified completely in order to define the quantity. For the rest 
of this book an exact differential will be denoted by dz whereas a 
non-exact differential will be denoted by <52. Finally, note that if <52 
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is not an exact differential and if only two variables are involved, a 
factor A (called the integration factor) may exist such that A Sz is 
an exact differential. 


2.2 Kinetic theory of heat 

Let us consider a system at a temperature T, consisting of N 
point-masses (molecules). According to the kinetic theory of heat, 
these molecules move randomly at all directions traversing rectilin- 
ear lines. This motion is called Brownian motion. Because of the 
complete randomness of this motion, the internal energies of the 
point-masses not only are not equal to each other, but they change 
in time. If, however, we calculate the mean internal energy, we will 
find that it remains constant in time. The kinetic theory of heat 
accepts that the mean internal energy of each point-mass, XJ, is 
proportional to the absolute temperature of the system (a formal 
definition of absolute temperature will come later; for now let us 
denote it by T), 

U = constant x T (2.7) 

Let us for a minute assume that N = 1. Then the point has 
only three degrees of freedom which here are called thermodynamic 
degrees of freedom and are equal to the number of independent vari- 
ables needed to completely define the energy of the point (unlike 
the degrees of freedom in Hamiltonian dynamics which are defined 
as the least number of independent variables that completely define 
the position of the point in state space) . The velocity v of the point 
can be written as 

2 2 , 2,2 

V =v x + v y + v z . 

Because we only assumed one point, then the total internal energy 
is equal to its kinetic energy. Thus, 


or 

U x 

where to each component corresponds one degree of freedom. 
According to the equal distribution of energy theorem, the mean 
kinetic energy of the point, U , is distributed equally to the three 
degrees of freedom i.e. U x = U y = U z . Accordingly, from equation 
(2.7) we can write that 

U i = AT, i = x,y,z 

where the constant A is a universal constant (i.e. it does not depend 
on the degrees of freedom or the type of the gas). We denote this 


Uy = 


m p Vy 


U; = 
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constant as k/2 where k is Boltzmann’s constant (k = 1.38 x 10 -23 
J K _1 ) (for a review of units, see Table A1 in the Appendix). There- 
fore, the mean kinetic energy of a point with three degrees of 
freedom is equal to 



( 2 . 8 ) 


m v v 2 3 kT 

2 ~~ ~T~ 


The theorem of equal energy distribution can be extended to N 
points. In this case, the degrees of freedom are 3N and 


or 


or 


n 2 
y- rn p vf 

^ 2 


SNkT 

2 


1 N 2 

y- m p vj 

N ^ 2 




(2.9) 


where { rnp 2 v ) is the average kinetic energy of all N points. Note that 
the above is true only if the points are considered as monatomic. 
If they are not, extra degrees of freedom are present that corre- 
spond to other motions such as rotation about the center of gravity, 
oscillation about the equilibrium positions, etc. 

The kinetic theory of heat has found many applications in the 
kinetic theory of ideal gases. An ideal gas is one for which the 
following apply: 


(a) the molecules move randomly in all directions and in such a 
way that the same number of molecules move in any direction; 

(b) during the motion the molecules do not exert forces except 
when they collide with each other or with the walls of the 
container. As such the motion of each molecule between two 
collisions is linear and of uniform speed; 

(c) the collisions between molecules are considered elastic. This 
is necessary because otherwise with each collision the kinetic 
energy of the molecules will be reduced thereby resulting in a 
temperature decrease. Also, a collision obeys the law of specular 
reflection (the angle of incidence equals the angle of reflection); 

(d) the sum of the volumes of the molecules is negligible compared 
with the volume of the container. 


Now let us consider a molecule of mass m p whose velocity is v and 
which is moving in a direction perpendicular to a wall (Figure 2.1). 
The molecule has a momentum P = m p v. Since we accept that the 
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Figure 2.1 

A molecule of mass m p 
moving with a velocity v 
and hitting a surface S. If 
this collision is assumed 
elastic and specular, then 
the change in momentum 
is 2 m p v. 



collision is elastic and specular, the magnitude of the momentum 
after the collision is —m p v. Thus the total change in momentum is 
m p t) — (— m p v) = 2m p v. According to Newton’s second law, F = 
dP/dt. If we consider N molecules occupying a volume V we can 
calculate the change in momentum dP of all molecules in the time 
interval dt by multiplying the change in momentum of one molecule 
(2m p v) by the number of molecules dN that hit a given area S on 
the wall, i.e. 


dP = 2m p vdN (2.10) 

Note that here we have assumed that all molecules have the same 
speed. The number dN of molecules hitting area S during dt is equal 
to the number of molecules which move to the right and which are 
included in a box with base S and length vdt. Since the motion 
is completely random, we can assume that | of the molecules will 
be moving to the right, | will be moving to the left, and | will be 
moving along the directions of the other two coordinates. Since the 
volume of the box is Svdt and the number of molecules per unit 
volume is N/V then the number of molecules inside the box is 

— Svdt. 

Accordingly, the number of molecules moving to the right and 
colliding with S is 

dN = ^ -Svdt (2.11) 


From equations (2.10) and (2.11) it follows that 
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Recalling the definition of pressure, p (pressure = force/area), and 
Newton’s second law we obtain 


P = 


1 N 

3 V 


777.pl> 2 . 


The above formula resulted by assuming that all molecules move 
with the same speed. This is not true, and because of that m p v 2 in 
the above equation should be replaced by the average of all points, 
(777.pi> 2 } . We thus arrive at 


P = 


IN 

3b 


(?77 p 1> 2 ). 


(2.12) 


It can easily be shown that by combining equations (2.9) and (2.12), 
we can derive an equation that includes all three state variables: 


pV = NkT. (2.13) 

Equation (2.13) provides the functional relationship of the equation 
of state f(p,V,T) and it is called the ideal gas law. More details 
follow in the next chapter. 



CHAPTER THREE 


Early experiments and laws 


At the end of Chapter 2 we derived theoretically the equation of 
state or the ideal gas law. This law was first derived experimentally. 
The relevant experiments provide many interesting insights about 
the properties of ideal gases and confirm the theory. As such a little 
discussion is necessary. 


3.1 The first law of Gay-Lussac 

Through experiments Gay-Lussac was able to show that, when pres- 
sure is constant, the increase in volume of an ideal gas, dV, is 
proportional to the volume V 0 that it has at a temperature (mea- 
sured in the Celsius scale) of 9 = 0°C and proportional to the 
temperature increase, dO : 


dV = aV 0 d9. (3.1) 

The coefficient a is called the volume coefficient of thermal expan- 
sion at a constant pressure and it has the value of 1/273 deg -1 for 
all gases. The physical meaning of a can be understood if we solve 
equation (3.1) for a: 

1 dV 
a ~ dfiVo' 

From the above equation it follows that if we increase the temper- 
ature of an ideal gas by 1 °C, while we keep the pressure constant, 
the volume will increase by 1/273 of the volume the gas occupies 
at 0°C. By integrating equation (3.1) we obtain the relationship 
between V and 9: 



V - V 0 = aV 0 9 


13 



14 


3 EARLY EXPERIMENTS AND LAWS 


Figure 3.1 

Graphical representation 
of the first law of 
Gay-Lussac. 



or 

V = V o (l + a0). (3.2) 

This is a linear relationship and its graph is shown in Figure 3.1. 


3.2 The second law of Gay-Lussac 


Again through experimentation Gay-Lussac was able to show that 
for a constant volume the increase in pressure of an ideal gas, dp, 
is proportional to the pressure p 0 that it has at a temperature of 
0 °C and proportional to the increase in temperature dO : 

dp = j3p 0 d6. 


The coefficient (3 is the pressure coefficient of thermal expansion at 
a constant volume and has the value 1/273 deg” 1 for all gases: 


P 


1 dp 
dO p 0 


The above formula indicates that an increase in temperature by 
1 °C (while V is constant) results in an increase in pressure by 
1/273 of the pressure the gas had at 0°C. 

As in the first law again it follows that (Figure 3.2) 


p = p 0 (l + PO). 


(3.3) 


Application 

The second law of Gay-Lussac can easily explain why, in the winter, 
heating a house at a much greater temperature than the outside 
temperature does not increase the pressure enough to break the 
windows. A difference of 20 °C between the inside and the out- 
side air increases the pressure inside the house by 7.3%. Glass can 
withstand such pressure changes easily. 


3.4 ANOTHER FORM OF THE GAY-LUSSAC LAWS 
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Figure 3.2 

Graphical representation 
of the second law of 
Gay-Lussac. 


P 



-273 


0 
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3.3 Absolute temperature 

From equation (3.2) it follows that if we extrapolate to 9 = —273 °C, 
then V=0. This means that if we were able to cool an ideal gas 
to —273 °C, while keeping the pressure constant, the volume would 
become zero. Similarly, from equation (3.3) it follows that if we 
were able to cool an ideal gas to the temperature of —273 °C, while 
keeping the volume constant, the pressure would become zero. This 
temperature of — 273 °C we call absolute zero. Up to now for the 
measurement of temperature the Celsius scale has been used, start- 
ing from the temperature “zero Celsius” . If we extend the Celsius 
scale to the absolute zero (i.e. —273 °C), then the temperature mea- 
sured from the absolute zero is called absolute temperature, T. This 
defines a new scale called the Kelvin scale (T = 273 + 9) . 


3.4 Another form of the Gay-Lussac laws 

Using the absolute temperature we can present the Gay-Lussac laws 
as follows. From Figure 3.3 it follows from the similarity of triangles 
ABC and AB'C' that his first law can be expressed as 


Similarly from Figure 3.4, his second law can be expressed as 


V T 

for p = constant, — = — 

F V' T' 


(3.4) 



(3.5) 


Verbally stated, the volumes of an ideal gas under constant pressure 
and the pressures under constant volume are proportional to the 
absolute temperatures. 
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Figure 3.3 

Expressing Gay-Lussac’s 
first law for two states C 
and C' . 



T 


r 


Figure 3.4 

Expressing Gay-Lussac’s 
second law for two states 
C and C'. 



V 


3.5 Boyle’s law 

While Gay-Lussac’s laws provide the change in volume or pressure 
as a function of temperature, Boyle’s law provides the change in 
pressure as the volume varies at a constant temperature. The law 
is expressed as follows: 

pV = p'V'. (3.6) 


3.6 Avogadro’s hypothesis 

The formal definition of a mole is the amount of substance which 
contains the same number of particles (atoms, molecules, ions, or 
electrons) as there are in 12 grams of 12 C. This number of atoms is 
equal to N = 6.022 x 10 23 and it is known as Avogadro’s number. 
It is obtained under the condition that one mole of carbon weighs 
12 grams. In general, a substance of a weight equal to its molecular 
weight contains one mole of the substance. Thus 25 grams of water 
contain 27/18 = 1.5 mole of water. In 1811 the Italian physicist 


3.7 THE IDEAL GAS LAW 
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and mathematician Amedeo Avogadro proposed his hypothesis that 
the volume of a gas, V, is directly proportional to the number of 
molecules of the gas, N', 


V = aN' 


where a is a constant. To relate to this hypothesis, simply imagine 
inflating a balloon. The more the air you pump into it, the bigger 
the volume. It follows that at constant temperature and pressure, 
equal volumes of gases contain the same number of molecules. For 
one mole N' = N. Since N is the same for one mole of any gas, 
Avogadro’s hypothesis (actually a law by now) can be stated simply 
as: A mole of any gas at constant temperature and pressure occupies 
the same volume. This volume for the standard state T 0 = 0°C, 
Po — 1 atmosphere has the value (see equation (2.13)) 

Vt 0 , P0 = 22.4 liters mol -1 
= 22 400 cm 3 mol” 1 . 


3.7 The ideal gas law 


Let us consider an ideal gas at a state p, V, T which is heated under 
constant volume to a state Pi,V, T' (Figure 3.5). Then, according 
to Gay-Lussac’s second law 

V 

Pi =p Y- 

If subsequently we keep the temperature constant and increase the 
volume to V' the gas goes to a state p'V'T'. Then, according to 
Boyle’s law, 

p'V' = Pl V. 


By combining the above two equations we obtain 


V 

p'V'=p-V 


or 


pV _ p'V ' 
~T ~ ~T r ' 


(3.7) 


This is called the Boyle-Gay-Lussac law and it indicates that if we 
change the pressure, the volume, and the temperature the value of 
PV/T remains constant. Thus, equation (3.7) can be written as 


pV = AT 


(3.8) 


where A is, according to equation (3.7), a constant depending on 
the type and mass of the gas. The dependence on the mass follows 
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Figure 3.5 p ▲ 

The three steps used to 
derive the ideal gas law 
experimentally. 



> 

V 


from the fact that if we consider twice the mass and keep the pres- 
sure and the temperature constant, then the volume will double 
(Avogadro’s hypothesis). Therefore, the value of pV/T (i.e. A) will 
double. Accordingly, we can write A = mR where m is the mass of 
the gas and R is a constant independent of the mass but dependent 
on the gas. This constant is called the specific gas constant. We can 
then write equation (3.8) as 

pV = mRT. (3.9) 

Since m = Nm p (where N is the total number of molecules each 
having mass m p ), equation (3.9) is identical to equation (2.13) with 
R = k/m p . Equation (3.9) is thus the equation of state or ideal 
gas law which is now derived experimentally. Equation (3.9) can be 
written as p = pRT where p = m/V is the mass density. By defining 
the specific volume a = lip, the ideal gas law takes the form 

pa = RT. 

If M denotes the molecular weight of the gas then the number of 
moles is n = m/M. It follows that 


pV = nMRT 


3.8 A LITTLE DISCUSSION ON THE IDEAL GAS LAW 
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or 

pV = nR*T (3.10) 


or 



T 


where R* = MR. Considering that R = k/m p it follows that 
R* = Mk/rrip = mk/nm p = Nk/n. Now recall that the number 
of molecules in one mole, N A , is equal to 6.022 x 10 23 (Avogadro’s 
number). Then, N = nN A . In this case R* = N A k , which is the 
product of two constants. This new constant is called the universal 
gas constant and its value is 8.3143 JK -1 mol -1 . 


3.8 A little discussion on the ideal gas law 


• Because the ideal gas law relates three variables (rather than 
two), we have to be careful when we interpret changes in one of 
the variables. For example, temperature increases when pressure 
increases but only if the volume (density) increases (decreases) 
or remains constant or decreases (increases) by a smaller amount 
compared with pressure. Let us consider the ideal gas law pV = 
mRT or p = pRT and differentiate it: 

pdV + V dp = mRdT 


dp = RT dp + RpdT. 


We can easily see that if dp > 0 and dV > 0 then dT > 0. Also, if 
dp > 0 and dV < 0 and \Vdp\ > \pdV\ then dT > 0. And so on! 
It follows that cold air is denser than warm air only if pressure 
remains constant or if its change does not offset the temperature 
difference. 

Similarly, if we consider the ideal gas law in the form 


R * 

pa = — -1 
y M 


or 


R* 

v = P„ T 


we can reason that moist air (in which some of the dry air of 
molecular weight 28 g mol -1 has been replaced by water vapor of 
molecular weight 18 g mol -1 ) is less dense than dry air at more 
or less the same pressure and temperature. This explains some 
interesting statistics in the American game of baseball where 
more home runs occur when the weather is hot and humid. In 
a warmer and more moist environment where p = constant, the 
density of air is smaller. Therefore, the ball has less resistance. 
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For a path of 400 feet or so the effect can be significant, resulting 
in higher chances for a home run. 

• Recall equations (2.9) and (3.2). According to equation (2.9) 
T = OK when v = 0. Because of this one might interpret the 
zero absolute temperature as the temperature where all motion 
in an ideal gas ceases. On the other hand, when equation (3.2) 
is extrapolated to the temperature at which V = 0 it results 
in ^solute = (— 4)°C = — 273 °C = OK. While this tempera- 
ture corresponds to vanishing volumes, it is not necessarily the 
temperature at which motion ceases. If both V = 0 and v = 0 
are true when T = OK, then pressure can be anything, but this 
does not invalidate the ideal gas law, both sides of which will 
equal zero. To avoid confusion, we must keep in mind that these 
issues emerge when these equations are extrapolated to states 
where one cannot assume that a gas behaves as an ideal gas. 
At those states, equations (2.9) and (3.2) simply do not mean 
much. 


3.9 Mixture of gases - Dalton’s law 


Consider a mixture of two gases occupying a volume V and con- 
sisting of Nx molecules of gas 1 and N 2 molecules of gas 2. The 
total pressure on the walls will be the result of all the collisions, 
i.e. the collisions by the molecules of gas 1 and the collisions by the 
molecules of gas 2. Thus we can write equation (2.13) as 


or 


(Nx + N 2 ) 
p= v kT 


P = 


NxkT N 2 kT 


V 


V 


The first term on the right-hand side of the above equation is 
exactly the pressure that gas 1 would create if all its Nx molecules 
occupied the volume V, i.e. the partial pressure of gas 1. The same 
is valid for the second term. We thus conclude that the total pres- 
sure is the sum of the partial pressures. This expresses Dalton’s 
law which states that for a mixture of K components, each one of 
which obeys the ideal gas law, the total pressure, p, exerted by the 
mixture is equal to the sum of the partial pressures which would 
be exerted by each gas if it alone occupied the entire volume at the 
temperature of the mixture, T : 


K 

p=^2pi- 

i = 1 
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If the volume of the mixture is V and the mass and molecular 
weight of the zth constituent are m; and M t respectively, then for 
each constituent 


Pi = 


M< V ' 


Applying Dalton’s law, we then can write that for a mixture 


or 


p=E 


R*T rrii 


P = 


R*T 

~v~ 


E 


rrii 

Mi' 


Since the total mass of the mixture m = 
equation becomes 


P = 


E K rm 

i = i m7 


V 


y K ■. 

Z-^i= 1 ‘ 


th e above 


or 


R*T 

P= 

a 


E K rrii 
i = 1 Mi 

E K 

i= i m i 


(3.11) 


In order for the mixture to obey the ideal gas law, it must satisfy 
the equation 


RT _ R*T 
a aM 


(3.12) 


where M is the mean molecular weight of the mixture. By compar- 
ing equations (3.11) and (3.12) we see that this would be possible 
if 


E k 


£ 


(3.13) 


This provides the proper way to compute the mean molecular 
weight for a mixture and indicates that the mixture follows the 
ideal gas law as long as there is no condensation (if there is, then 
some rrii s may not remain constant). 


Examples 

(3.1) Determine the mean molecular weight of dry air. 

For our planet, the lowest 25 km of the atmosphere is made 
up almost entirely by nitrogen (N 2 ), oxygen (0 2 ), argon (A), 
and carbon dioxide (C0 2 ) (75.51%, 23.14%, 1.3%, and 0.05% 
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by mass, respectively). Thus, for dry air the mean molecular 
weight is 

— _ to N 2 + m 02 + m A + m c o 2 

M ~ m N 2 , rnp 2 mA m C Q 2 
m N2 _r m 02 _r m a _r M C o 2 

or 

— 75.51 + 23.16 + 1.3 + 0.05 , , 

M — 75.51 i 23.14 . 1.3 , O-OiT § m0 ^ 

28.02 ’ + ~ 32.0 39.94 44.01 


or 


M = 28.97 gmol' 1 


or 


M = 0.02897 kg mob 1 . 


It follows that the specific gas constant for dry air is 
R d = R*/M = 287 J kg -1 K -1 . 


(3.2) Determine the mean molecular weight of a mixture of dry air 
saturated with water vapor at 0°C and I atmosphere pressure. 
The partial pressure of water vapor at 0°C is 6. 1 I mb. 

The mixture consists of dry air and water vapor. If p d is the 
pressure due to dry air and p v the pressure due to water 
vapor, then the pressure of the mixture is p = p d + p v . If 
we recall that the number of moles n = m/M, then we can 
write equation (3.13) as 


M = 


n i M i Eg 

EL 


, TliA/li ^ ^ Tli 
1 1 


n 


i= 1 


(3.14) 


where rii and n are the corresponding number of moles of the 
constituents of the mixture and the total number of moles 
in the mixture. Assuming that both dry air and water vapor 
are ideal gases we have that 


pV = nR*T for the mixture 

PiV = HiR*T for each of the constituents. 


From these two equations it follows that 

p n 
Pi n , ' 

Combining the above relationship and equation (3. 14) yields 

K 

M —Mi. 
i = 1 P 
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In our case i = 1,2. Thus, 

M = —M,. + —M 2 

P P 

or changing 1 — > cl (dry air) and 2 — > v (water vapor) 

M = —M d + —M v 
P P 

or 

M = + — M v 

P P 

or 

M = 28.9 g mol -1 . 

The above estimated mean molecular weight is not very 
much different from the mean molecular weight of dry air. 
The difference becomes more apparent for higher temper- 
atures. For example for T = 35 °C the partial pressure of 
water vapor is 57.6 mb and M = 28.3 g mol -1 . 


(3.3) Two containers A and B of volumes V A = 800 cm 3 and V B = 
600 cm 3 , respectively, are connected with a tube that closes and 
opens by a hinge. The containers are filled with a gas under 
pressures of 1 000 mb and 800 mb, respectively. If we open the 
connection what will the final pressure be in each of the con- 
tainers? Assume that the temperature remains constant. 

Once the connection is open, each gas will expand to fill the 
total volume V = V A +V B thereby equalizing the difference in 
pressure between the two containers. Thus the final pressure 
in each container will be the same. Let us first assume that 
container B is empty. Since the gas in container A expands 
at a constant temperature we have (Boyle’s law) 


p f Vf =PiVi 


where i and / stand for initial and final. It follows that 
Pf = 571 mb. 

Similarly if we assume that container A is empty the final 
pressure in both containers will be 

p f = —— = 343 mb. 
v f 

Since Pf(p'f) is the pressure the gas in A(B) would exert if it 
occupied the total volume V A + V B ,Pf and p' f can be consid- 
ered as partial pressures of two gases. Then from Dalton’s 
law it follows that the final pressure in each container should 
be 914 mb. 
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(3.4) If Boyle had observed y/ pV = constant, what would the equa- 
tion of state for an ideal gas be? In this case calculate the 
temperature of a sample of nitrogen which has a pressure of 
800 mb and a specific volume of I200cm 3 g _l . How does this 
temperature differ from that estimated when the correct law 
pV = constant is observed? Can you explain the difference? 

If we follow the procedure to arrive at equation (3.7) but 
with the new Boyle’s law, we have that 

T 

Pi =P~f 


y/p'V' = y/piV. 


By combining the above equations we have that 


or 


yfp'V' = y/p 



VFV 

Vt 


VpV 

Vt 


= constant. 


As we know, if we keep the temperature and pressure con- 
stant and double the mass, the volume will double. As such, 
here as well the ratio y/pV/y/T will double. Then we can 
write that y/pV/y/T = mR where m is the total mass and 
R is the specific gas constant. Thus, the ideal gas law in this 
case will be 


Solving for T we have that 


ypf = 1.143 


or 


T = 1.07K. 

Under normal circumstances (i.e. when pa = ( R*/M)T ), 
we find that T = 323.3 K which makes more sense. The 
tremendous difference is due to the fact that the square root 
in yjpa = constant introduces two corrections to the ideal gas 
law (both pressure and temperature appear under a square 
root) having the net result of significantly reducing T . 
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Problems 

(3.1) What is the mass of dry air occupying a room of dimensions 
3x5x4 m at p=l atmosphere and T = 20 °C? (72.3 kg) 

(3.2) Graph the relationship V = f(T) p = constant from absolute 
zero up to high temperatures, for two samples of the same 
gas which at 0°C occupy volumes of 1000 and 2000 cm 3 , 
respectively. 

(3.3) In a 2-D coordinate system with axes the absolute temper- 
ature and volume, graph (a) an isobaric (constant pressure) 
change of one mole of an ideal gas for p = 1 atmosphere, 
and (b) the same when p = 2 atmospheres. 

(3.4) Determine the molecular weight of the Venusian atmosphere 
assuming that it consists of 95% C0 2 and 5% N 2 by volume. 
What is the gas constant for 1 kg of such an atmosphere? 
(43.2 g mol- 1 , 192.5 J kg" 1 K" 1 ) 

(3.5) If p = 1 atmosphere and T = 0°C how many molecules are 
there in 1cm 3 of dry air? (2.6884 x 10 19 molecules) 

(3.6) Given the two states p, V, T and p', V',T', define on a (p, V) 
diagram the state pi , V, T' that was used to arrive at the 
ideal gas law (section 3.7). 

(3.7) An ideal gas of p,V,T undergoes the following successive 
changes: (1) it is warmed under a constant pressure until 
its volume doubles, (2) it is warmed under constant volume 
until its pressure doubles, and (3) it expands isothermally 
until its pressure returns to p. Calculate in each case the 
values of p, V, T and plot all three changes in a (p, V) 
diagram. 




CHAPTER FOUR 


The first law of 
thermodynamics 


4.1 Work 

As we have already mentioned, in atmospheric thermodynamics we 
will be dealing with equilibrium states of air. If a system (par- 
cel of air, for example) is at equilibrium with its environment no 
changes take place in either of them. We can imagine the “shape” 
of the parcel remaining unchanged in time. If the pressure of the 
surroundings changes, then the force associated with the pressure 
change will disturb the parcel thereby forcing it away from equilib- 
rium. In order for the parcel to adjust to the pressure changes of 
the surroundings, the parcel will either contract or expand. If the 
parcel expands we say that the parcel performs work on the envi- 
ronment and if the parcel contracts we say that the environment 
performs work on the parcel. By definition, if the volume change is 
dV then the incremental work done, dW , is 

dW = pdV. 

Accordingly, when the system changes from an initial state i to a 
final state / the total work done, either by the system or on the 
system, is 



The above equation indicates that the work done is given by an area 
in a (p, V) diagram (Figure 4.1) [or a (p, a) diagram, where a is the 
specific volume]. If dV > 0 (the system expands) it follows that 
W > 0 and if dV < 0 (system contracts) it follows that W < 0. 
Thus, positive work corresponds to work done by the system on 
the environment and negative work corresponds to work done to 
the system by the environment. 

Now let us consider a situation where the system expands 
through a reversible transformation from i to / and then contracts 
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Figure 4. 1 

The shaded area gives the 
work done when a system 
changes from an initial 
state i to a final state /. 


Figure 4.2 

The shaded area gives the 
work done during a cyclic 
reversible transformation. 



P 



from f to i along exactly the same path in the ( p , V) diagram. 
Then, the total work done will be 


W = / pdV + / pdV = / pdV — / pdV = 0. 


If, however, the system expands from i to f and then contracts 
from f to i but along a different reversible transformation (different 
path), then the total work done would be (see Figure 4.2) 


w = 

[ f pdV 

+ 

[ pdV 


.J i 

1 

J f 


= area under curve 1 — area under curve 2 


— 7 ^ 0 


where A ili2 i is the area enclosed by the two paths. It follows that 
§ dW = j> pdV ^ 0 and thus dW is not an exact differential, which 
means that work is not a state function. As such it depends on the 
particular way the system goes from i to /. Because of this from 
now on we will denote the incremental change in work as SW not 
as dW. 


4.2 DEFINITION OF ENERGY 


29 


If we write dV = dAds (where dA is an area element and ds a 
distance element) we have that SW = pdAds = Fds = Fvdt or 
that 


or 


or 


or 


where v denotes the velocity of the parcel and K is the kinetic 
energy of the parcel. It follows that the work done and the kinetic 
energy are related. The last equation indicates that one way by 
which a thermodynamic system can exchange energy with its envi- 
ronment is by performing work. The other is through transfer of 
heat. More on that will follow soon. According to (4.1) the units 
for work are those of energy. Thus, the unit for work in the MKS 
system is the joule which is defined as J = N m where the newton 
N = kg ms -2 . In the cgs system, the unit is the erg which is defined 
as erg = dyncm where dyn = gems -2 . It follows that 1 joule = 
10 7 erg. 


5W 

dt 


= Fv 


dv 


SW 

—r~ = m—v 
dt dt 


SW 


1 


dt dt V 2 
SW dK 


7 / 


7 i 


(4.1) 


4.2 Definition of energy 

The first law of thermodynamics expresses the principle of con- 
servation of energy for thermodynamical systems. The idea here 
is that energy cannot be created or destroyed. It can only change 
from one form to another. In this regard if during a transforma- 
tion the energy of the system increases by some amount, then this 
amount is equal to the amount of energy the system receives from 
its surroundings in some other form. 

Let us consider a closed system contained in an adiabatic enclosure. 
In this case the energy of the system, U, is equal to the sum of the 
potential and kinetic energy of all its molecules. The sum of the 
energies of all molecules depends on the state of the system at a 
given moment (i.e. on the values p,V,T) but obviously is indepen- 
dent of past states. It follows that the internal energy of the system 
depends on the state in which it exists but not on the way it arrived 
at that state. Thus, in a transformation i — > /, A U = Uf — [/;. 
Then, for a cyclic transformation, j> A U = 0 which means that for 
an infinitesimal process dU is an exact differential. If no external 
forces are acting upon the system (i.e. the system is at equilibrium 
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Figure 4.3 

The work done by an 
external force on a 
system that is taken 
adiabatically from a 
reference state O or a 
reference state O' to a 
state A. In this case the 
work done depends on 
the initial and final states, 
not on the particular path 
from an initial state to a 
final state. 



with its environment), then its energy remains the same (A U = 0). 
If some external force acts upon the system taking it from state i 
to state / then the conservation of energy principle dictates that 

U f -Ui = -IT ad (4.2) 

where — IT ad is the work done adiabatically by the external force 
on the system (+TF ad is the work done by the system). Note that 
equation (4.2) implies that the work done in this case depends only 
on states i and /, not on the particular transformation i — > /. This 
is contrary to what we proved previously - that the work done 
is not an exact differential and that it depends on the particular 
transformation i — > /. The discrepancy is due to the fact that here 
we are dealing with adiabatic transformations. Only in this case 
does the work done not depend on the particular transformation. 
Assume now a standard state O of this system where by definition 
Uo = 0 and a state A where U a yf 0. Assume further that through 
some external influences we take the system from state O to state 
A. Then according to equation (4.2) the energy of state A is 

U A = -Wf. (4.3) 

Apparently, equation (4.3) (and hence any definition of energy) 
depends on the standard state O. If instead of O we had chosen 
another standard state O' we would have obtained U' A . In this case 
(see Figure 4.3) we would have (since IT ad depends on initial and 
final states only) 

—W A d = -wf - W$ 

or 

U A = U' A + Uo' 
or 

U a -U' a = Uq> 

where Uo’ is the energy of the standard state O' and thus is a 
constant. It follows that had we chosen a different standard state, 
U' A would differ by a constant (additive) amount. This constant is 
an essential feature of the concept of energy which has no effect on 
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the final results when we are working with differences (changes) of 
energy and not with actual amounts. 


4.3 Equivalence between heat and work done 

Consider a sample of water at a state i where the temperature is 
Tj and at a state / where the temperature is Tf > T). In how many 
different ways can we take the system from i to /? Obviously we 
can simply warm up the system. In this case, and as long as Tf 
is low enough for evaporation to be negligible, the volume of the 
system remains unchanged and therefore no work is done by the 
external forces. This is the first way. The second way is to raise 
the temperature via friction. Just imagine a set of paddles rotating 
inside the system. In this case, we will observe that as long as the 
paddles continue to rotate the temperature of the water increases. 
Since the water resists the motion of the paddles we must perform 
mechanical work in order to keep the paddles moving until Tf is 
reached. 

Thus, the work performed by the system depends on whether we 
go from i to / by means of the first or the second way. If we 
assume that the principle of the conservation of energy holds, then 
we must admit that the energy transmitted to the water in the 
form of mechanical work from the rotating paddles is transmitted 
to the water in the first way in a non-mechanical form called heat. 
It follows that heat and (mechanical) work are equivalent and two 
different aspects of the same thing: energy. In the first way, where 
no work is performed, the change of energy can be expressed as 

AU = Qw= o 

where Qw = o is the amount of heat received. In the second way (the 
experiment can be thought of as taking place in some adiabatic 
enclosure) the change in energy can be expressed as 

AU = -W ad . 

When both work and heat are allowed to be exchanged, we replace 
the above equations by the general expression 

A U + W=Q. (4.4) 

The quantity Q w =o required to raise lg of water from 14.5 °C to 
15.5 °C is by definition equal to 1 cal. In this case W ad turns out to 
be equal to 4.185 J. This value is known as the mechanical equiv- 
alent of heat. Note that since 6W is not an exact differential, SQ 
is not either. Thus, for an infinitesimal process the above equation 
takes the form 


dU + SW = 6Q. 


(4.5) 
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Note that for a cyclic transformation dU = 0 and thus W = Q. 
Equation (4.5) is the mathematical expression of the first law of 
thermodynamics. Since 6W = pdV we can write the first law as: 


or for unit mass 


dU + pdV = 6Q 

(4.6) 

du + pda = 5q. 

(4.7) 


4.4 Thermal capacities 

Since dU is an exact differential we can write it as: 

dT 




or as: 


or as: 


dU = 


dU = 


f dU 
\dT 

dU 

dp 


dT + 


fdU 

\dp 


dp 


dp - 


fdU 


dV. 


By combining equations (4.6) and (4.8) we get 


fdU 

\dT 


dT- 


p + 


fdU 

Vay 


dV = 5Q 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


By combining equations (4.6) and (4.9) and using the fact that 


dV 


we arrive at 

fdU 
lar 


dV 

dT 


dp 


dT - 


dV = ITT d P + dT 


dV\ 


dT) , 


dU 

dp 


dV 

T P \dp 


By combining (4.6) and (4.10) we get 


dp = SQ. 
(4.12) 


dU \ w 

dj) dP ' 


dU\ 

~dV) 


dV = 6Q. 


(4.13) 


We define as thermal capacity, C, the ratio SQ/dT and we denote 
it as C v if the thermal capacity is measured at a constant V and as 
C p if it is measured at constant p. According to this definition if in 
a process at constant volume there is no change of the physical and 
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chemical state of a homogeneous system, then the heat absorbed is 
proportional to the variation in temperature 

SQ = C v dT 


or 


5Q = c v mdT 


where to is the mass of the system and c v is the specific heat 
capacity at constant volume (i.e. c v = C v /m). Another specific 
heat is the molar specific heat c Vm = C v / n where n is the number 
of moles. Similarly, at constant pressure 

SQ = C p dT 


or 


SQ = c p mdT. 


From equation (4.11) if we set dV = 0 we obtain 


_ SQ_(dU\ 
v - ' dTj 


and 


dT 


c v — 


du 

dT 


From equation (4.12) if we set dp = 0 we obtain 

~P 


^ 


dT 


dV 

dT 


(4.14) 


By defining here a new state function, the enthalpy, H = U + pV 
(h = u + pa), it follows from the above equation that 


and 



(4.15) 


Again here we can define the molar specific heat at constant pres- 
sure as c pm = C p /n. Also, since pV = nR*T, for ideal gases the 
enthalpy is a function of two state variables (U, T) and as such it 
can be expressed as an exact differential. Note that for isobaric pro- 
cesses dH = dU + pdV = SQ , i.e. the change in enthalpy equals the 
amount of heat absorbed or given away. 
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Figure 4.4 

The experimental design 
of Joule’s proof that 

U = U (T). 



4.5 More on the relation between U and T (Joule’s law) 

As we have stated earlier, any thermodynamical quantity can be 
expressed as a function of two state variables. As such we can 
assume that U = f(V,T). In Chapter 2 we discussed the fact that 
the kinetic theory of heat accepts that U is proportional to the 
absolute temperature. This indicates that U is a function of only 
one state variable, U = f(T). An experimental proof of this was 
provided by Joule, who performed the following experiment. He 
constructed a calorimeter and placed in it a container having two 
chambers, A and B, separated from each other by a stopcock (see 
Figure 4.4). He filled A with an ideal gas and evacuated B. After 
the whole system settled to a thermal equilibrium (as indicated 
by the thermometer in the calorimeter), Joule opened the stopcock 
thus allowing the ideal gas to flow from A to B until pressure every- 
where equalized. He observed that this did not cause any change in 
the temperature reading. From the first law we thus have Q = 0 or 
from equation (4.4) 

A U = —W. 

Now note that since the volume of the system (the container) does 
not change, the system performs no work. As such W = 0 and 
A U = 0. But, you may ask, how about the gas inside chamber 
A?. Its volume increases as it expands to fill container B. What 
is the effect of that? Since there was no temperature change and 
no variation in energy during the process we must conclude that a 
variation in volume at a constant temperature produces no variation 
in energy. In other words U ^ f(V) which means that U is only a 
function of T : 

U = U(T). (4.16) 


Later in problem 6.2 you will be asked to prove this mathematically. 
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More on thermal capacities 

Because of equation (4.16), ( 8U/dT) v = dU/dT. As such we can 
write equation (4.14) as 


Cy 


dU 

dT 


du 

° v = dT ' 

Similarly, H = U+pV = U + nR*T = H{T) and 


(4.17) 


r _dH 
p dT 


dh 


Cp dT' 


(4.18) 


In the atmospheric range of temperatures, c v and c p are nearly 
constant. From equation (4.17) we then have that 


U = 


C v dT + constant 


or 

U = C V T + constant. 

Since only differences in internal energy (and thus in enthalpy) 
are physically relevant, we can set the integration constant to zero 
and obtain 

U W CyT. 

Similarly H « C p T,u « c v T, and h « c p T. Equations (4.17) and 

(4.18) can be combined to obtain 


C p — C v = nR* (4.19) 

which implies that 

c p Cy R , c pm Cy m R . (4.20) 

Recall now from Chapter 2 (equation (2.8)) that for a monatomic 
gas consisting of N point-masses the total internal energy is 

U = ^ NkT . 

It follows that 

° v= % = \ Nk = \ mR = l nMR = 2 nRt ' 
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and 

3 p* 



Cv m ^ ^ 

, Cy 

Similarly, 

„ dH 

dU 



~ ~dT 

and 


5 


c = 

2 *’ 


3 n . 
2 171 


c p 2 ^ ' 

For a diatomic gas these values change to 

Cy 2 U ^ 7 ^Vm 2 ^ 5 2 ^ 


and 


C p = - nR *, 



7 

2 


i?. 


Note that often in the literature the values of C v and c Vm or C p 
and c pm are assumed the same. This is correct only in the case of 
one mole. From the above it follows that 


C p 

Cy 


Cpm 

Cy m 


Cp 

Cy 


We will denote these ratios as 7 . 

Dry air is considered a diatomic gas. According to the above 
relationships we estimate that for dry air 


7a = 1-4 

c vd = 718 J kg " 1 K " 1 = 171 cal kg -1 K -1 
Cp d = 1005 J kg " 1 KT 1 = 240 cal kg -1 K " 1 
R d = 287 J kg -1 K _1 . 

According to the above definitions we can write the first law as: 

C v dT + pdV = SQ (4.21) 


or 


c v dT + pda = 5q 

Considering that C v = C p — nR* and that pV = nR*T , equation 
(4.21) becomes 

( C p - nR*)dT + ( nR*dT - Vdp) = 5Q 


or 

C p dT - Vdp = SQ. (4.22) 

For convenience Table 4.1 summarizes all the different expressions 
of the first law. 
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Table 4.1. Expressions of the first law 


For a gas of mass m For unit mass 


dU + SW = SQ 
dU + pdV = SQ 
C v dT + pdV = SQ 
C p dT — V dp = SQ 


du + Sir = Sq 
du + pda = Sq 
CydT + pda = Sq 
CpdT — adp = Sq 


4.6 Consequences of the first law 

Let us look at the form of the first law in the following special cases. 


, , , - . . T=constant n 

• Isothermal transformations: i — > j 

During such transformations dT = 0 and as such from equation 
(4.21) it follows that SQ = 
is 

Q = 

or 

o 

or 

V f 

Q = nR*T In -f. 

* i 

• Isochoric transformations: i — > j 

In this case dV = 0 and from the first law it follows that 


SW. Then the amount of heat exchanged 
[ f SW= [ f pdV 


= nR*T / 


SQ = dU = C v dT 


or 


AU = Q = C v 



dT = C v (T f 


Ti). 


• Isobaric transformations: 


p=constant 


/ 


or 


In this case it follows that 
SQ = C p dT 
SW = pdV 


or 


Q = Cp{T f - Ti), 
W = p(V f -Vi). 



38 


4 THE FIRST LAW OF THERMODYNAMICS 


and 

A U = C p (T f -T i )-p(V f -V i ) 


Cyclic transformations 


As we know, in cyclic transformations j> dU = 0. Then, from the 
first law it follows that 



• Adiabatic transformations - Poisson’s relations 


During an adiabatic process there is no exchange of heat between 
a system and its surroundings by virtue of temperature differences 
between them. Thus, SQ = 0. Then the first law can be expressed 
as 

dU = ~SW 


or 


or 


or 


or using pV = m.RT 


or 


C v dT = - pdV 

CpdT = Vdp 

dT _ V dp 
~T ~ ~C V ~T 

dT mR dp 

~t = ~c; 7 

dT nR* dp 
T = 


or using C v — C v = nR* 


or 


dT 


\ dp 

~T ~ ' 

\ C r 

) P 

dT 

fl- 1 ) 

1 

~T ~ 

V 7/ 

V ' 


Integrating equation (4.23) yields 


1 


or 


In T = I 1 ) In p + ln(constant) 


In X = In p t + ln(constant) 


(4.23) 
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Figure 4.5 

The various relationships 
between adiabats, 
isotherms, isobars, and 
isochores in partial and in 
the complete state spaces. 


P 



(a) 





(c) 




or 


T = constant x 


-y — 1 
p ■y 


or 


1 — 7 

Tp i = constant. 


(4.24) 


Through the ideal gas law we can then derive the following 
equivalent expressions of equation (4.24): 

TV 7-1 = constant (4-25) 

pV 1 = constant. (4-26) 

Equations (4.24)-(4.26) express Poisson’s relations for adiabatic 
processes. Since for dry air 7 = 1.4 > 0 equation (4.25) proves 
something that we all more or less know: that as an air parcel 
rises and expands (i.e. V increases), its temperature must decrease. 
For adiabatic expansion, since 6Q is zero this decrease is due to 
expansion work that the parcel performs on the environment. 
From the ideal gas law it follows that isotherms are given by the 
equation 

pV = constant. 


Thus on a (p, V ) diagram isotherms are equilateral hyperbolas 
(Figure 4.5(a); an equilateral hyperbola is one whose asymptotes 
are perpendicular to each other). From equation (4.26) we see that 
the adiabats on a (p, V ) diagram will also be equilateral hyper- 
bolas but since 7 > 1 the adiabats will be steeper than the 
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isotherms. Obviously in a ( p , V) diagram isobaric and isochoric 
transformations are represented by straight lines on the p and V 
axis, respectively. On a ( p , T ) diagram the isochores are given by 
the equation 

rr nR*T 

V = = constant 

P 


or 


Tp 1 = constant 


which is the equation of a straight line. From equation (4.24) it 
also follows that in a ( p , T) diagram the adiabats are again hyper- 
bolas (Figure 4.5(b)). On a (T,V) diagram (Figure 4.5(c)) the 
isobars are straight lines (as before TV -1 = constant) and the adi- 
abats are hyperbolas (see equation (4.25)). The complete ( p,V,T ) 
space and transformations are shown in Figure 4.5(d). Since they 
require no heat exchange between a parcel and its environment (i.e. 
they require perfect isolation) , adiabatic processes are idealizations 
that never happen in nature. However, for processes that happen 
quickly, energy transfer is not very significant. Such processes are 
thus adiabatic to a good approximation. 


• Polytropic transformations 

During an isothermal transformation the system can exchange 
heat with its environment but its temperature remains constant. 
This, however, is an ideal situation. In reality the equiliza- 
tion of temperatures does not happen instantaneously. Simi- 
larly, an adiabatic transformation is an ideal situation since it 
requires perfect heat insulation. The same comments can be 
made for isobaric and isochoric transformations. Thus, the rela- 
tions pV = constant, pV 1 = constant, Tp -1 = constant, and 
TV -1 = constant are special cases of the general relationship 
pV v = constant where can take on any value. Such transforma- 
tions are called polytropic transformations (in Greek polytropic 
means multi-behavioral). If we assume that rj = 7 /e where e 
is a constant, we can work our way backwards from equation 

1 — 7? 

pV v = constant or Tp = constant to arrive at 

dT = / _ e\ dp 
T V 7 ) P‘ 

Considering that 7 = the above equation can be expressed 
as 

dT f C p — CVe\ dp 
~T = ^ ~C V ) J 
or 

dT ( C p — C v — C v e + C v \ dp 
~T = ^ C, ) J 
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or 

dT ( C p — C v \ dp C v { 1 — e) dp 

~T = V c p JJ + 7 7 

or 

dT nR* dp (1 — e) dp 

T C p p + 7 p 

or 

dT V dp (1 — e) dp 

~T = C~ P Y + 7 7 

or 

C p dT = Vdp + ( 1 ~ CpT dp 
IV 
or 

C v dT = Vdp - C iTzl 1 77 dp _ 

7 nR * 

Then, according to the first law, the second term on the right- 
hand side is the heat exchanged, which is zero if e = 1 (i.e. 77 = 7, 
adiabatic process), greater than zero (i.e. the system absorbs heat 
from the environment) if 77 < 7 (e > 1) and dp < 0, or 77 > 7 
(e < 1) and dp > 0, and less than zero (i.e. the system gives 
away heat to the environment) if 77 > 7 and dp < 0, or 77 < 7 and 
dp > 0. From the above equation we also find that for e = 0 (i.e. 
77 = 00) nR* dT = Vdp , which owing to the ideal gas law implies 
dV = 0. For e = 7 (i.e. 77 = 1), it reduces to C p dT = 0, while for 
e = 00 (i.e. 77 = 0) it dictates that dp = 0. As such a polytropic 
process pV 17 = constant reduces to an isobaric, to an isothermal, 
to an adiabatic, and to isochoric process for 77 = 0, 1,7, and 00, 
respectively. 

Having said this we must add that in the atmosphere over a 
rather wide range of motions the timescale for an air parcel (our 
system) needed to adjust to environmental changes of pressure, 
and to perform work, is short compared with the corresponding 
timescale of heat transfer. For example, above the boundary layer 
and outside the clouds the timescale of heat transfer is about two 
weeks whereas the timescale of displacements that affect a parcel 
is of the order of hours to a day. It thus appears that adiabatic 
approximations are good approximations for many atmospheric 
phenomena. 

• Dry adiabatic lapse rate 

From equation (4.24) we have 

7— 1 

T = constant -pi. 

By taking the logarithmic derivative of the above equation we 
have 

7 — 1 

din T = d In p 

7 
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or 


or 


dT 7 — 1 dp 
T 7 P 

1 dT 7 — 1 1 dp 
T dz 7 p dz 


(4.27) 


For large-scale motions the hydrostatic approximation states 
that pressure gradient force balances the force due to gravity. 
Therefore, 

77 = -P9- (4-28) 


For a parcel of air rising in the atmosphere, equation (4.27) is 
valid as long as its ascent is an adiabatic process. For the same 
parcel equation (4.28) is valid only if dp/dz experienced by the 
parcel is equal to that of the large scale (surroundings). We 
assume this is more or less true but in equation (4.28) p rep- 
resents the density of the surroundings which from the ideal gas 
law is 


P = 


P 

RT S 


(4.29) 


where T s denotes the temperature of the surroundings. In this 
case by combining equations (4.27)-(4.29) we have that 


dT _ _ 7 - 1 g T 
dz 7 RT S 

Recall that by definition an adiabatic ascent or descent is one in 
which there is no exchange of energy with its surroundings by 
virtue of temperature differences between them. Accordingly, if 
the parcel remains dry (i.e. no condensation takes place) we may 
define from the above equation the dry adiabatic lapse rate, T d , as 
the atmospheric temperature profile such that the temperature 
of the parcel is always at the temperature of its surroundings. 
For such a profile T = T s and as such for dry air 

T d = —- 7 — = = — = 9.8°Ckm -1 . (4.30) 

dz 7 R d c pA 

This value is larger than the observed average decrease of envi- 
ronmental temperature with altitude. The difference is mainly 
due to the fact that in the derivation of T d we neglected the 
presence of water vapor. When vapor is present we have new 
lapse rates, which we will discuss in a later chapter. 

Recalling that h = u + pa and that h = c pd T we can write 
equation (4.30) as 


dh 
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or 

~~ 7 ~(h + gz) = 0 (4.31) 

dz 

where h + gz is defined as the dry static energy. The above 
equations indicate (1) that enthalpy of a parcel decreases as it 
rises adiabatically because it is doing gravitational work, and (2) 
that the static energy (the sum of enthalpy and the gravitational 
potential energy) is conserved in an adiabatic motion (ascent or 
descent). Enthalpy is thus the specific internal energy plus the 
term pa that accounts for the work done by the parcel on the 
surroundings. 


• Potential temperature 

Let us assume that there is a variable 9 defined by the relation 
9 = ATp~ 0 where A and (3 are constants. By taking the loga- 
rithm on both sides we get In 9 = In A + In T — (3 In p. Then by 
differentiating we arrive at 


d In 9 = d In T — f3d In p. 


Consider now the first law in the form c p dT—adp = 5q and divide 
it by T . Then with the help of the ideal gas we obtain 

d In T — — d In p = . 

Cp CpJ. 


For f3 = R/c p the above two equations combine to give 

d In 9 = (4.32) 

c p T 

For adiabatic processes 5q = 0 and thus d In 9 = 0. As such we 
should expect that in adiabatic processes there should exist a 
measure 9 which is conserved. This measure can be defined as 
follows. 

From equation (4.24) we have 


v. k 

p k p^ 

where k = (7— l)/j = l—c v /c p = R/c v = 0.286 (for dry air). The 
state (T 0 ,p 0 ) can be taken as a reference state. As such we can 
choose p 0 = 1000 mb and denote the corresponding temperature 
as T 0 . The above equation is expressed as 


T 0 = T 



k 


or 


T 0 =p k 0 Tp~ k . 


(4.33) 
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If we compare the above expression to 9 = ATp ~ 13 we see that 
T 0 = 9 for A = pg and /? = /;. We call 9 the potential tempera- 
ture and we regard it as the temperature a parcel will have if it is 
compressed or expanded adiabatically from any state (T,p) to the 
1000 mb level. It follows that the potential temperature remains 
invariant during an adiabatic process. Thus, under adiabatic con- 
ditions, 9 can be used as a tracer of air motion. On timescales for 
which a parcel can be considered as adiabatic, constant values of 
9 track the motion of this parcel, which takes place on a surface 
of constant 9. From equation (4.33) it is clear that the distri- 
bution of 9 in the atmosphere depends on the distribution of T 
and p. In the atmosphere dp/dz dT/dz (just think that 5 km 
above the surface pressure decreases on the average from about 
1000 mb to 500 mb while the temperature decreases from 288 K to 
238 K). Therefore, surfaces of constant 9 tend to be like isobaric 
surfaces. Before, we saw that for an ascent (descent) the parcel’s 
temperature must decrease (increase) owing to the work done by 
the parcel (by the surroundings). With the definition of potential 
temperature we can extend this statement to read as follows: dur- 
ing an adiabatic ascent (descent) the parcel’s temperature must 
decrease (increase) but in such a proportion as to preserve the 
parcel’s potential temperature. Note that equation (4.32) indi- 
cates that for non-adiabatic processes the change in the potential 
temperature is a direct measure of the heat exchanged between 
the parcel and its environment. As a result a non-adiabatic par- 
cel will drift across potential temperature surfaces in proportion 
to the net amount of heat exchanged with its environment. 


Examples 

(4.1) One mole of a gas expands from a volume of 10 liters and a 
temperature of 300 K to (a) a volume of 14 liters and a temper- 
ature of 300 K. and (b) a volume of 14 liters and a temperature 
of 290 K. What is the work done by the gas on the environment 
in each case? 


(a) According to the definition of work 


W = 



nR*T 

V 


dV. 


Since the temperature remains the same the above 
equation gives 


W = nR*T In 2 = 839J. 
Vi 
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(b) In this case T does not remain constant. As such it 
cannot be taken outside the integral and the above 
equation cannot be applied. 


When T does not remain constant, in order to calculate the 
work done we must have the explicit function describing the 
path from i to / in the (p, V) diagram, not just i and /. An 
approximation to this would be to define a temperature T 
or a pressure p that satisfied the relationship 


"' /2 TdV - 


'v i 


V 


fV 2 


' V 1 


dV 

V 


or the relationship 


rV 2 2 

/ pdV = p I dV 

/vi Jv i 


Since in our case no clues are given about the function 
describing the path from * to / we are free to make (rea- 
sonable) assumptions. We can thus assume that i and / 
are connected by a straight line. From the data we can 
estimate that pi = 249 420 Pa and Pf = 172 200 Pa. From 
Figure 4.6 it is easy to see that the work done is the area of 
the trapezoid ifVfVi. Because of the equality of the shaded 
triangles this area is equal to the area of the rectangle 
defined by the lines p = 0, p = p, V = Vi, V = V f where 
p = pi +p 2 /2. Thus 


W = p(V f - V) = 843 J. 


Figure 4.6 

Diagram for example 4. 1 . 
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Figure 4.7 

Diagram for example 4.2. 


P 


Pi 


PI 



7= constant 


V; Vf 


V 


(4.2) Consider an initial, i, and a final, /, state of a gas on an isotherm 
with Vf > Vi . The gas goes from i to / via an isobaric expansion 
and subsequent isochoric cooling; sketch on a (p, V) diagram 
the complete transformation and find the work done. If instead 
the gas goes from i to / via an isochoric cooling followed by an 
isobaric expansion, will the work done be different from before? 

In the first case the gas goes from i to m (isobaric from V z to 
Vf) and then from m to f (isochoric from p t to Pf) (Figure 
4.7). The total work done from i — > f will be the sum of the 
work done from i -4 m and from m — > / 

r v f r v f 

W 1 = pdV + / pdV 
J Vi JVf 

= Pi [ dv =Pi(V f - Vi). 

JVi 

In the second case the gas goes from i to m! (isochoric from 
Pi to Pf) and then from m' to / (isobaric from V to Vf). As 
before we have that the total work from i — > f is 

r Vi f v f 

W 2 = I pdV + / pdV 

JVi JVi 

= p f [ 1 dV=pf(Vf -Vi). 

JVi 

Since p t ^ p f it follows that W X ^W 2 - 

(4.3) Calculate Q , W, and A U for the following processes 

(Figure 4.8); 

(a) isothermal reversible compression from state i = (p i: Vi) 
to state / = (p f , V f ); 


EXAMPLES 


47 


Figure 4.8 

Diagram describing the 
processes in example 4.3. 



(b) adiabatic reversible compression from state i = (p i: V,j to 
state to = (pf,V m ) and a subsequent isobaric reversible 
compression to state / = (pf, V}); 

(c) reversible isochoric increase of the temperature from 
state i = ( Pi, Vi ) to state to' = (pf,Vi) and a subse- 
quent reversible isobaric decrease in temperature to state 
f = (Pf,V f ). 

Express all answers in terms of Pi,Pf, and T, where T is the 

temperature of the isotherm in (a). 

(a) During this transformation T t = T f = T , AT = 0 and 
A U = 0. Then from the first law it follows that 

Q=W= [ f pdV = nR*T 

J i 

= nR*T In 1 = nR*T In — . 

Vi p f 

(b) This transformation consists of two branches (i — > m 
and to — > /). From (a), the initial and final temper- 
ature is the same, i.e. At/ = 0. Thus, again we have 
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Q = W with 


At/ — A t/i_> m + AU m ^ f 

Q Qi—¥m T Qm—if 


w = W^ m + w m ^ f . 

Since the branch i — > m is adiabatic it follows that 
Qi^m = 0 and 

W^ m = -A U^ m = -C v (T m - T 4 ). (4.34) 

T m and T) are related via the equation 

r iP r^ = T mP tf, 

where p m = Pf- It follows that 


T m = T i [ — \ =T 


Then from equation (4.34) we obtain 


W^ m = -C V T -1 . (4.35) 

Branch m — > f is isobaric. Thus, 

w m ^ f = f pdV = p f (V f -V m ). (4.36) 

J m 

V rn and Vi are related via Poisson’s equation as follows 
ViV?=p m Vl- 

Since p m = Pf the above equation becomes 

(YnY = ?L 

\vj Pf 


v = V [ — 

v m v i \ 


Then equation (4.36) can be written as 


Wm-t-f = Pf V f ~V, [ — 
\Pf 


= PfVf 1 - w" (— 

Vf \Pf 
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= nR*T 


= nR*T 


1 _El(El 

Pi \Pt 

i- 


l/'Y 


1 - 


(4.37) 


From equations (4.35) and (4.37) it follows that overall 

\ / n\ 

Q = W = — C V T 


* ) -1 


+ nR*T 


= CAT 


1 - 


Pf 


i-l Bl 

Pf 

1 Y 

Pi ' 


Pf 


(c) Here also the transformation consists of two branches, 
i — > m! and m’ — > f . Branch i — > m' is isochoric, so 

Wi^ m f = 0 . 

Branch m' — > f is isobaric, which gives 

= [ f pdv = p f (V f - v m .) = Pf(V f - HO 

J m' 


W m ^ f =p f \ ) 1 -- 
v f 


= nR*T 1 - 


or 

K 

r i 

Pl 
Pi 

Here again T, t = Tf = T, AT = 0 and A U = 0. 
Thus, overall 

Pf 


Q = W = nR*T 1 - ^ . 

Pi. 


(4.4) This following problem is due to Iribarne and Godson (1973). 
Students should try to solve it before looking at the solution 
here. Figure 4.9 shows an insulated box with two compartments 
A and B. Both compartments contain a monatomic ideal gas and 
they are separated by a wall that on one hand does not allow 
any heat through it but on the other hand is flexible enough to 
ensure equalization of pressure in both compartments.The initial 
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Figure 4.9 

Diagram for example 4.4. 



conditions for both compartments are the same: T t = 273 K, 
Vi = 1000cm 3 , and ^ = I atmosphere (1013 x l0 2 Pa). Then, 
by means of an electrical resistance, heat is supplied to the gas 
in A until its pressure becomes ten times its initial pressure. 
Estimate: (a) the final temperature of the gas in B, (b) the work 
performed on the gas in B, (c) the final temperature of the gas 
in A, and (d) the heat absorbed by the gas in A. 


(a) The heat that is supplied to gas in A results in an 
increase in all its state variables. As it expands, it com- 
presses the gas in B but it does not provide it with any 
of its heat because of the insulating wall between them. 
Since the wall is flexible the gas in B is compressed 
adiabatically until its pressure equals that inside com- 
partment A. As such p B f = 10 atmospheres. From 
Poisson’s equation we can then solve for T B f 

TmPsl =T Bf p B } 

Since the gas is monatomic 7 = 1.666. Thus, 

T B /=(—) T Bi ss 686 K. 

\PBfJ 

Having T Bf we can use the other form of Poisson’s 
equation to find V B f\ 

TsiV^ 1 = TafVSj 1 


or 


V B f = 



0.25 V Bi 


1/7-1 

V Bi = 
= 250 cm 3 . 



V B 


1.5 
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(b) The work performed on the gas in B is given by 

W = -A U = —C v (T Bf - T m ) = -413 C y J. 

As we know, for a monatomic gas C v = \nR* . The 
number of moles n can be found from the ideal gas law 
and the initial conditions. It follows that 

n = ^ * = 0.0446 mol. 

R*T t 

Accordingly, C v = 0.556 JK -1 and 
W ft -230 J 


or 

W ft -55 cal. 


(c) 


Since the final volume in B is 250 cm 3 it follows that 
the final volume in A is 1750 cm 3 . By applying the ideal 
gas law for the initial and final conditions in A we find 
that 


Ta S = 


PAfV A f ^ 

T r -L Ai 
PAi VAi 


= 4777 K. 


(d) The heat absorbed by the gas in A is given from the 
first law 

Q = AU + W. 

Here W = 55 cal (i.e. the opposite of the work done on 
gas B), and A U = C v (T Af - T Ai ) = 2504 J = 598 cal. 
Thus 

Q = 598 + 55 = 653 cal. 


Problems 

(4.1) Commercial aircraft fly near 200 mb where typically the 
outside temperature is — 60 °C. (a) Calculate the tempera- 
ture of air if compressed adiabatically to the cabin pressure 
of 1000 mb. (b) How much specific heat must be added 
or removed (isobarically) to maintain the cabin at 25 °C? 
Consider the air as dry air. (337.5 K, 9.5 cal g _1 has to be 
removed) 

(4.2) A sample of 100 g of dry air has an initial temperature of 
270 K and pressure 900 mb. During an isobaric process heat 
is added and the volume expands by 20% of its initial vol- 
ume. Estimate: (a) the final temperature of the air,(b) the 
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amount of heat added, and (c) the work done against the 
environment. (324 K, 5427 J, 1550 J) 

(4.3) An air parcel moves from p. t to Pf. If its initial temperature 
is T t find: (a) the specific work done by or on the parcel 
if its change happens adiabatically, (b) the specific work if 
its change happens isothermally, and (c) the change in the 
parcel’s potential temperature in (a) and (b). 

(4.4) One mole of dry air has an initial state T = 273 K and 
p = 1 atmosphere. It undergoes a process in which its 
volume becomes four times its initial volume at 400 mb. If 
air is considered an ideal gas and if the process obeys the 
law pV v = constant, estimate (a) the value of rp (b) the 
final temperature, (c) the change in internal energy, the work 
done, and the heat exchanged between the air and its sur- 
roundings. Look at your results carefully and elaborate on 
what kind of process would generate results like these. Is 
this process a realistic process? (0.67, 431 K, 3286 J, 3989 J, 
7275 J) 

(4.5) Assume that you are on the top of a mountain at an altitude 
of 15 000 feet and there are no clouds above or below. If the 
temperature is —12 °C what will the temperature be at 3500 
feet altitude? (22.4 °C) 

(4.6) A parcel of dry air has a volume 10 liters, a temperature 
27 °C, and a pressure 1 atmosphere. The parcel (a) is com- 
pressed isothermally until its volume becomes 2 liters, and 
(b) expands adiabatically until its volume becomes 10 liters. 
Describe these changes graphically in a (p. V ) diagram. On 
it mark the values of p, V, T for each case. 

(4.7) A sample of dry air has a temperature of 300 K, a volume 
of 3 liters, and a pressure of 4 atmospheres. The air then 
undergoes the following changes: (a) it is warmed under a 
constant pressure to 500 K, (b) it is cooled under constant 
volume to 250 K, (c) it is cooled under constant pressure 
to 150 K, and (d) it is warmed under constant volume to 
300 K. (1) Describe graphically on a (p,V) diagram each of 
the changes and mark on it the final values for pressure and 
volume after each change. (2) Calculate the total work done. 
(405.2 J) 

(4.8) For an ideal gas what on a (p, V) diagram describes the lines 
for which U = constant? 

(4.9) An ideal gas undergoes two transformations from state 1 to 
state 2 which follow the two curves in Figure 4.10. In which 
of the transformations is (a) the change in internal energy 
greater? (b) the heat absorbed greater? 

(4.10) One kilogram of dry air is warmed from 20 °C to 70 °C under 
a constant pressure of 1 atmosphere. Calculate: (a) the heat 
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Figure 4. 1 0 

Diagram for Problem 4.9. 



V 


absorbed by the gas, (b) the work done, and (c) the change 
in internal energy. (50 250 J, 14 350 J, 35 900 J) 

(4.11) Hydrogen under p = 5 atmospheres and T = 20 °C is 
warmed isobarically until its volume increases from 1 to 2 
liters. Calculate the absorbed heat. (1773.4 J) 

(4.12) Inside a cylinder closed by a piston there is 88 g of an 
unknown diatomic gas of a temperature 0 °C. The gas is then 
compressed adiabatically until its volume becomes equal to 
1/10 of the initial volume. The change in internal energy is 
17158 J. Identify the gas. (Carbon dioxide) 

(4.13) Provide an example that proves that heat transfer is an 
irreversible process. 

(4.14) At T = 0 °C and p = 1000 mb, 1 g of dry air receives during 
an isochoric process an amount of heat. It is then observed 
that its pressure increases by 50 mb. What is the change in 
the temperature of the air and what is the amount of heat 
absorbed? (13.6 °C, 2.33 cal) 

(4.15) Calculate the change in specific energy of an adiabatic 
parcel whose speed changes from 10ms _1 to 25 ms -1 . 
(-262.5 J kg -1 ) 



CHAPTER FIVE 


The second law of 
thermodynamics 


The first law of thermodynamics arose from the conservation of 
energy principle. The first law, even though it implies that we can- 
not create or destroy energy, places no limits on how energy can be 
transformed from one form to another. Thus, on the basis of the 
first law, heat can be transformed into work, work into heat, work 
can be done at the expense of internal energy, and so on. However, 
if no other laws existed the first law would allow certain phenom- 
ena to happen that never happen in reality. For example, consider a 
heavy body falling on the ground. We will observe that during the 
impact the body will warm. The opposite phenomenon according 
to which a body at rest on the ground begins to rise by itself while 
it is cooling is impossible. Similarly, no engine has yet been built 
which, for instance, would receive heat from the sea, transform it 
to work, and then set a ship in motion. Both the above examples 
are not in disagreement with the first law since the work would be 
done at the expense of the internal energy of the soil or the sea. 
The impossibility of these phenomena is due to the second law of 
thermodynamics, often hailed as the supreme law of nature. We 
will start our discussion of this law with the following example. 


5.1 The Carnot cycle 

The Carnot cycle is a thermal engine. A thermal engine is one 
that receives from some source an amount of heat, part of which it 
transforms into work. In thermodynamical terms the Carnot cycle 
is a cyclic transformation consisting of the following four steps 
(Figure 5.1): 

Step 1: reversible isothermal expansion (1 — ► 2), T = r I\ = 
constant 

Step 2: reversible adiabatic expansion (2 — »• 3), 9 = 9 1 = constant 
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Figure 5.1 

Steps in a Carnot cycle. 



Step 3: reversible isothermal compression (3 — > 4), T = T 2 = 
constant 

Step 4: reversible adiabatic compression (4 — > 1), 0 = d 2 = 
constant 

with T 2 <T l , 6 2 < 9 1 . 

Let us now calculate the work done and heat absorbed when an 
ideal gas is subjected to this cycle. The change in internal energy 
is zero because the overall transformation is a cyclic one. We will 
consider each step separately. 

Step I 

AT = 0 

r V2 r V2 dv v 2 

W 12 = pdV = nR*T \ / — = nR*T x In -± 

Jv ! Jv 1 V bi 

A U 12 = C V AT = 0 
Q 12 = bFi2 

Since V 2 > V 2 it follows that W 12 ,Q 12 > 0 indicating that during 
step 1 work was done by the gas and heat was absorbed by the gas 
from the heat source at temperature Ti. 

Step 2 

Q 23 = 0 

A U 23 = C v {T 2 - TO < 0 
W 23 = -A U 23 = -C v (T 2 - TO > 0 


AT = 0 



nR*T 2 



Step 3 
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At/34 = 0 

Q34 = W 34 

Since V 4 < V 3 it follows that W 34 , Q 34 < 0 indicating that now work 
is done on the gas and heat is given away by the gas to the source 
at temperature T 2 . 

Step 4 

Q 41 = 0 

A U 41 = Cy(T ! - T 2 ) > 0 

W 41 = -At /41 = -C v (T 4 - T 2 ) < 0 

From the above it follows that the total work done is 

W = W 12 + W 23 + W 34 + W 41 

= nR*T ! In ^ + nR*T 2 In ^ (5.1) 

Vi V 3 

Because transformations 2 — > 3 and 4 — > 1 are adiabatic the 
following relations apply: 

TxV ^ 1 = T 2 Vr X 


and 

T x vr x = T 2 V2~ x . 

By dividing the above two equations we have 

= V 3 

v 4 v 4 - 

Accordingly, equation (5.1) can be written as 

V 2 V 2 

W = nR*T 4 In — - nR*T 2 In — 

V) V) 

or 

W = nR * ^ln ^ (T, - T 2 ). 

The total amount of heat absorbed during this cycle is 

Q = Q 12 + Q34 = ll'd 2 + W'34 


(5.2) 


or 

Q = W = nR* (hi jQ (7\ - T 2 ) (5.3) 

For simplicity we will now denote the amount of heat absorbed by 
the gas from the heat source at T) (i.e. Qi 2 ) as Q 4 and the amount 
of heat given by the gas to the source at T 2 (i.e. Q 34 ) as Q 2 . Since 
W is the area enclosed by the cycle it follows that W > 0. Then 
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from equation (5.3) it follows that Q > 0. Since W = Q 1 + Q 2 and 
Q 2 is negative it follows that only part of the heat absorbed by 
the gas at the source T) (higher temperature) can be transformed 
into work. The other part is surrendered to the source at T 2 (lower 
temperature). We define the efficiency of the Carnot cycle, r/, as 
the ratio between the work done and the heat absorbed by the gas 
at the source T). 


Q i + Q 2 _ , Q 2 
Qi Q~i 


(Qi > 0 ,Q 2 < 0 ). 


By considering the relations derived for Q 1 and 
the above ratio as 


rj = 1 + 


nR*T 2 In £ 

V3 

In || 


Q 2 


we can write 


or 


77 = 1 — 


T 2 

T] 


(5.4) 


5.2 Lessons learned from the Carnot cycle 

• According to equation (5.4) the thermodynamic efficiency of the 
cycle depends only on the temperatures of the two heat sources 
and it becomes zero if T 2 = Ti (i.e. if the two sources collapse into 
one source). As such we conclude that doing work via a thermal 
engine operating at only one source of heat is impossible. This is 
known as Kelvin’s postulate and is an expression of the second 
law of thermodynamics. Another way of stating this postulate 
is that it is impossible to construct an engine which transforms 
heat into work without surrendering some heat to a source at a 
lower temperature. This explains why stones on the ground do 
not rise in the air or why no machine has been discovered which 
by taking heat from the sea or the air sets in motion a ship or 
a car. According to the second law, apart from the warm heat 
source (sea or air), a second colder heat source would be required 
which we do not really have. I use the word really because some 
differences in temperature do exist in both the oceans and the 
atmosphere. For example, the temperature of the surface air is 
around 35 K warmer than the air at 5 km, and the temperature 
at the surface of the oceans is about 20 K warmer than at a depth 
of 1 km. In both cases, however, it turns out that the Carnot effi- 
ciency is about 10%. This is a low efficiency given the tremendous 
expenses that would be required to build such engines. On the 
other hand if the atmosphere, in the horizontal, is thought of as 
a Carnot cycle operating at T) = equatorial temperatures and 
T 2 — polar temperatures we can explain why in the winter the 
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circulation is stronger. Recall from above that 


r) = 1 + 


Q 2 
Qi 


= l - 


Ti 

Ti 


or 

Q i + Q 2 _ 21 — t 2 

Qi “ Ti 

If we think of summer as the limiting case where Tj — T 2 — F 0, 
then during summer Q 2 — F — Qi and 77 — F 0. On the other hand if 
we think of winter as the limiting case where T 2 << Ti, then 
during winter ij — F 1 and Q 2 — F 0. It follows that the effi- 
ciency of this “machine” is higher in the winter than in the 
summer. Accordingly, the amount of heat given by the atmo- 
sphere to the cold source is smaller in the winter, than in the 
summer. Since the amount of heat absorbed by the atmosphere 
from the warm source is more or less constant, this means that 
there is more heat “remaining” in the atmosphere in the win- 
ter than in the summer that can be transformed into kinetic 
energy. 

Apart from the impossibility of transformations whose only final 
result is to transform into work heat extracted from only one 
source of heat, there is another type of transformation that is 
impossible. Such transformations refer to transfer of heat from 
a cold body to a warm body without work being done. This 
phenomenon is allowed by the first law (since heat lost by the 
cold body will be exactly equal to the heat gained by the warm 
body) but has never been observed. Of course transformations 
that result in heat flow from a colder body to a warmer body 
are allowed if work is done (such as in refrigerators). This leads 
to another expression of the second law of thermodynamics. A 
transformation whose only final result is to transfer heat from a 
body at a given temperature to a body at a higher temperature 
is impossible. This is known as Clausius’s postulate. It is easy 
to show that both postulates are equivalent and two expressions 
of the same law. To do that we shall prove that if Clausius was 
wrong then Kelvin would be wrong and vice versa. If Kelvin was 
wrong and we could transform into work heat extracted from only 
one heat source, then we could transform this work (by means 
of friction) into heat and with this heat raise the temperature 
of a given body of a higher initial temperature. This would be 
a violation of Clausius’s postulate. Now if Clausius was wrong 
and we could transfer a certain amount Q 1 from a source at 
T 2 to a source at Ti (T 2 < Ti) in such a way that no other 
change in the state of the system occurred, then with the help 
of a Carnot cycle we could absorb this amount and produce an 
amount of work, W. Since the source T x receives and gives up 
the same amount of heat (Qi) it suffers no final change. But 
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what we have just described is a process whose only final result 
is to transform into work the heat extracted from a source at the 
same temperature T 2 throughout. This is contrary to Kelvin’s 
postulate. 

• From the definition of efficiency it follows that 


Tl 

T\ 


-i , Q 2 1 t 2 

= -tt (Qi>0,Q 2 <0). 
bn 


Based on the above relationship an important application of 
Carnot cycle is that it offers a definition of the absolute tempera- 
ture (Kelvin) scale based on pure thermodynamic arguments. For 
this we only need to have a value for a given temperature. An 
example for such value could be the melting point of ice. As such 
if we choose for the melting point of ice the value Tj = 273 °C 
every other temperature T 2 can be defined via a Carnot cycle 
operating between heat source Tj (consisting of melting ice) and 
a heat source T 2 . If we then measure Q 1 and Q 2 we can apply 
the above relationship to calculate T 2 . This thermodynamic def- 
inition of the Kelvin scale allows us to present a thermodynamic 
definition of absolute zero: the absolute zero is the temperature 
of a heat source (assumed to be the cold source of a Carnot 
cycle) which absorbs no heat (i.e. Q 2 = 0) even though the cycle 
produces work. 

• In reality natural processes are irreversible processes. However, 
if some processes happen very slowly and the heat losses are 
small enough not to affect the temperature (which supposes that 
friction is also negligible), then the results will not differ from 
those expected from perfectly reversible processes. Such processes 
are often called partially reversible. If on the other hand fric- 
tion and heat flow are important then the degree of reversibility 
decreases. In nature the degree of reversibility varies between 
almost perfectly reversible and perfectly irreversible. 

If in a thermal engine we increase the degree of irreversibility 
(for example, by increasing friction or direct leakages of heat), 
then the work from a given amount of heat will decrease. If 
the thermal losses keep on increasing we will end up with a 
perfectly irreversible engine. It follows that the efficiency of a 
thermal engine can take on values from zero (corresponding to 
the worst thermal engine operating at perfectly irreversible condi- 
tions) to some maximum value (corresponding to the best engine 
operating at perfectly reversible conditions). The Carnot cycle 
is a perfectly reversible cycle and as such it gives the maximum 
value the efficiency can take. This expresses the Carnot theorem 
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which states that it is impossible to construct an engine oper- 
ating between two heat sources which would have an efficiency 
greater than the efficiency of a Carnot cycle operating between 
the same heat sources. However, this maximum value will not 
be equal to one as this would require that R = OK which is 
practically impossible. Recall from section 5.1 that 


Qi = nR*T i In 


C 2 

n 


and 


It follows that 


Q 2 = nR*T 2 In 

*3 


Q i 
R 


= nR* In 


V 2 

R 


and 


Q 2 p* , V 2 
- = -nR In - 


or 

^ + ^ = 0 (Q!>0,Q 2 <0)- (5.5) 

i 1 ±2 

Now by considering equation (5.4) we can, for the reversible 
Carnot cycle, rewrite equation (5.5) as 


Q 1 ^/rev Q 2 

T x ~T 2 + R 


Since Q 2 /T 2 is less than zero it follows that Qii] iev /(Ti — T 2 ) 
must be greater than zero. Since their sum must equal zero 
and, as we discussed above, the value of is maximum for per- 
fectly reversible processes, it follows that for irreversible processes 
where rj becomes smaller 

Q 1 ^7irrev t Q 2 ^ 

R~T 2 R ■ 

This proves that 



where the equality applies to reversible processes. The above 
argument can be extended to any cyclic transformation in which 
the system receives heat from or surrenders heat to a set of tem- 
peratures R,R, . . . ,T n (Fermi, 1936). If the amounts of heat 
exchanged are Qi, Q 2) • • • , Qn (positive when heat is given to 
the system, negative otherwise) it can be shown that 
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where again the equality applies to reversible processes. In the 
limit where TV — > oo the above equation takes the form 


/ 


SQ 

T 


< 0 . 


(5.6) 


We may conclude that with the aid of the Carnot cycle we can 
introduce a quantity, which for reversible processes is an exact 
differential. This quantity, 5Q/T, is defined as dS, where S is a 
new state function called entropy. Since 5Q/T is, for reversible 
processes, an exact differential it follows that for any reversible 
transformation i — > f 

f6 ^=A S=S / -S i . (5.7) 


In other words the change in entropy depends only on the initial 
and final states, not on the particular transformation. This result 
is fundamental in the formulation of the second law. 

Since S' is a state function it follows that entropy changes are 
due to changes in both temperature and volume (or pressure). 
Using the first law we have that for a reversible process 


or 


or 


or 


or 


SQ dT dV 
dS = — = C v — + p— 


dT 


dV 


dS = C v —— + nR* — — 


T 


V 


dS _ dT nR * dV 

c y ~ ~Y + c7~u 


dS 
c v 


dT 

~T 


+ (7- 1) 


dV 

U 


Sf = Si + Cy In 


( Tfvr 1 ) 


(5.8) 


Equation (5.8) provides the change in entropy of an ideal gas 
(since the ideal gas law was used) as a function of the initial and 
final temperature and volume. This relation cannot, therefore, 
be used for liquids or solids. For liquids or solids C v = C p = C 
and one can show (but we will not go into this proof here) that 
in this case the entropy change is dominated by the temperature 
change. As such for liquids or solids and for reversible processes 

S f -Si = C ln^f 


(5.9) 



5.3 MORE ON ENTROPY 


63 


5.3 More on entropy 

As in the case of energy, we will choose an arbritary equilibrium 
state O as the standard state where by definition So is zero. Then 
we will define the entropy S A of an equilibrium state A as 


where the integral is taken over a reversible transformation. It can 
easily be shown that if instead of a standard state O we had chosen 
a standard state O' then the entropy of state A would differ from 
the original one by an additive constant: from the above equation 
we have 


Since O' is fixed, S 0 ' is some constant. The entropy is thus defined 
except for an additive constant. As in the case of energy this does 
not present a problem as long as we are dealing with differences, 
not actual values of entropy. Now let us consider the cyclic trans- 


versible transformation. As we know, for any cycle j> SQ/T < 0. As 
such 




From equation (5.7) it follows that 

S' A = S A - So . 


or 



formation shown in Figure 5.2. The path i — > f represents a 
reversible transformation and the path / i represents an irre- 



or 



Figure 5.2 


f 


A cyclic transformation 
consisting of a reversible 
part from i to f(R) and 
an irreversible part from 

/ to i(I). 


R 


I 
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Since i 


f is reversible we have that 
rf SQ 


s f ~s t = 


T 


Combining the last two equations leads to 

r SQ 


S f - & 


T 


< 0 


or 


SQ 

T 


< St - S f . 


By interchanging i and / we conclude that for any process 

" S 4<s,- Sl 


or 


or 


T 


as > !‘4 


T 


- T 


(5.10) 

(5.11) 

(5.12) 


where the equality applies to reversible processes only. Equation 
(5.12) is the general expression of the second law of thermody- 
namics and it indicates that the upper bound to the heat that 
can be absorbed by the system during a given change is SQ = 
TdS. For a completely isolated system (SQ = 0) equation (5.10) 
translates to 

S f > Si. (5.13) 


This is an important conclusion as it indicates that for any spon- 
taneous irreversible transformation (i.e. one not related to external 
influences) occurring in an isolated system, the final entropy is 
greater than the initial entropy. It follows that when an isolated 
system has attained a state of maximum entropy it cannot undergo 
any further transformation because any change would decrease 
its entropy, and this is not allowed by equation (5.13). Thus the 
state of maximum entropy is a state of stable equilibrium. Note 
that if the isolated system consists of a number of microscopic 
subsystems it would be possible for some of the subsystems to 
reduce their entropy but only at the expense of the rest of the 
subsystems, whose entropy must increase enough that the over- 
all entropy of the system increases. A perfect isolated system is 
our universe (assuming that no other universes exist with which 
our universe can interact). It follows that the entropy of our uni- 
verse increases in time. Equations (5.10)-(5.12) as well as equation 
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(5.13) are the basic mathematical expressions of the second law of 
thermodynamics . 


5.4 Special forms of the second law 

Given the above definitions of the second law we can derive its 
expression for special cases such as the following. 


• Finite isothermal transformations (A U = 0) 

From equation (5.11) it follows that in this case 


A S>- 


i r f 


6Q 


or 


or 


AS> 2 
AS>f. 


• Adiabatic transformations 

Using equation (5.12) we obtain 


(IS > 0. 


(5.14) 


(5.15) 


• Isentropic transformations 

An isentropic transformation is one during which the entropy 
does not change. In this case it is clear that (recall equation 
(5.12)) 

5Q < 0. (5.16) 

Note that according to equation (5.15) a reversible adiabatic 
process is isentropic. 

• Isochoric transformations 

From the first law we have that when dV = 0 then 6Q = C v dT. 
As such for isochoric transformations 

i S>c v f 

or 

AS’ > C v In . (5.17) 

• Isobaric transformations 

In this case 

SQ = C p dT 

and it follows that 

dT 
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or 

AS > C p In . (5.18) 

i 

From equation (5.14) it follows that irreversible work increases 
a system’s entropy. From equation (5.17) it follows that in the 
absence of work the change in entropy depends on the relation 
between T f and 7). 


5.5 Combining the first and second laws 

Consider the following form of the first law: 

5Q = CpdT — Vdp. 

Combining the above equation with dS > yields 
TdS > C p dT - Vdp. 

Recalling that C p = dH/dT we can reduce the above equation to 



TdS >dH- Vdp 


or 

dH < TdS + Vdp. 

(5.19) 

Similarly, from 

SQ = C v dT + pdV 


we arrive at 

TdS >dU + 5W 


or 

dU < TdS - pdV. 

(5.20) 

It is often convenient to introduce two new functions: the Helmholtz 
function F = U — TS and the Gibbs function G = H — TS = U + 
pV — TS. Since S = S(T, V) and U = U (T) it follows that both 
F and G are state functions (i.e. F = F(T, V),G = G(T, V)), and 
thus are exact differentials. The advantage of these functions is that 
they can be used to express equations (5.19) and (5.20) in a form 
where the pairs ( T,p ) and (T,V) appear as the independent vari- 
ables instead of the pairs (S,p) and (5, V). In this case equations 
(5.20) and (5.19) can be written as: 


dF < - SdT - pdV 

(5.21) 

and 

dG < —SdT + Vdp. 

(5.22) 


The interpretation of the above functions is that for isothermal 
processes dF < — SQ or dF < — 5W which makes F the energy 
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available for conversion into work. The usefulness of G becomes 
clearer for isothermal-isobaric transformations. During such trans- 
formations (which apply in phase transitions, for example water to 
vapor) dG = 0 and therefore G is conserved. Relations (5.19)-(5.22) 
are often referred to as the fundamental relations. 


5.6 Some consequences of the second law 

Thermodynamic definition of temperature 

As we know, the integral ('/ 5Q depends on the path from i to /, 
not just the states i and /. In this chapter we showed that the 
integral // 6Q/T depends only on the initial and final states, not 
on the particular path from i to /. As such 1/T is the integration 
factor that makes SQ an exact differential. It follows that in ther- 
modynamical terms temperature is the inverse of the integration 
factor of the differential of heat in reversible processes. 


The statistical nature of thermodynamics 

The change in entropy during a reversible cooling at constant 
pressure of m grams of dry air from Tj to T 2 (Xj > T 2 ) is 


A S = 



SQ 
T ' 


Since p = constant we have that SQ = mc pd dT where c pd is the 
specific heat capacity of dry air. It follows that 


AS = 



dT 

niCpd ^ 


or assuming that c pd does not vary significantly with temperature 


T 2 

AS = St 2 - S Tl = mc pd In — . 

-L l 

Since m and c pd are greater than zero and In fff <0 it follows 
that AS < 0. Thus, the entropy decreases during the cooling. This 
decrease does not violate the second law. Even though the dry air 
decreases its entropy the environment that provided the heat suffers 
a positive change so that the overall entropy change of the system 
is positive. We can extend the above arguments to show that as 
the temperature goes down so does the entropy (or as temperature 
goes up so does the entropy). The question is why and what is the 
physical meaning of this. Kinetic theory provides an explanation 
for the above. According to the theory, at very low temperatures 
the motion of the molecules is very slow with the molecules being 
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Figure 5.3 

The situation in (a) 
represents greater order 
and a lower number of 
complexions than in (b) 
where the number of 
complexions and disorder 
have increased. 



e 

N, = N N 2 = 0 

• ' . A 

. ‘ • B 


N 1 = N/2 N 2 = N/2 


more or less uniformly distributed in space. This is a picture of a 
high order, with complete order (motionless molecules uniformly 
spaced) presumably achieved at absolute zero. As the temperature 
increases the motion of the molecules increases and the order is soon 
destroyed. If we combine this with the above result it follows that we 
can associate entropy with order or disorder. Accordingly, decreas- 
ing entropy implies increasing order and increasing entropy implies 
increasing disorder. This association justifies the word entropy, 
which in Greek means the “inner behavior” (eisTpon La) . Now if 
such an association is correct we should be able to find a measure 
that defines disorder which should be mathematically related to 
entropy. 

Let us consider a volume V consisting of two compartments A and 
B , with TV total particles moving around. Let us denote with TVj 
the number of particles in compartment A and with TV 2 the number 
of particles in compartment B. Statistics tells us that the number 
of different ways, P, in which TVi particles will be found in A and 
N 2 particles in B is 


TV! 

NANA-' 


(5.23) 


Now consider two cases as shown in Figure 5.3. In the first case all 
particles are found in A (Figure 5.3(a)), and in the second case the 
particles are all over the box (Figure 5.3(b)). The first case (N-i,= 
TV, TV 2 = 0) obviously represents a state of more order than the 
second where it can be assumed that TVj = TV 2 = TV/2. Applying 
equation (5.23) to both cases yields 


TV! 

tvToT 


= l 


Case 1: 
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Case 2: 


N/ 2! N/ 2! ' 

We thus see that we have a measure P (often called the number 
of complexions) that increases as disorder increases. If we consider 
the first case as the initial state and the second case as the final 
state we can conclude (since i — > f is an irreversible process) that 
as the entropy increases P increases. So P is a possible candidate 
to relate disorder to entropy. But how can we find their functional 
relationship? 

The derivation of this relationship (one of the most beautiful yet 
simplest ones) is due to Boltzmann (Fermi, 1936). From the above 
we saw that as S increases, P increases and vice versa. So let us 
start with a general relationship of the form 

S=f(P). 

Next, consider a system consisting of two subsystems and let Si 
and S 2 be the entropies and Pi and P 2 the corresponding number 
of complexions. Then, Si = /(Pi) and S 2 = /(P 2 ). For the whole 
system we can then write that 

S = Si + s 2 


and 


P = P 1 P 2 


or 

/(P 1 P 2 ) = /(Pi) + /(P 2 ) 

i.e. the function / obeys the functional equation 
f(xy) = f(x) + f(y). 

Since the above equation is true for all values of x and y we may 
take y = 1 + e with eCl. Then we can write it as 

f(x + xe) = f(x) + f(l + e). (5.24) 

Expanding both sides of the above equation using Taylor’s theorem 
and neglecting terms higher than first order we get 

f(x) + xef'(x) = f(x) + /( 1) + e/'(l) (5.25) 

where /' indicates the first derivative. For e = 0 from equation 
(5.24) we have that /( 1) = 0. Therefore, equation (5.25) reduces to 

xef{x) = e/'(l) 


or 


Xf Or) = /'( 1). 
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Since /'( 1) is a derivative evaluated at a value of 1 it is a constant. 
Thus 


xf'(x) = k 


or 


or 


or 


/'M = k - 

X 

f(x) = j * dx 
f{x) = k In x + constant 


or by changing x — > P, 


S = k In P + constant. (5.26) 


This equation relates order and entropy and constitutes a funda- 
mental relationship connecting statistical mechanics and thermo- 
dynamics, thereby giving thermodynamics a statistical character. 


Entropy and potential temperature 

Recall that the potential temperature 6 is given by 

q ^ rp ^1000^ fl/Cp 

By taking the logarithms of both sides we obtain 

D JD 

In 6 = In T 4 In 1000 In p 


or 

c p d In 6 = c p d In T — Rd In p. 
From the first law we have 

5Q = CpdT — Vdp 


(5.27) 


or 


or 


6Q _ dT dp 

rj~i ^ P rj-t * rj\ 


$Q n dT dp 

—— = Cp — — mK — 

T p T p 


or (for reversible processes) 


dT dp 

- o — 

T p 


(5.28) 
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Figure 5.4 

A hypothetical upper-level 
flow and the trajectory 
followed by a parcel. 



where ds is the specific entropy. Combining (5.28) and (5.27) 
leads to 


ds = c p d In 6 


(5.29) 


or 


s = c p In 9 + constant. 


(5.30) 


Thus, except for an additive constant, the specific entropy of a sys- 
tem is given by the logarithm of the potential temperature. When 
9 remains constant the entropy remains constant. It follows that a 
reversible adiabatic process is isentropic. For an irreversible adia- 
batic process ds > 0 ,d0 = 0. In this case, the increase in entropy 
comes from irreversible work (for example, dissipation of kinetic 
energy due to friction). It follows that an isentropic process is 
adiabatic, but an adiabatic process may not be isentropic. 

Atmospheric motions 

Figure 5.4 shows a hypothetical upper-level flow close to some lat- 
itude A. Along this latitude the radiative-equilibrium temperature 
T R e(A) is assumed to be the same. Such thermal structure can 
be achieved if the motion is completely zonal (i.e. parallel to lati- 
tude circles), since air parcels then have infinite time to adjust to 
local thermal equilibrium. If at position 1 a parcel is allowed to 
move with the flow two possibilities exist. If the motion is very 
slow the parcel has time to come into equilibrium with its sur- 
roundings, so its temperature differs from T RE only infinitesimally. 
Therefore, we can assume that heat transfer along the trajectory 
of the parcel is reversible. Then for positions 1 and 2 we have 
that 






2 


2 
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If we assume that 5q depends only on the parcel’s temperature (for 
example, Sq = Tf(T)) then we can write the above equation as 

c p d In 9 = J df(T) = /(T 2 ) - /(TO. 

Since T 2 = T) = T RE the above relationship reduces to 

c p d In 0 = 0 
or 

9 2 = Ox- 

Thus, the parcel is restored to its initial thermodynamic state when 
it returns to the original latitude. Since 9 depends on p and p 
depends on height it follows that under the slow motion condi- 
tions, successive crossings of the initial latitude result in zero net 
vertical motion and zero heat transfer. 

On the other hand if the motion is fast enough that reversibility 
cannot be assumed, then 

Cpd In 9 < ds ^ 0 

and therefore 9 2 ^ 9\. In this case when the parcel returns to the 
initial latitude it is displaced vertically and the net heat trans- 
fer is not zero. Vertical motion resulting from irreversibility in 
the atmosphere plays an important role in maintaining meridional 
circulations, heat and moisture transfer, etc. 




Examples 

(5.1) Imagine a volume V divided into two smaller volumes V) and 
V 2 by a partition. Assume that m i grams of an ideal gas at 
temperature T) occupy volume V) and m 2 grams of another 
ideal gas at temperature T 2 occupy volume V 2 . Then the partition 
is removed. If the specific heat capacities at constant volume of 
the gases are c Vl and c V2 , find (a) the limiting temperature of the 
mixture, and (b) the change in entropy for each gas and for the 
system as a whole. 

(a) We assume that while the two gases exchange heat 
between them and perform work, the system as a whole 
does not exchange heat with its surroundings and its 
volume remains constant (V = Vi + V 2 ). In other words 
for the whole system the process of mixing is adiabatic, 
with no work being done. As such, AU = 0. The total 
internal energy before mixing is C Vl T x + C V2 T 2 and 
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after the mixing it is C Vl T + Cy 2 T (here we have made 
use of the approximation U ~ C V T+ constant and 
then neglected the constant). Since A U = 0 it follows 
that 

C Vl T\ + C \/ 2 7 2 — Cy x T + Cy 2 T 
which yields 


T = 


a 


Vi 


c Vl + a 


-Ti 


v 2 


a 


v 2 


c Vl + a 


V 2 


(5.31) 


Considering that C Vl = c Vl nii and Cy 2 = Cy 2 m 2 the 
above equation becomes 

ji _ c vi m i j i , Cy 2 m 2 

c Vl m i + Cy 2 m 2 Cy^mx + Cy 2 m 2 


(b) From equation (5.8) and the additive property of 
entropy it follows that for each gas the total entropy 
change is the sum of the entropy change due to 
the temperature change of the gas alone and the 
entropy change due to the volume change of the gas 
alone, i.e. 


AS 1 ! = A S 1T + A S 1V 
AS 2 — AS 2 t AS 2 y 


where the lower limits for AS 1T , AS 1Vl AS 2T , and 
A S 2V are given by (recall equation 5.8) 


ASit — Cy 1 In — 

^ i 

T 
T 2 


AS 2 t — Cy 2 In 


A5 1V = C Vl In 


Vi + V 2 
Vi 


71-1 


= rriiRi In 


Vi + V 2 
V 1 


and 


A S 2 y — Cy 2 In 


( V1 + V2 

V v 2 


72-1 


m 2 R 2 


In 


Id + V 2 

v 2 


where R x and R 2 are the specific gas constants of the 
two gases. For the system as a whole the lower limit 
for the total entropy change is 


AS = A S 1 + AS 2 


or 


AS= (AS 1T + AS 2T ) + {AS 1V + AS 2V ) = AS T + ASy 
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or 


A S = 


c Vl m i 



+ c V2 m 2 In 


T' 

T 2 _ 


(5.32) 


m 1 R 1 


In 


V! + V 2 
Vi 


+ m 2 R 2 In 


V x +V 2 

^2 


In equation (5.32) the second term on the right-hand 
side (ASV) is necessarily greater than zero. The first 
term (A S T ) is either positive or zero when Tj = T 2 (see 
problem 5.13). As such for the whole system AS > 0 
which according to (5.12) indicates irreversible process. 
This proves that mixing is irreversible. 


(5.2) Calculate the change in air pressure if the specific entropy 
decreases by 0.05 J g 1 Kr 1 and the air temperature decreases 
by 5%. 

Assuming that the process is reversible we have 


dT 


dS = C p — - —dp 


V 

T C 


= c„A_ nfi .* 

p T p 


Since we do not have any information on the mass of the 
air involved in the process we will write the above per unit 
mass: 


or 


ds = 


dp 

P 


c 


dT dp 

P rp H 

T p 

c p dT ds 
R^T ~ Id 


or considering the air as dry 


In Pf 1005 /0-95rA (-50) 

Pi 287 \ T.i ) 287 


or 

Pf = 0.994 Pi 


The pressure will decrease by a small amount of about 0.6%. 


(5.3) In example 4.3, what is the change in entropy in each of the 
three transformations? 

(a) The first transformation is an isothermal reversible 
compression. Thus, 
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A S = 




Q = nR* In = nR * In — . 

T Vi p f 


(b) The second transformation consists of two branches. 
The first one is a reversible adiabatic transforma- 
tion. Thus, the entropy change in the first branch is 
zero. The second branch is an isobaric reversible com- 
pression. As such the entropy change of the second 
transformation is (recall that since and Tf lie on 
the same isothermal T t = Tf = T) 


A S = 


rf SQ_ 

rr i P 


rf dT_ T f _ T, 
C p In C p In 

m m. m. 


= C„ In [ ^ 
Pi 


= - nR * In ^ 


Pi 


= nR* In — 
Pf 


(c) The third transformation also consists of two branches: 
one reversible isochoric and one reversible isobaric. The 
total change in entropy is the sum of the changes in 
entropy in the two branches 


AS = A Si + AS 2 



SQ 

T 



a 


dT 


= C v In 


T / 
T 

1 l 


+ C p In 


T / 


During this transformation PiVi = nR*T.i and PfV t = 
nR*T rn t. Thus 

T — ELt 
Pi 


If we substitute the above value of T m / in the previous 
equation we obtain 


AS = C v In — + In — . 

Pi U p f 


Since Tf = T t = T it follows that 

AS = C v In ^ + C_ In ^ 
Pi Pf 
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or 

AS = nR* In — . 

Pf 

Thus, in all transformations the change in entropy is the 
same. 


Problems 

(5.1) Starting from the first law show that SQ/T is an exact 
differential. 

(5.2) Derive the so-called Maxwell relations for reversible pro- 
cesses 



(5.3) Find the lower limit for the change in specific entropy of dry 
air at a given temperature and pressure which is allowed to 
expand freely into an insulated chamber to assume twice its 
original volume. (200 J K -1 kg” 1 ) 

(5.4) What is the form of isenthalpic lines in a (p, V) diagram? 

(5.5) Dry air at T= 0 °C and p= 1 atmosphere is compressed isen- 
tropically until its pressure becomes 10 atmospheres. Find 
the final temperature. (527.5 K) 

(5.6) In a Carnot cycle the area for each transformation corre- 
sponds to the following work (a) isothermal expansion: 
31 165 J, (b) adiabatic expansion: 21 517 J, (c) isothermal 
compression: 24 282 J, (d) adiabatic compression: 21 517 J. 
Calculate (1) the heat provided, per cycle, to the gas by 
the warm source, and (2) the efficiency of the cycle. (6883 J, 
0 . 22 ) 

(5.7) One mole of ideal gas occupies 5 liters at 273 K. If it expands 
into a vacuum until its volume becomes 20 liters, what is the 
lower limit for the change in entropy and the upper limit for 
the change in Gibb’s function? (11.5 JK -1 , — 3146 J) 

(5.8) A fixed volume 2V is divided into two equal volumes V 
separated by a diaphragm. Each volume contains the same 
amount of the same gas at the same temperature T. Using 
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Boltzmann’s expression for entropy, show that after the 
diaphragm is removed the entropy change is zero. How 
will you approach this problem if the temperature in the 
two volumes is different? Use the approximation In N\ = 
N In N - N. 

(5.9) In a thermodynamic process a parcel of dry air is lifted 
in the atmosphere such that its pressure decreases from 
1000 to 800 mb while its temperature remains constant. 
Calculate the change in the specific entropy of the parcel. 
(64 Jkg -1 K -1 ) 

(5.10) During a process a parcel of dry air descends from 900 
to 950 mb and its specific entropy decreases by 30 Jkg -1 
K” 1 . If its initial temperature is 273 K what is (a) its 
final temperature and (b) its final potential temperature. 
(269 K, 273 K) 

(5.11) Show that for flows on a constant pressure level the relative 
change (i.e. percent change of original value) in temperature 
is equal to the relative change in potential temperature. 

(5.12) Cooler air moves over a warmer surface. What will happen 
to the surface pressure? (Hint: the condition of problem 5.11 
does not apply here.) 

(5.13) Show that the first term on the right-hand side of equation 
(5.32) is greater than zero except when T ) = T 2 when it is 
zero. (Hint: define b = C Vl /C V2 and x = T 2 /T 1 and find the 
first and second derivatives of /( x) = A S T /C Vl .) 

(5.14) Try to explain to your friends why the second law of ther- 
modynamics implies a positive arrow of time (i.e. time goes 
only forwards). 

(5.15) Consider two isolated systems existing very far from each 
other. The same irreversible transformation occurs in both 
of them. Will the change in entropy be the same in both of 
them? 

(5. 16) In an isentropic process the specific volume of a sample of dry 
air increases from 300cm 3 g _1 to 500cm 3 g _1 . If the initial 
temperature was 300 K, what are the final values of T and 
pi (244.5 K, 1404 mb) 

(5.17) Prove geometrically (in a p, V diagram) equation (5.2). 




CHAPTER SIX 


Water and its 
transformations 


The fundamental equations derived in the previous chapter can only 
be applied to closed systems (i.e. to systems that do not exchange 
mass) that are homogeneous (i.e. they involve just one phase). In 
such cases we do not need to specify how thermodynamic functions 
depend on the composition of the system. Only two independent 
variables (T and p , or p and V, or p and T ) need to be known. 
Since the total mass ( m ) remains constant, if we know the values 
of the extensive variables per unit mass we can extend the equa- 
tions to any mass by multiplying by in or by n (the number of 
moles). 

A heterogeneous system involves more than a single phase. In this 
case we are concerned with the conditions of internal equilibrium 
between the phases. Even if the heterogeneous system is assumed to 
be (as a whole) a closed system the phases constitute homogeneous 
but open “subsystems” which can exchange mass between them. In 
this case the fundamental equations must be modified to include 
extra terms to account for the mass exchanges. These extra terms 
involve a function p called the chemical potential, p = p(T,p). 
We will not go into the details of defining p\ we will only accept 
that in the case of open systems something else must be included 
to account for the heterogeneity of the system. In this book we 
are concerned with a heterogeneous system that involves dry air 
(N 2 , 0 2 , C0 2 , Ar) and water, with the water existing in vapor 
and possibly one of the condensed phases (water or ice). Then one 
component of the system is the dry air (which is assumed to remain 
always unchanged and in a gaseous state) and the other component 
is the water substance which can exist in two phases. Thus dry air 
is a closed system whereas the two phases are two open systems. 
The system composed of dry air, water vapor, liquid water, and ice 
is also of interest but it is quite unstable and as a result it cannot 
exist in equilibrium. 
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6.1 Thermodynamic properties of water 

In this chapter we will put aside dry air and we will concentrate 
on the one-component heterogeneous system “water” comprised of 
vapor and one of the condensed phases. Water vapor, like dry air, 
can be treated to a good approximation as an ideal gas. Accordingly, 
water vapor obeys the equation of state and if it existed only by 
itself its states would be determined by two independent variables. 
When, however, water vapor coexists with liquid water or ice things 
become a little complicated (or less complicated depending on the 
point of view). In this case, the mixture does not constitute an ideal 
gas and the equations applicable to ideal gases do not apply. One 
phase requires two independent variables ( p and T) and the other 
phase two more {p' and T'). Then, equilibrium between the phases 
will require that p = p' and T = T' . However, since the mass of 
each open subsystem does not remain constant another criterion for 
equilibrium must be considered. This criterion is given by p = p' 
and expresses equilibrium in the terms responsible for the mass 
exchange. Accordingly, for equilibrium between the two phases the 
following three constraints must be satisfied: 

P = P' 

T = T 

p = p'. 

The first two equations reduce the number of independent variables 
to two (say p and T). Then, the third equation p — p(p,T) = 
p'(p', T') = p'(p, T) reduces the number of independent variables by 
one more. Since p(p,T ) = p'(p,T ) it follows that p = f{T). Thus, 
if we fix the temperature at which the phases exist in equilibrium 
the value of pressure also becomes fixed. This defines curves along 
which equilibrium between two phases can exist. We conclude that 
in the case of a one-component system involving two phases in 
equilibrium the number of independent variables (otherwise known 
as degrees of freedom) is one, not four. 

Let us now consider the case of a one-component system involving 
all its phases (vapor, liquid, solid). Following the same arguments 
as above the following constraints must be satisfied for equilibrium 

p = p' = p" 


p = p' = p”. 

Then the first two equations reduce the number of independent 
variables from six to two. The first equality of the third equation 
reduces them by one more and the second equality by yet one more. 
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Thus, we remain with zero independent variables, meaning that all 
values are fixed, which implies that coexistence of all phases at equi- 
librium requires just a point, called the triple point, in the p, V, T 
state space. It follows that the greater the number of phases the 
fewer the degrees of freedom. The following formula suggested by 
Gibbs generalizes the above as it gives the number of independent 
state variables, N, for a heterogeneous system involving C different 
and non-reactive components in a total of P phases: 


N = C + 2 — P. (6.1) 

Figure 6.1 shows the triple point for water as well as the curves 
p = f(T) representing the phase transition equilibria. We call 
the three equilibria vaporization (gas -O- liquid), fusion (liquid O 
ice), and sublimation (gas -o- solid). Along the curves for vaporiza- 
tion and sublimation, vapor is in equilibrium with water and ice, 
respectively. As such these curves provide the equilibrium vapor 
pressure for water and ice. Vapor in equilibrium with a condensed 
phase is often called saturated and the corresponding equilibrium 
vapor pressure is called saturation vapor pressure. There is no real 
distinction between the two terms. Note that the extension of the 
vaporization curve to temperatures below the triple point corre- 
sponds to supercooled water, with vapor pressure over supercooled 
water being greater than that over ice. This is a metastable equi- 
librium where supercooled water and water vapor coexist. Here the 
system of supercooled water and water vapor may be stable with 
respect to small changes in temperature and pressure but intro- 
duction of ice in the system makes it unstable and as a result water 
freezes. The triple point corresponds to p t = 6.11 mb and T = 273 K 
(more accurately 273.16 K). Because of the different densities of 
liquid water, ice, and water vapor, the specific volumes for liquid 


Figure 6. 1 

Phase-transition equilibria 
for water. 
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water, a w , ice, a i? and water vapor, a v , at the triple point are: 

a„ = 1.000 x 10 -3 m 3 kg -1 
a ; = 1.091 x 10 -3 m 3 kg -1 
a v = 206 m 3 kg -1 . 


An interesting feature of Figure 6.1 is that the vaporization curve 
ends at point C where the temperature (T c ) is 374 °C and the pres- 
sure ('Pc) is approximately 2.21 x 10 5 mb. Beyond this critical point 
there is no line to separate liquid from vapor. Otherwise stated, 
beyond point C there is no discontinuity between the liquid and 
the gaseous phase. This indicates that beyond point C we cannot 
differentiate between the liquid and the vapor phase. In order to 
understand the significance of critical point C let us consider point 
A representing the temperature and the pressure of a sample of 
water vapor below critical point C . If we decrease the volume while 
keeping the temperature constant the pressure increases until it 
becomes equal to the pressure at point A 1 where it will condense. 
Similarly, if we keep the pressure constant and cool the gas until its 
temperature becomes equal to the temperature corresponding to 
point A 2 , the gas will condense. Obviously the gas will condense if 
we perform any intermediate process described by line AA 3 . Thus, 
there are three ways by which the gas can become liquid. We either 
have to cool it at constant pressure or compress it at constant 
temperature or simultaneously compress it and cool it. All three 
ways require crossing the equilibrium line. The fact that there is 
an end to the vaporization curve means that above T c we cannot 
liquefy a gas by compressing it at a constant temperature and that 
above pc we cannot liquefy a gas by cooling it at a constant pres- 
sure. Because in our atmosphere temperature and pressure are well 
below the critical point, water vapor can condense. For the other 
gases in the atmosphere (N 2 , 0 2 , Ar) the critical temperatures are 
very low (— 147 °C, — 119 °C, — 122 °C, respectively) and so they do 
not condense. 

Let us now consider phase changes along isotherms in the ( p , V) 
domain. Such a domain (shown in Figure 6.2) is called an Amagat- 
Anclrews diagram. We start with a sample of water vapor at a state 
corresponding to point A in Figure 6.1 (i.e. at a temperature X) and 
pressure p 1 greater than the triple point temperature T t and pres- 
sure p t . If we compress the vapor isothermally the pressure increases 
until we reach point A 1 . At this point liquid water and water vapor 
coexist in equilibrium. This means that some water vapor has con- 
densed to form liquid water and that liquid and vapor are in equilib- 
rium. Now what will happen if we further compress the vapor? Since 
we are on the equilibrium curve the pressure of the vapor depends 
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Figure 6.2 

An Amagat-Andrews 
diagram. 



only on temperature. As such, since the transformation is isother- 
mal, as long as liquid and vapor coexist a further decrease in volume 
cannot change the pressure. Do not get confused here. The ideal gas 
law does not apply to the mixture. Thus, if we keep on compress- 
ing the vapor we will observe that more vapor condenses to form 
liquid water (with liquid and vapor being always in equilibrium) 
until we reach point B where all vapor is condensed. This part will 
be represented in the (p, V) diagram by the horizontal (constant 
p) stretch A 1 B. A further (small) compression causes the liquid 
water’s pressure to increase rapidly (owing to the low compress- 
ibility of liquids) . This part is indicated by the segment BB 1 . The 
complete isotherm is depicted in Figure 6.2 by the line AA 1 BB 1 . 
For higher temperatures we observe that the horizontal stretch of 
the isotherms decreases until it reduces to a point. This is the point 
C corresponding to the critical state T Cl pc,V c - The isotherms 
above the critical temperature are monotonic decreasing functions 
with no discontinuity. They go over into equilateral hyperbolas 
which is how ideal gases behave. For temperatures below Ti the 
horizontal stretch of the isotherms increases. A solid line in Fig- 
ure 6.2 connects the beginning and end points of the horizontal 
stretches. This line together with the critical isotherm and the 
isotherm corresponding to the triple point temperature ( T t ) par- 
tition the diagram into six regions: vapor, gas, liquid water, ice, 
liquid water plus vapor, and ice plus vapor. Here the distinction 
between gas and vapor reflects only the fact that above T c no 
two-phase discontinuity for condensation exists. 


6.2 Equilibrium phase transformations - latent heat 

In the case of a homogeneous system undergoing an isobaric 
transformation that does not result in a phase change the heat 
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exchanged is proportional to the temperature change (the 
proportionality constant is C p ). In the case of a heterogeneous sys- 
tem involving two phases we say that when the two phases are in 
equilibrium a fixed temperature implies a fixed pressure, meaning 
that isobaric transformations are also isothermal transformations. 
Thus, even though the masses of the two phases change, the tem- 
perature of the system does not (recall the horizontal stretches in 
Figure 6.2). In such cases, the amount of heat exchanged depends 
only on mass changes (which modify the internal energy of each sub- 
system) and on the work done owing to possible volume changes 
(the volume can change because the mixture does not obey the ideal 
gas law). By definition, the latent heat, L, of a transformation is 
the heat absorbed (or given away) by the system during an isobaric 
phase transition 

L = 6Q P = constant. 

Recall that enthalpy is defined as H = U + pV. Differentiation of 
this equation gives dH = dU +pdV+Vdp. Since the change of phase 
takes place at a constant pressure it follows that dH = dU + pdV. 
The right-hand side of this equation is according to the first law 
equal to 5Q. It follows that 


L = dH. (6.2) 

Thus, the latent heat of a phase change is the change in enthalpy 
during the transformation. The latent heats of vaporization (liquid 
FT vapor), fusion (ice -FT liquid), and sublimation (ice -FT vapor) 
are denoted by L V ,L { , and L a , respectively. Note that the above 
latent heats are positive when during the transformation heat is 
absorbed (liquid — > vapor, ice — > liquid, ice — > vapor) and negative 
when heat is released (vapor — > liquid, liquid — > ice, vapor — > ice). 
L v ,L f , and L s are also called the enthalpy of vaporization, fusion, 
and sublimation, respectively. According to the above, 

L v = H v — 7T W = U v — U„ + p wv (K — V w ) 

U = - H, = U w - U, + p wi (y w - VC) (6.3) 

L s = H v -Hi = U v - U,+p vi (V v - Vi) 

where here p stands for the equilibrium pressure between two 
phases. At the triple point p wv = p vi = p wi = 6.11 mb. It follows 
that at the triple point the specific latent heats l v ,l c , and l s satisfy 
the relationship l a = l f + l v . Obviously, the latent heats depend on 
temperature (see Table A. 3). Nevertheless, in the range of temper- 
atures observed in the troposphere they do not vary significantly. 
Because of that it is often assumed that the latent heats are inde- 
pendent of temperature. It is easy to show that the first law in 
the case of an isobaric change of phase can be written as dU = 
L — pdV. For fusion dV « 0 and for vaporization and sublimation 
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dV = Kapor-Miquid water or ice « Kapor ■ Thus, the first law reduces to 


dU = L — mR v T for vaporization or sublimation (6.4) 
where R v is the gas constant for water vapor. 

6.3 The Clausius-Clapeyron (C-C) equation 

Before we proceed and to avoid confusion with too many subscripts 
we need to make a change in our notation. From now on we shall 
represent the water vapor pressure as e. Furthermore, during phase 
transitions the equilibrium (saturation) vapor pressure over water 
will be denoted as e sw and equilibrium vapor pressure over ice will 
be denoted as e si . 

For an ideal gas we saw how the equation of state p = f(V,T ) 
relates changes of pressure, temperature, and volume. In the case 
of a heterogeneous system involving two phases we saw that p = 
e s = f(T). Does this indicate that a relationship similar to the 
equation of state for ideal gases exists for heterogeneous systems? 
Let us consider the transformation described by isothermal 7\ in 
Figure 6.2 and focus our attention on segment A 1 B. Going from 
Ax to B a change of phase takes place (vapor — > liquid). As we 
discussed before, between A x and B vapor is in equilibrium with 
liquid water and p = e sw = f{T). Because of that any other state 
variable or function is a function of temperature only. Let us assume 
that at point C (Figure 6.3) the mass of the liquid phase is to w and 
the mass of the vapor phase is m v . Then the total volume and 
internal energy of the system at point C are 


dU = L 


for fusion 


V = m w a w + m v a v 
U = m w u w + m v u v 


Figure 6.3 

The isotherm T± of Figure 
6.2. For any given change 
C — > C' the change in 
volume is dV. This causes 
an amount dm of liquid 
water to evaporate. 


P (mb) 



V 
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where a w ,a v ,u w ,u v are the specific volumes and internal energies 
of the two phases. Now assume that the system moves from C to 
C' . This change corresponds to a change in volume dV and causes 
an amount dm of liquid to evaporate. The volume at point C' will 
be 

V + dV = (rn w — dm)a„ + (?n v + dm)a v 
or using the previous equation 


dV = ( a v — a w )dm. 


Similarly 


dU = (u v — u„)dm. 


Recall from equation (6.3) that during this change 


(6.5) 

( 6 . 6 ) 


U v U V.' C sw (tt v Ct w ) I y 


(6.7) 


Dividing equation (6.6) by (6.5) yields 

dU _ u v — u w 
dV a v — a w 

or using equation (6.7) 


dU _ K 
dV a v — a w 


( 6 . 8 ) 


Here is as good a place as any to restate that while general defini- 
tions such as the first law (equation (4.5)), the second law (equation 
(5.12)), enthalpy ( H = U+pV ), specific heat (equation (4.14)), etc. 
are valid for any system, variations or expressions derived using the 
ideal gas law (for example, equations (4.22) and (4.18)) are not valid 
for any system but ideal gases. According to problem 6.1, for any 
system 



(6.9) 


Also, according to problem 6.2, for an ideal gas (dU/dV) T = 0 
(Joule’s law). However, water in equilibrium with its vapor is not 
an ideal gas and U and e sw are functions of temperature only. Then 
since the transformation is isothermal we can write that 


(dU\ _ dU 
\WJ T - dV 


and 



de 


sw 


dT ' 


Accordingly, equation (6.9) becomes 


dU rpdes-w 
dV = ~dT 
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Combining this equation with equation (6.8) yields 


T 


dc 

u, '-'SW 

~dT 


6 S - 


- e a 


or 


dT T(a v — a w ) 


( 6 . 10 ) 


Equation (6.10) is called the Clausius-Clapeyron (C-C) equation. 
In its general form, 


de s l 
~dT = TAa 


( 6 . 11 ) 


it relates the equilibrium pressure between two phases to the tem- 
perature of the heterogeneous system. Here l is the specific latent 
heat (or specific enthalpy) of the change of phase and A a is the 
difference in the specific volumes between the two phases at tem- 
perature T. It is for heterogeneous systems what the equation of 
state is for ideal gases. 


6.4 Approximations and consequences of the C-C 
equation 

• Temperature dependence of enthalpy of vaporization 

By differentiating equation (6.3) we have 

dL v _ dH v dH„ _ 

dT dT dT pv pw 

where C pv and C pw are the heat capacities at constant pressure 
for vapor and liquid water, respectively. Since L v is the enthalpy 
difference between the vapor and liquid water phases at equilibrium, 
it depends only on temperature. As such the above equation can 
be written as 

^ = C pv - C pw . (6.12) 

For a wide range of temperatures (— 20 °C to 30 °C), C pv and C pw 
vary very little (about 1%). Because of that we can consider them 
independent of temperature and integrate equation (6.12) to get 

L v = L v0 + ((7 pv — C pw )(T — T 0 ) 


or 


l v — lv o + ( c pv — c pw )(T — T 0 ) (6.13) 

where l v0 is the specific latent heat of vaporization at the ref- 
erence state T 0 . For the reference state T 0 = 273 K with l v0 = 
2.5 x 10 6 Jkg- 1 ,^ = 1850 J kg " 1 K- 1 , and c pw = 4218 J kg " 1 K" 1 , 
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equation (6.13) provides a very good approximation of l v over the 
range from — 20 °C to 30 °C. 

We can repeat the above for sublimation and fusion to arrive at 


or 


dT a _ n _ n 
dT { ^ / p v 

dLf _ „ 
dT 


h — ho + (c P v — c pi )(T — Tg) (6.14) 

h = ho + ( c pw — c pi )(T ~ T 0 ). 


Recall that the specific heat for constant volume is c v = Sq/dT. 
Since water vapor is considered an ideal gas (recall Joule’s law), it 
follows that 

Sq du v 
Cvv= dT = ~dT 

For liquid water (not an ideal gas) 


Cvw — 


Sq_ 

dT 


du w 

~dT 


da w 

} ~dT' 


Since a w varies very little with T the above equation reduces to 


Cy w 


du w 

~dT' 


Similarly, for ice c vi = duJdT. Now since h = u + pa it follows that 
for water vapor c pv = dh v /dT = du v /dT + i? v . Therefore for water 
vapor c pv ^ c Vv . Liquid water and ice do not behave as ideal gases. 
Therefore pa ^ RT. In this case 

_ dh„ _ du„ da w dp ~ du w 

Cpw ~ It ~ hr + p lr + a ™df ~ It 

(because dp = 0 and da w /dT « 0). Thus, c pw ~ c Vw = c„. Similarly 
c pi « Cyi = Ci. Using the values of c pv ,c Vv ,c w , and c ; at T = 
0 °C (1850, 1390, 4218, and 2106 J kg -1 K -1 , respectively) and using 
equations (6.13) and (6.14) we find that 


dl v 

dT 


-2368 J kg” 1 K" 1 


dl s 

dT 


—256 J kg -1 K _1 


^«2112Jkg- 1 K- 1 . 

These rates are small compared with the values of l v ,l sl and l { at 
0°C (2.5xl0 6 , 2.834xl0 6 , and 0.334xl0 6 Jkg- 1 K _1 , respectively). 
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We can thus conclude that all latent heats vary very little with 
temperature. As such the above equations justify approximating 
l v , l a , and If as constants. An improvement to this approximation is 
given by equations (6.13) and (6.14) which express them as linear 
functions of T. 


• Temperature dependence of equilibrium (saturation) vapor 
pressure 


In the case of vaporization, a v 
approximated as 

dc 

u,c 'SW 

~dT 


a w and C-C equation can be 


Ta v 


(6.15) 


By combining equations (6.15) and (6.13) and considering vapor as 
an ideal gas (i.e. e sw a v = R V T ) we obtain 


1 de sw l v o ((-'p w Cp V )To 


dT 


R V T 2 


R V T 


or 


de s . 


’ e S w at To '-aw 


1-vQ (Cp-w c pv ) T 0 ( 

f T dT c - c 

f \^p w '-'pV 

f T dT 

R v J 

To T 2 R v J 

~p~ 
To 1 


or 


I e sw l v o "T ( ('jtv: c pv )T 0 f 1 1 \ 

e s o R v \T 0 T ) 


c — c T 

W '-'pV 


iL 


T 0 


Considering as the reference state the triple point where T 0 = 
0°C, e s0 = 6.11 mb, l v0 = 2.5 x 10 6 Jkg-\ c pv = 1850 4 kg- 1 
K _1 , and c pw = 4218 Jkg _1 K _1 the above equation reduces to 


In ^ = 6808 
CsO 



1 


T 


T 


- - 5.09 In — 


T 0 


or 

( 6808 \ 

53.49 — 5.09 InT J . (6.16) 

The above equation provides the relation between e sw and T. 
Another way to establish such a relation would be to integrate equa- 
tion (6.15) assuming that l v is constant and independent of T. This 
would result in an approximation to the above equation given by 


1 ^SW 

^sO 


R v Tq R v T 


or 


6.11 exp 


^19.83 



(6.17) 


In both equations (6.16) and (6.17) e sw is in mb and T is in Kelvin. 
In problem 6.3 you are asked to show that equations (6.16) and 
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(6.17) are nearly identical for temperatures in the range — 20 °C to 
30 °C. 

The corresponding relation to equations (6.16) and (6.17) for subli- 
mation can be derived following the same procedure (here again 
o v a;) with l s = 2.834 x 10 6 Jkg _1 and Ci = 2106 J kg -1 
K -1 . They are: 


6293 

e si = 6.11 exp ( 26.16 — 0.555 In T 


and 


6.11 exp 


^22.49 



(6.18) 


(6.19) 


The C-C equation and its simplified versions describe the thermo- 
dynamic state of water vapor and liquid water (or water vapor and 
ice) when the two phases are in equilibrium. When the system’s 
pressure is greater than the equilibrium pressure, equilibrium is 
approached by condensation (sublimation) of water vapor. This 
reduces the amount of vapor in the system thereby reducing the 
vapor pressure. If the system’s pressure is lower than the equilib- 
rium pressure, equilibrium is approached by water (ice) evaporating 
(sublimating) . Because of the particular exponential dependence of 
e sw and e si on T in equations (6.17) and (6.19) it follows that for 
higher temperatures e sw is greater than for lower temperatures. As 
such more water vapor can exist at higher temperatures. For exam- 
ple, for T = 300 K from equation (6.17) we have e sw = 36 mb. For 
T = 273°C, e sw = 6 mb. This is a substantial difference. Note that 
this does necessarily imply that warmer air can hold more water 
vapor than colder air as it is often wrongly stated. All the above 
arguments and formulas were derived in the absence of air! We only 
considered two phases with water evaporating (or vapor condens- 
ing) into a vacuum. If you think about it this is consistent with 
Dalton’s law which would indicate that in a mixture of air and 
water vapor the partial pressures are independent of the mixture. 
Therefore, keep in mind that when we say that “air at a given tem- 
perature is saturated” we do not imply that the air holds as much 
vapor as it can hold. We strictly mean that at that temperature 
the amount of vapor that can exist is maximum regardless of the 
presence of air. 

The reason for the relation between equilibrium vapor pressure 
and temperature arises from purely kinetic reasons. At equilibrium 
evaporation is zero (strictly speaking it is the net evaporation = 
evaporation minus condensation that is zero as molecules inside 
water have as much chance to evaporate as vapor molecules have 
to condense). At higher temperatures the molecules inside the liq- 
uid acquire greater speeds and their chance to escape increases. As 
a result (net) evaporation increases. 
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The exponential dependence of e sw (or e si ) on T has important con- 
sequences for our climate system. Tropical warm water can transfer 
more water vapor into the atmosphere than colder extratropical 
water. Most of this water vapor is then precipitated as a result of 
organized convection in the tropics thereby releasing large amounts 
of heat. This heat can then be converted to work and generate 
kinetic energy which helps maintain the general circulation against 
frictional dissipation. 


• Changes in the melting and boiling points for water 


At boiling point (T 
1.673 m 3 kg -1 , and a w 
it follows that 


= 100°C),/ V = 2.26 x 10 6 Jkg-\ a v = 
= 0. 00104m 3 kg -1 . Then from equation (6.15) 

^ = 3621 PaK” 1 


or 


df> 

u, '-'SW 

~dT 


0.03575 atm K" 1 . 


Accordingly, a decrease in pressure by 0.03575 atm (36.2 mb) lowers 
the boiling point by 1°C. Note that if we evaluate the approximate 
relationship (6.16) at Tj = 373 K and T 2 = 372 K, we will find 
that e sw (Tj) « 1003 mb and e sw (T 2 ) « 969 mb, or Ae sw « —34 mb. 
If we use instead equation (6.17), we will find that e sw (Ti) ss 
1237 mb, e aw (T 2 ) ss 1190 mb, and Ae sw ss —47 mb. The differences 
are due to the fact that in deriving the approximate expressions 
we have assumed that the specific heats and latent heats are either 
constant or linear functions of T. While in the range of temper- 
ature relevant to weather these approximations may be good, for 
temperatures outside this range the differences may be significant. 
At melting point (T = 0°C), l { = 0.334 x 10 6 Jkg _1 , a w = 
1.00013 x 10 -3 m 3 kg -1 , and a; = 1.0907 x lO^nPkg -1 . Substi- 
tuting these values into equation (6.15) yields 


dpwi 

dT 


— 13 503 800 PaK -1 


or 

= — 133.3atmK _1 . 
dl 

This result indicates that lowering the melting point by just 1°C will 
require a tremendous increase in pressure. This is actually rather 
anomalous as in most cases the melting point increases with increas- 
ing pressure. The strange behavior of water and ice is due to the 
fact that ice is less dense than water whereas in most cases the solid 
is denser than the liquid. The fact that the melting point of ice is 
lowered by increasing pressure is very important to geophysics as 
it explains the motion of glaciers. When a mass of ice encounters 
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a rock on the glacier bed the high pressure of the ice against the 
rock lowers the melting point of ice causing the ice to melt and 
then to flow around the rock. Once behind the rock the pressure 
is restored and ice freezes again. In this way ice propagates around 
obstacles. 


Examples 

(6.1) Calculate the work done, the amount of heat absorbed, and the 
change in the internal energy during (a) melting of I g of ice to 
water at a temperature of 0°C and under a constant pressure 
of I atmosphere, and (b) evaporation of I g of water to water 
vapor at a temperature of I00°C and under a constant pressure 
of I atmosphere. 

(a) When the effect of absorbing an amount Q at a con- 
stant pressure is a change in the physical state of the system 
(which also occurs at constant temperature), then Q is pro- 
portional to the mass that undergoes the change. In our case 
Q = mlf where Zf is the specific latent heat of fusion of 0 °C. 
For ice Z f = 79.7calg -1 . Thus, 

Q = Z f TO = 79.7 cal. 

During melting the volume changes. The corresponding work 
done is 

/♦water /“water 

W= pdV = P dV= Plater - ^ ce ) • 

J ice J ice 

The density of ice at 0°C is 0.917 gem -3 , so lg occupies a 
volume of 1.09 cm 3 . The density of water at 0 °C is 1 g cm -3 . 
Accordingly, lg of water occupies a volume of 1.0cm 3 . It 
follows that 

W = 1013 x 10 2 x (1 - 1.09) x 10 -6 J 


or 

W = —0.00912 J 


or 

W = -0.0022 cal. 

The work done is negative because the volume of the result- 
ing water is smaller than the volume of ice. From the above 
we then find that 


AZ7 = Q — W = 79.7022 cal. 
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(b) In this case at 100 °C, l v = 539 cal g -1 . Thus, 

Q = l v m = 539 cal. 

Water at 100 °C has a density of 0.958 gem -3 which for lg 
gives a volume of 1.04 cm 3 . The volume of 1 g of water vapor 
at 100 °C is 1673cm 3 . As such 

W = P (V v apor - V water ) = 1013 x 10 2 x (1673 - 1.04) x 10 -6 
= 169.37 J = 40.47 cal. 


It then follows that 


AU = 498.53 cal. 

Note that in both cases the internal energy increases even 
though the temperature remains constant. This is due to 
the fact that the degrees of freedom increase. The degrees 
of freedom of the molecules in H 2 0 is much greater in the 
liquid phase than in the ice phase and much greater in water 
vapor than liquid water. 


(6.2) Can raindrops grow from cloud droplets through condensation 
alone? 

In order to address this question we have to make some rea- 
sonable assumptions. Cloud droplets form over nuclei and 
begin to grow according to the available water vapor. If the 
vapor pressure is greater than the equilibrium vapor pres- 
sure droplets grow by condensation. So we need to get an 
idea of how much vapor is available and how many droplets 
are going to compete for this vapor. Let us assume that 
the concentration of nuclei is 100 cm -3 which corresponds 
to rather clean air. Now we need to estimate the amount 
of vapor that will be available during the process of grow- 
ing by condensation. We can get a gross estimate of this 
amount by assuming that at the level where clouds form (say 
1000 m) the air is saturated and the surface temperature is 
about 30 °C. If this air ascends adiabatically, the tempera- 
ture decreases. As the temperature decreases the equilibrium 
(saturation) vapor pressure decreases. This means that not 
as much vapor is needed for equilibrium or that the vapor 
pressure becomes greater than the equilibrium vapor pres- 
sure. Thus, the droplets begin to grow. If the air ascends to 
the maximum level (say the top of the troposphere at about 
16 km), then we can estimate the equilibrium vapor pressure 
at that level by assuming that no water vapor condenses 
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using the equation 


In 


^SW 

tm 


19.83 


5417 

~f~’ 


where T= 30 - (16 x 9.8) = -126.8 °C = 146.2 K. It follows 
that at this level e aw = 0. At cloud base where the temper- 
ature is ss 20.2 °C, e sw ss 23.7 mb. As such in this scenario 
the equivalent of 23.7 mb in water vapor will be available 
to the droplets for growth. How much water vapor is that? 
From the ideal gas law we find that for V = 1 m 3 

m v = = 0.0175 kg. 

Assuming that a typical raindrop has a radius of 0.5 mm it 
follows that the mass of a raindrop is m rain dro P = V/p ~ 
0.00052 g. Given the total amount of water vapor available 
this translates to ~ 33 650 raindrops. However, in 1 m 3 of air 
there are 10 s “competitors” with equal rights. Simply, none 
of them will grow to be a raindrop by condensation only. 
Other processes must also take place (collision, coalescence, 
for example). 


(6.3) Consider a sample of 2 moles of supercooled water (i.e. liquid 
water existing at environmental temperatures less than 0 °C). 
Assume that the temperature of the surroundings is — 10 °C 
and that the supercooled water freezes. In this process the 
latent heat of fusion (freezing) is lost to the surroundings and 
finally ice and surroundings return to the original temperature. 
Calculate AU, AH, and AS for both supercooled water and 
surroundings. 

Two moles of water equal 36 g. If we neglect volume and 
pressure variations during phase changes between water and 
ice it follows from the first law and the definition of enthalpy 
that for both supercooled water and surroundings 

AH = AU = Q. 

There is a correct and a wrong way to calculate Q and AS 
during this process. If we calculate the amount of heat given 
away from the supercooled water as it freezes at —10 °C, we 
find that Q w = -mZ f (-10°C) = -0.036 x 0.312 x 10 6 J = 
—11 232 J. Then, if we use equation (5.7) we will find that 
A S w = J SQ/T = — 42.7JK -1 . The amount of heat released 
is gained by the surroundings. Thus, Q s = 11 232 and AS s = 
42.7 J K _1 . It follows that the change in entropy of the whole 
system (surroundings + supercooled water) is zero. But there 
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was a change in the system, which, according to the sec- 
ond law, should increase the entropy of the system. The 
reason for this violation is that we have indirectly assumed 
that the spontaneous freezing is a reversible process, which 
is not. 

A correct way would be the following: Start with the super- 
cooled water at —10 °C and warm it reversibly to 0 °C (step 
1), freeze it reversibly at 0°C (step 2) and then cool it 
reversibly to — 10 °C (step 3). In this scenario we have 


Step I 

Q 1 = mc pw AT= 1518.5 J 
r T 2 d t 

ASi = / mc pw — = 5.67 JK -1 
Jt-l 1 

Step 2 

Q 2 = —mlf(0 °C) = -12 013.2 J 

A S 2 = — = -44.0 JK" 1 
T 2 

Step 3 

Q 3 = mc pi AT = — 758.2 J 
f Tl dT 

A S 3 = / mc pi — = -2.83 JK -1 
Jt 2 1 

In the above calculations we have considered that Ti = 
263K,T 2 = 273 K, c pw = c w = 4218 J kg -1 K _1 , c pi = 
Ci = 2106 J kg -1 K- 1 , and Z f (0°C) = 0.3337 x lO^kg" 1 
K _1 . It follows that the amount of heat released by the 
supercooled water is 

Qw = Qi + Q 2 + Qs = -11 252.9 J, 

and the change in entropy of the supercooled water is 

A S w = AS 3 + AS 2 + AS 3 = -41.16 JK" 1 . 

The surroundings (T s = — 10 °C) gain the amount that the 
freezing of supercooled water releases. Thus, 

A U s = A H s = Q S = — Q w = 11 252.9 J. 

At the same time their change of entropy is 

AS s = / ^ = ^= 42.79 JK' 1 . 

J s -^s 
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It follows that the change in the entropy of the whole system 
(surroundings + supercooled water) is 

AS = AS s + A S„ = 1.63 J KT 1 > 0. 


Problems 

(6.1) Prove that for any system 



(Hint: express dU using equation (2.1) and substitute into 
the first law. Then consider the definition of dS = SQ/T and 
the fact that dS is an exact differential, perform the partial 
differentiations and collect terms.) 

(6.2) Using the equation derived in problem 6.1, show that for an 
ideal gas U is only a function of temperature and does not 
depend on V (Joule’s law). 

(6.3) By plotting e sw against T show that in the range — 20 °C to 
30 °C equations (6.16) and (6.17) are virtually identical. If 
we were very picky about choosing one over the other which 
would you choose? 

(6.4) Inside an open container is 500 g of ice with temperature 
T = — 20 °C. Heat is then provided to the ice at a rate of 
100 cal min -1 for 700 min. Plot (a) the relationship between 
temperature (in °C) and time (in seconds), (b) the relation- 
ship between the heat absorbed by the ice and time, (c) 
the relationship between the heat absorbed and tempera- 
ture. Neglect the heat capacity of the container and assume 
that If = constant = 79.7 cal g -1 and c ; = constant = 
0.503 calg -1 K -1 . 

(6.5) A closed box of volume 2 m 3 and temperature 120 °C con- 
tains 1 m 3 of water saturated with vapor. Calculate (a) the 
mass of water vapor, (b) the mass of water vapor that must 
slowly escape in order for the temperature inside the box to 
decrease to 100 °C while maintaining saturated conditions. 
(1.43 kg, 0.71 kg) 

(6.6) An open container of volume 2 liters contains dry air and 
some quantity of water at standard conditions. After equi- 
librium is achieved the container is closed and heated to 
a temperature of 100 °C. What is the least amount of 
water that must have been added if after the heating the 
air is saturated? What is the final pressure inside the box? 
Use equation (6.16). (1.16 g, 2.35 atmospheres) 

(6.7) A substance obeys the equation of state pV 12 = 10 9 T 11 . 
A measurement of its internal capacity inside a container 
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Figure 6.4 

Illustration for problem 
6 . 10 . 


Ice 



having a constant volume of 100 liters shows that under 
these conditions the thermal capacity is constant and equal 
to O.lcalK -1 . Express the energy and the entropy of the 
system as function of T and V. (Hint: Is the substance an 
ideal gas?) 

(6.8) How would you have to alter your assumptions in example 
6.2 in order to make cloud droplets grow to raindrops by 
condensation only? 

(6.9) Calculate the change in entropy when 5 kg of water at 10 °C 
is raised to 100 °C and then converted into steam at that 
temperature. At 100 °C the latent heat of vaporization of 
water is 2.253 x 10 6 Jkg -1 . In the above range of tempera- 
tures, the specific heat capacity under constant volume for 
water is nearly independent of temperature with a value of 
about 4.18 x 10 3 Jkg^K- 1 . (> 35 972JKT 1 ) 

(6.10) A block of ice rests as shown in Figure 6.4. A wire on which 
we attach a weight surrounds the ice. Explain why the wire 
will go through the ice block without the ice breaking into 
two pieces. 

(6.11) The warm source of a Carnot cycle has a temperature of 
100 °C while the cold source consists of melting ice. When 
the cycle operates for 1 h it is observed that 1 ton of ice 
has melted. Calculate (1) the heat absorbed in the cycle 
from the warm source (2) the heat given away in the cycle 
to the cold source. The latent heat of fusion at 100 °C is 
0.3337 x 10® Jkg" 1 . (333.7 x 10® J, 244.2 x 10 6 J) 
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Moist air 


Our atmosphere is basically a two-component system. One compo- 
nent is dry air and the other is water existing in vapor and possibly 
one of the condensed phases (liquid water or ice) . According to Dal- 
ton’s law, in a mixture of ideal gases each gas can be assumed to 
behave as if the other gases were absent. As such, in a mixture of dry 
air, water vapor, and a condensed phase, the “water” system (water 
vapor + condensed phase) can be treated as being independent 
of the dry air. In this case, the concepts developed in the previ- 
ous chapter (for the one-component heterogeneous system “water” ) 
are valid for the two-component heterogeneous system “dry air + 
water”. We will call the system consisting of dry air and water 
vapor “moist air” and it can be unsaturated or saturated with water 
vapor. Since liquid water is absent, moist air is a two-component 
system with one phase present. As such, according to equation (6.1) 
we need three state variables to specify the system’s state. Usually 
these variables are taken to be pressure, temperature, and a new 
variable (to be defined soon) called mixing ratio. If the condensed 
phase is present and in equilibrium with the vapor phase, then two 
variables are needed (typically temperature and pressure). 

For clarity in our notation we will use the subscripts d, w, v to indi- 
cate dry air, liquid water, and water vapor, respectively. The only 
exception will be for the vapor pressure, which we will denote sim- 
ply as e. Variables with no subscripts will correspond to a mixture 
of dry air and water vapor. 

The following table gives various systems, their type, and the 
variables needed to specify their state (see equation 6.1). 
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System 

Type 

Variables needed to 
specify system’s state 

Dry air and 
water vapor 

Heterogeneous - 
two components, 
one phase 

3 ( p , T mixing ratio) 

Dry air, water vapor, 
and liquid water 

Heterogeneous - 
two components, 
two phases 

2 (P,T) 

Dry air, water vapor, 
liquid water, and ice 

Heterogeneous - 
two components, 
three phases 

1 ( T ) 

Liquid water and 
water vapor 

Heterogeneous - 
one component, 
two phases 

1 ( T ) 

Liquid water, water 
vapor, and ice 

Heterogeneous - 
one component, 
three phases 

0 

Water vapor 

Homogeneous - 
one component, 
one phase 

2 (P,T) 

Dry air 

Homogeneous - 
one component, 
one phase 

2 (p,T) 


7.1 Measures and description of moist air 

7.1.1 Humidity variables 

In a sample of moist air, dry air and water vapor have the same 
temperature T and occupy the same volume V. Thus, for water 
vapor 

eV = m^RjT 
or 

e = p v R v T 

where e is the vapor pressure, p v is the density of water vapor, and 
R v is the specific gas constant for water vapor. Since the molecu- 
lar weight of water M v = 18.01 gmol -1 , the value of l? v is equal 
to R*/M v = 461.5 J kg -1 K -1 . Since R* = R d M d it follows that 
R d M d = R V M V . Thus, 
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M v 


e=TT = TT" = 0.622. 

R v M,i 

We define the specific humidity, q, and mixing ratio, w, as 


(7.1) 


q = 



m v 

m 


and 



m v 

m d 


(7.2) 


where m = m d +m v is the total mass of the mixture (p v is also called 
the absolute humidity). Since p v = eM v /R*T, p d = p d M d /R*T, 
and p d = p — e it follows that 


e 

w = e . 

p — e 

At saturation the saturation mixing ratio is 

e s 

w s = e 

P~e s 


(7.3) 


(7.4) 


where e s is either the equilibrium (saturation) vapor pressure with 
respect to liquid water (e sw ) or with respect to ice (e si ). In general 
e s , e <C p which reduces the above two equations to 

( r\ 

w~e-, w s ss e — . (7-5) 

P P 

Since 1/q = (m d + m v )/m v = (r?z d /m„) + 1 it follows that 


or 


w = 




w 


q = 


w 

1 + w 


(7.6) 


In the atmosphere both w and q are very small ( w , q <C 1). For this 
reason we can always assume that w « q. 

The relative humidity, r, is defined as 


. _ mv 

m vs 


where m v is the mass of water vapor in a sample of moist air of 
volume V and m vs is the mass of water vapor the sample would 
have had if it were saturated. Because of the ideal gas law we can 
write r as 


(7.7) 


Using (7.5) it follows that approximately 


w 

w s 
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7. 1 .2 Mean molecular weight of moist air and other quantities 

According to equation (3.13) we have that the mean molecular 
weight of moist air is 

— _ m d + m v 

M md mv ' 

M d ' M v 

We can manipulate the above equation by writing it as 


to d m v W 1 

M d M v ) lm d + m. 


M M d \m d + m 


1 1 m d 

AT M d m d + to v 


TO v /m d 

M v /M d 


1 _ 1 1 
W ~ M^l + w 


(-7) 


M M d V 1 + w e 


1 + w w 


1 — q ’ 


1 1 / q 


F = i 1+ (7- 1 > 


= = ITT (1 + 0.61q) . 

M M d y 

Then for moist air the equation of state is 

pa = R uunst T = =T = ^(1 + 0.61q)T, 
= f? d (l + 0.61q)T. 

The above equation defines the virtual temperature, 


T v i r t — (1 + 0.61q)T, 
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which can be interpreted as the temperature of dry air having the 
same values of p and a as the moist air. Otherwise stated, virtual 
temperature is the temperature that air of a given pressure and 
volume (or density) would have if the air were completely free of 
water vapor. Since q in reality is always greater than zero, T virt 
is always greater than T. The above equation also gives the gas 
constant of the mixture: 

R= (l + 0.61q)i? d . (7.10) 

Similarly, we can define the specific heat capacities of moist air. 
Suppose that a sample of such air receives, at a constant pressure, 
an amount of heat SQ which increases its temperature by dT. Some 
of this amount (i5Q d ) is absorbed by the dry air and some (SQ V ) by 
the water vapor. In this case we can write that 

5Q = SQ d + 5Q V 


or 


or 


or 


SQ = m d 6q d + m v 5q v 


X TO d r , m v r 

dq = Sq d H Sq v 

m d +m v m d + m v 


Sq = (1 - q)Sq d + qSq v . 

Recalling the definition of specific heat capacities we then have 

,, Jdd , Sq v 


or 


or 


or 


c p — (1 — q)c pd + q c p 


Cp £pd 


1 _|_ [ ypy — 1 ) q 

Cpd 


c p = c pd (l + 0.87q) » Cp d (l + 0.87u>). 


Similarly 


Cy — (1 — q)cy d + qcy v 

which leads to 

Cy = Cya(l T 0.97q) ~ Cy d ( 1 T 0.97w). 

From equations (7.11) and (7.12) we also have 

_ Cp _ c pd (1 + 0.87q) 

Cy Cy d (1 + 0.97q) 


(7.11) 


(7.12) 
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which because q <C 1 can be approximated with 
7 ss 7 d(l + 0.87q)(l — 0.97q) 
or after ignoring the second-order term 

7»7d(l-0.1q), (7.13) 

where 7 d = 1.4. Similarly, 

k _ 7-1 _ R _ R* (1 + 0.61q) 

7 c p c pd (1 + 0.87q) 

« fc d (l + 0.61q)(l — 0.87q) 
wfcd(l-0.26q), (7.14) 

where fc d = 0.286. Note that similar expressions in the presence of 
liquid water cannot be derived because liquid water + vapor does 
not constitute an ideal gas. 


7.2 Processes in the atmosphere 

In this section we will discuss scenarios that are relevant to basic 
processes in the atmosphere. To make the discussion flow effec- 
tively we will start with processes happening at a constant pressure 
level (for example, at the surface). After that we will proceed with 
processes involving ascent, first of unsaturated air and then of sat- 
urated air. Finally, we will discuss processes that involve mixing of 
two air masses in the horizontal and in the vertical. The laws and 
formulas describing each process will be derived analytically. Subse- 
quently, each process will be described visually in a (p, T) diagram. 
As we will see later, better diagrams can be designed. However, for 
now the simple and familiar (p, T) diagram will be used for illus- 
trating the various processes. In the following presentation when a 
condensed phase is present, we will assume that it is liquid water. 
The same theory holds when the condensed phase is ice as long as 
the appropriate substitutions are made (for example e si for e sw , l s 
for l v , Ci for c w , etc.). 


7.2.1 Isobaric cooling - dew and frost temperatures 

Consider a parcel of unsaturated moist air and assume it is a 
closed system. In this case q and w remain constant. If this parcel 
begins to cool at a constant pressure the vapor pressure e will also 
remain constant as long as condensation does not occur (recall that 
e = wp/(w + e). However, the equilibrium vapor pressure will not 
remain constant. As we saw in the last chapter both e sw and e si are 
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strongly dependent on T, and as T decreases so does e s . Since e s 
decreases it follows that the relative humidity increases and satura- 
tion is approached. The temperature at which saturation is reached 
is called the dew point temperature, T dew , if saturation is reached 
with respect to liquid water or the frost point temperature, T { , if 
saturation is reached with respect to ice. 

These temperatures derive their name from the fact that if the par- 
cel cools isobarically below these temperatures droplets (or ice) will 
form. However, keep in mind that for condensation (or sublima- 
tion) to happen at temperatures just below saturation point either 
condensation nuclei or solid surfaces must be present. Spontaneous 
condensation in the absence of nuclei or surfaces needs extremely 
high relative humidities - supersaturation - which cannot easily 
occur in the atmosphere. 

Isobaric cooling can happen naturally because of radiation cooling 
of the surface or the air itself. In this case the ground can lose 
enough heat to cool below the temperature at which atmospheric 
water vapor condenses. If this happens at temperatures greater than 
0°C water vapor in the layer very close to the surface condenses 
on the surface and becomes visible as dew. If it happens at tem- 
peratures less than 0 °C the product is frost. If the cooling is very 
strong or if the air itself cools by radiation, the layers of air above 
the surface may also cool below their saturation point leading to 
the formation of radiation fog. Alternatively, isobaric cooling can be 
achieved when air moves in the horizontal ( p w constant) and trav- 
els over colder regions. Then saturation may occur as the air cools 
by heat conduction to the surface. This leads to advection fogs. 

At T dew vapor pressure becomes the equilibrium (saturation) vapor 
pressure (and to a very good approximation the mixing ratio 
becomes the saturation mixing ratio). As such if e is the partial 
vapor pressure at T, the dew point must satisfy the equation 

e sw (T dew ) = e. (7-15) 

This is shown graphically in Figure 7.1 where point B is assumed 
to have a temperature T and vapor pressure e. To approach 
saturation isobarically we must draw a line from B perpendicu- 
lar to the pressure axis. Point C corresponds to a temperature 
where the water vapor pressure is its equilibrium vapor pressure. 
This point thus defines T dew . In this diagram AT = e sw (T) and 
BT = CT dew = e = e sw (T dew ). Obviously, since the equilibrium 
vapor pressure curve is known the relative humidity and dew point 
temperature of any point in the diagram can easily be found 
(r = TB/TA, T dew = OT — BC). If only the temperature is 
known then we need to know T dew in order to estimate the rela- 
tive humidity and the absolute humidity. From the figure it follows 
that neither relative humidity nor the dew point spread (or dew 
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Figure 7.1 

Graphical definition of the 
dew point temperature, 
T dew . Clearly, 
e s „(T d ew) = e(T). 



point depression), T — T dew , are direct measures of the actual water 
vapor amount. Both describe how far we are from the equilibrium 
curve but similar values for r or T — T dew correspond to different e 
and thus to different p v (absolute humidity). Note that, depending 
on where a given point in the (e, T) diagram is, T dew (or T f ) may or 
may not be defined. For example, in Figure 7.2, given point P only 
Tdew can be defined. For point P', T dew = T/ = triple point tem- 
perature. For point P" and in the presence of supercooled water, 
both T and T dew can be defined with T" > T dew . If supercooled 
water is not present, T dew cannot be defined. 

It can be shown that the dew point temperature of a rising parcel 
decreases. If we denote as T dew ,o and T deWjZ the dew point temper- 
atures at the surface and at some higher level z, respectively, we 
have that 

e sw (7d eW: o) 6o P\- 

Csw(Tdew,z) R v p v T z . 

As a parcel rises and expands, both its temperature and vapor 
density decrease (because the vapor is distributed over a larger 


Figure 7.2 

Depending on where a 
given state is in a ( p , T) 
diagram, T dew or Tf may 
or may not be defined. 
For example, for state P 
only T dew can be defined. 



7.2 PROCESSES IN THE ATMOSPHERE 


107 


volume). Thus e 0 > e z and e aw (T deWi0 ) > e sw (T dew>z ). Then from 
equation (6.17), it follows that T deW:Z < T dew0 . 


Estimating T dew 

Since e = e aw (T dew ) the C-C equation can be written as 


or 


de 


dT d( 


l v e 

TD rp2 
1 ' v L dew 


1-ydT 

~e R T 2 ' 

'-'SW w 


Considering Z v ss constant and integrating from T dew (where 
e S w(T dew ) = e) to T (where e sw (T) = e sw ) we obtain 


6 / 

1 '-'SW •'V 

e i? v 


T — T dl 
TT dm 


or 


- ln r = TT 

Tlv 


/ T — T dew \ 
V TT dew / 


or 


T - T dew = —R v TT dew lnr/l v . 


(7.16) 


Equation (7.16) provides T dew from temperature and relative 
humidity. Using Z v « 2.501 x 10 6 J kg -1 and i? v = 461.51 J kg -1 K _1 , 
the above equation becomes 


T - T dew = -1.845 x 10~ 4 TT dew ln r. 


Following the same procedure we can obtain (for l s = 2.8345 J kg 4 ) 
T-T { = -1.628 x 10 _4 TT f ln r. 


Estimating the dew point temperature can be visualized in the fol- 
lowing (p, T) diagram (Figure 7.3) where pressure decreases along 
the y-axis. 

In this diagram a couple of dry adiabatic (solid) and a couple of 
saturation mixing ratio (broken) lines are shown. The dry adiabats 
can be calculated from Poisson’s equation 

ln 6 = ln T + k d ln 1000 — k d ln p 

by setting ln 9 = constant: 

ln p = — ln T + constant. 
k d 
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Figure 7.3 

Illustration of the 
estimation procedure of 
the dew point 
temperature. 



In a (p,T) coordinate system, this is not a straight line. The sat- 
uration mixing ratio lines can be calculated from equations (7.5) 
and (6.17) 

ee sw (T) 


or 


or 


or for some w BW 


6. lie 


w s . 


exp 19.83 — 


5417 \ 
T ) 
5417 


In p ss In h 19.83 — 


= constant 


In + 


B 

T 


where A = In 6.11e/u> sw + 19.83 and B = —5417. In a (p, T) dia- 
gram this is not a straight line either. For this reason both the 
dry adiabats and the saturation mixing ratio lines are depicted as 
nonlinear curves. Since the purpose here is to illustrate the process 
graphically, these curves are realistic but by no means exact. Let 
us now assume that an air parcel is at a pressure level p and has 
a temperature T. Let us further assume that its mixing ratio is 
lOgkg -1 . The parcel is not saturated because at this (jp,T) state 
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the saturation mixing ratio is 15gkg _1 . In order to estimate the 
dew point, we then have to move, as the arrow shows, to the left 
(lower temperatures) (but remain at the same pressure level) until 
we “meet” the lOgkg” 1 saturation mixing ratio line. At that tem- 
perature, the mixing ratio of the parcel is equal to the saturation 
mixing ratio and thus the parcel becomes saturated. That temper- 
ature is then the parcel’s dew point temperature. Note that if from 
(p, T) we follow the dry adiabat 6 1 down to 1000 mb we obtain the 
potential temperature of the parcel, which will be greater than its 
temperature, T. 


7.2.2 Adiabatic isobaric processes - wet-bulb temperatures 

Let us now suppose that our system is unsaturated moist air and 
liquid water and that the mixture has a temperature T and a com- 
bined water mixing ratio w t = ( m v + m w )/m d . Since this is an 
unstable situation the system will tend to approach equilibrium. 
This will be achieved by water evaporating. If we assume that the 
system is closed and the process is taking place adiabatically (i.e. 
the process is also isenthalpic), then the liquid water must draw 
heat for the evaporation from the unsaturated moist air. Since no 
other processes are taking place, this decrease in temperature must 
be related to the change in the mixing ratio of moist air due to 
evaporating liquid water. Note that, since at any time during this 
process the system differs slightly from saturation, the process is 
spontaneous and irreversible. If there is enough liquid water, sat- 
uration will be reached. The temperature at saturation is called 
the isobaric wet-bulb temperature, T w . An important atmospheric 
process described by this type of transformation is cooling by evap- 
oration of rain at a given level. The system consists of a certain air 
mass plus water that evaporates into it from the rain falling through 
it. If instead of water we had ice present, then we would define the 
isobaric ice-bulb temperature, T ; . Note that T w is not the same as 
T dew since the two cooling processes are different. Now, if we follow 
the reverse procedure, we could define the temperature which air 
would achieve if all its vapor were to condense in an adiabatic iso- 
baric process and with the produced water being withdrawn. This 
limiting temperature is called the isobaric equivalent temperature, 
T ei . Note, however, that the definition of this temperature is not a 
“legal” one. The reason is that we assume that the reverse of an 
irreversible process can happen, which is not possible. 


Derivation of the wet-bulb temperature 

Consider a closed heterogeneous system consisting of dry air, water 
vapor, and liquid water. If this system undergoes a change, then 
the change in enthalpy is 
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AH = H f - Hi 

where the subscripts i and / refer to the initial and final states, 
respectively. We can write the above equation in terms of the three 
components 


AH — (H d f + H v f + H w f) — (H di + H vi + H wi ) 


or 


AH = ( m d h df + m v fh v f + m^fh^f) - ( m d h di + m vi h vi + m wi /i wi ). 

(7.17) 

If we recall the definition of l v ( l v = h v — h w ) and consider that the 
total amount of water m t = m v + m w is conserved (i.e. m vi + m WI = 
m v f + m w f = m t ), then the above equation becomes 

AH = m d (h d f - h di ) + m t (h wf - h wi ) + m vf l vf - m vi l vi . (7.18) 

If we now assume that in the range of meteorological interest c pd 
and c w are nearly independent of temperature, then we can write 
(recall equation (4.18)) that 

f h d i C pd (Tf if) ) 

hv.- f h w i Cv,{Tf Ti) . 

In addition, if X) is not greatly different from Tf we can assume 
that 

l =1 t —l 

Then equation (7.18) becomes 


AH = m d c pd (Tf - Ti) + m t c w (T f - T)) + (m v/ - m vi )l v 
or 

AH/m. d = (c pd + — c w ^j (' T f - T t ) + ) l - 

or 

AH/m d = (c pd + w t Cw)(T f - T)) + (w f - w t )l v 
or 

A H/m d = (c p d + w t c w )AT + l v Aw. (7-19) 

Note that w t (the total water mixing ratio) is a different con- 
stant for each system and varies only with the concentration of 
liquid water. As such it does not depend on temperature. Since we 
have also assumed that c pd ,c v ,l v ^ f(T) we can write the above 
equation in differential form: 


dH/m d = (c pd + w t c v )dT + l v dw. 


(7.20) 
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Equation (7.20) is a general expression for the change in enthalpy of 
a closed heterogeneous system of dry air, water vapor, and liquid 
water during any process. If we restrict ourselves to isenthalpic 
processes only, then equation (7.20) becomes 

(c p d + w t c w )(Ti - T f ) = (w f - Wi)l v . (7-21) 

Equation (7.21) indicates that in an adiabatic isobaric process, 
where moist air cools from evaporating liquid water, the change 
in temperature relates to the change in the mixing ratio. It also 
tells us that Tf depends on w t with Tf becoming minimum when 
w t — > 0. Since m t = m v + m w cannot be zero (if it were neither 
the mixing ratio nor temperature would change), this condition 
requires that —> oo. Of course this cannot be true but to a good 
approximation ?rid m t . We define this limiting temperature as 
the wet-bulb temperature, T w , which is determined by setting in 
equation (7.21) T t = T,Tf = T w , Wi = w, Wf = w s w , and w t = 0: 


Cpd(.T T w ) (^Csw w)l v 


or 

T„ + —w sw =T+ —w (7.22) 

Cpd Cpd 

where w s „ is the saturation mixing ratio over liquid water at 
temperature T w . As expected, since w sw > w,T w < T. 


Estimation of the isobaric equivalent temperature 

In order to estimate this temperature (T ei ), the only thing we have 
to do is to set in equation (7.21) 'uy = w t = w (the mixing ratio of 
the moist air), Wf = 0, T) = T, and Tf = T ei . Then we obtain 

T ei = T + — . (7.23) 

Cpd T mc w 

Obviously T ei > T. Also, since for typical values of w (around 
10gkg -1 ) l v /(cp a + wc w ) fn 2500 it follows that T ei > T virt . Neither 
the wet bulb nor the isobaric equivalent temperatures is estimated 
graphically. They are calculated from equations (7.22) and (7.23), 
respectively. 


Relation between wet-bulb and dew point temperatures 

If we employ the general approximation w ~ ee/p we can write 
equation (7.21) as 

c p (Ti - T f ) = ^(e f - ei)l v 
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where c p = c pd + w t c„. Setting T) = T, T f = T w , and e f = e sw (T w ) 
yields 

e» = e sw (T w ) — j^-(T — T w ). 


From the definition of dew point temperature we have that = 
e sw (Tdew)- Then, the above equation becomes 


^'sw(Tdew) 


,(T W )-^(T-T W ). 

t v e 


Since T — T„ > 0 it follows that e sw (T dew ) < e sw (T w ). Then, from 
equation (6.17) it follows that 

Tdew < T w . 

Are you keeping up with all these temperatures? Up to now we have 
discussed five different temperatures: T, temperature; T virt , virtual 
temperature; T dew , dew point temperature; T w , isobaric wet-bulb 
temperature; and T ei , isobaric equivalent temperature. They are 
related as follows: 


r r r T <r'T 1 <'T 7 - 'T ■ 

dew ^ 1 w ^ 1 ^ 1 virt ^ ei • 


Brace yourself. There are more! 


(7.24) 


7.2.3 Adiabatic expansion (or compression) of unsaturated 
moist air 


Recall Poisson’s equation for an adiabatic expansion or compression 
of an ideal gas from p, T to p', V 



which for p' = 1000 mb defined the potential temperature, 9, 


9 = T 



(7.25) 


Poisson’s equation is valid for any ideal gas. As such it is valid for 
moist air provided that the appropriate k is used and that con- 
densation does not take place. We can then define the potential 
temperature of unsaturated moist air as 

/iooo\ fcd(1 -°- 26q) 


p 


(7.26) 
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If we substitute T from equation (7.25) applied to dry air into 
equation (7.26) we obtain (assuming that p A ~ p) 


Sm = Od 



— fcd 0 . 26 q 


or 



(7.27) 


As we have already discussed, in our atmosphere q <C 1 which 
indicates that to a very good approximation 


6„, ss 0 A 


As in the case of virtual temperature, here again we can define the 
virtual potential temperature , 9 V irt , by substituting T virt for T in 
equation (7.25): 


Q rri 

^virt vir 



(7.28) 


This will be the temperature dry air will acquire if it were to expand 
or compress from a level (T virt ,p) to the 1000 mb level. Since T virt > 
T it follows that 9 vilt > 9. 

Note that even though q and w are very small (of the order of 
0.01) we should not jump to the conclusion that all variables or 
parameters of unsaturated moist air are for all practical purposes 
the same as those of dry air. If we consider a typical value of q = 
w = 0.01 we find that differences in virtual temperature of the 
order of 2-3 °C are not uncommon. Such differences are often not 
negligible. On the other hand, as we saw above, differences in moist 
potential temperature are indeed negligible (of the order of 0.1 °C). 
Consequently, from now on the symbol 9 will be used to indicate 
the potential temperature of both dry and unsaturated moist air. 
Similarly to the dry adiabatic lapse rate we may define here the 
moist (unsaturated) adiabatic lapse rate as 


or 


r = 9_ = 9 

m c p c pd (l + 0.87w) 

r " = TT^' sra(1 ~°- 87u ' ) - 


(7.29) 


7.2.4 Reaching saturation by adiabatic ascent 

As an unsaturated moist air parcel rises adiabatically, it cools. Dur- 
ing the ascent (and as long as no condensation occurs) its q or w 
remains the same but the equilibrium vapor pressure decreases. At 
the same time the vapor pressure e also decreases. This can be 
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easily seen by recalling that e = wp/(w + e). As w and e remain 
constant e/p remains constant. Since during the ascent p decreases, 
e must decrease if e/p must remain constant. Accordingly, when 
unsaturated moist air ascends adiabatically, its relative humidity 
(r = e/e s ) increases only if e decreases at a slower rate than e s . 
If this condition is not met clouds will not form in ascent but they 
could still form in descent. However, in our experience, clouds never 
do form as a result of descending air, so this condition is apparently 
satisfied. But what is the mathematics and physics behind it? 

If we log-differentiate equation (7.7) we get 


d In r = d In e — d In e sw . 


(7.30) 


1 — O' 

For moist air, Poisson’s equation tells us that Tp = constant. 

1 — O' 

Since e/p remains constant, it follows that Te = constant. Then, 
again by log-differentiating, we obtain 


d In T = 


7-1 


d In e. 


7 


(7.31) 


Combining (7.30) and (7.31) with the C-C equation yields 


d In r = d In T - J^dT. (7.32) 

7-1 R V T 2 v ' 

The first term on the right-hand side of the above equation accounts 
for the change in relative humidity due to the decrease in e while 
the second accounts for the change in relative humidity due to 
the temperature decrease. As one term is positive and the other 
negative, the net result can be either positive or negative. This 
indicates that adiabatic expansion can lead to an increase or to a 
decrease in r, which in turn implies that the opposite (i.e. adiabatic 
compression) can also lead to a decrease or an increase in r. This 
opens the road to bizarre scenarios where clouds can actually form 
on descent. If we rewrite equation (7.32) as 


dr _ r ( 7 K \ 

dT~ T \ 7 — 1 ~~ R^f) ’ 


we see that if 


7 - 1 > R V T ’ 


(7.33) 


then dr < 0 when dT < 0 and dr > 0 when dT > 0, which would 
indicate that clouds can form on descent. Assuming the typical 
values for 7 ss 1.4, l v = 2.5 x 10 6 Jkg' 1 , and R v = 461.5 J kg -1 K _1 
we find that inequality (7.33) is satisfied only when T > 1550 K. 
Obviously, such a situation does not exist on our planet and clouds 
can only form on ascent. However, all it takes is an exotic planet 
whose atmosphere has a condensable gas with smaller l v and all 
sorts of things may happen. 
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Now that we have settled this issue, we can expect that ascend- 
ing moist unsaturated air will eventually reach a relative humidity 
of 100%. Immediately after that condensation will take place in 
order to sustain equilibrium conditions. The level at which satu- 
ration of ascending unsaturated moist parcels is achieved is called 
the lifting condensation level (LCL). The temperature at this level, 
T LC l, is called the saturation temperature. By integrating equa- 
tion (7.32) from an initial state T, T dew ,r to the final state ?lcl — 
Tdew.LCL, r-LCL = 1 we obtain 


— In r = 


7 

7-1 


In 


2lcl 

T 


U_ ( l 

Rv \7lcl 



(7.34) 


Since r = e/e sw (T) and e = e aw (T dew ) we have from the C-C 
equation 


c.sw (7 (i ew ) 6.11 exp 


U 

273 R v 


R v Td, 


and 


e sw (T) 


6.11 exp 


A 

273 R v 



Dividing the above two equations and taking logs yields 


— In r = 


R v 


/ T — T dew \ 
\ TT dew J 


Accordingly, equation (7.34) becomes 


r r j 1 

dew 

TT 

x L dew 


A In 


r LC L 
~T~ + 




where 



(7.35) 


Equation (7.35) can be solved numerically to provide the saturation 
temperature T LCL from T and T dew . If you would rather not solve 
this equation you may use the approximation provided by Bolton 
(1980) 

Tec l = , 1 lnr + 55, (7.36) 

T— 55 2840 

where T is in K, r = w/w s „, u; sw = e[e sw /(p - e sw )], e sw = 
6.11exp(19.83 - 5417/T), and T 

-^dew — i? v TT dew In r/l v . 

A parcel that rises from an initial level where the temperature is T 
to a level where the temperature is T LCL has expanded. Therefore, 
the vapor density at LCL is smaller than that at the initial level. 
From the ideal gas law it then follows that e(T LCL ) < e(T). Because 
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Tlcl = Td ew ,LCL and e(T) = e aw (T dew ) this inequality is equivalent 
to e sw (T LC L) < e sw (T dew ). Then from equation (6.17) it follows that 
T LC l < T dew . Thus, relation (7.24) can be extended to include yet 
one more temperature: 

Tlcl < T dew < T w < T < T v < T ei . (7.37) 


Estimating the height of LCL 

The height of LCL for a given parcel (zlcl) depends only on its 
temperature and relative humidity. At z LC l the parcel’s tempera- 
ture equals its dew point: T LC l = T dewLC L- As the parcel is lifted 
from a reference level z 0 its temperature decreases according to 

T(z) = T 0 - r d (z - zo) (7.38) 


where T 0 is the temperature of the parcel of z 0 and we have assumed 
that r, n T d . Similarly, its dew point temperature decreases 
according to 

T dew (z) = T dew , 0 - r dev ,(z - Zq). (7.39) 

Here both T d and / dew are considered constant. From equations 
(7.38) and (7.39) we have that 

TlCL = T 0 — /d(^LCL — Zq) 


or since T LCL 


T de w,LCL T dew Q 
T de w,LCL 7 


f dew ( ^LCL Zq) 


^LCL — Zq — 


To — T dew ,o 

r — r 


(7.40) 


In the above equation we know everything except for / dew . In order 
to estimate r de „ we proceed as follows. The equation that defines 
dew point is e = e sw (T dew ). Differentiation with respect to z gives 


or 


de sw (T dew ) de 

dz dz 


^SW (T dew ) dT dew de 

dT dew dz dz 


(7.41) 


Now we know that e = wp/(w + e). If we assume that w remains 
constant until condensation occurs, we have 


de w dp 
dz w + e dz 

or 


de _ edp _ e sw (r dew ) dp 
dz p dz p dz 


(7.42) 
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Combining equations (7.41) and (7.42) gives 

1 de sw (T dew ) dT dew 1 dp 

e sw (T dew ) dT dew dz pdz' 

After considering the C-C equation we can reduce the above to 

l v dl dew 1 dp , j 

dz pdz' 


Since the ascent is assumed to be adiabatic, we have from Poisson’s 
equation that 


or 


7 — 1 

d In T = d In p 

7 

dT 7 — 1 dp 
Tdz 7 pdz 


(7.44) 


From equations (7.43) and (7.44) it follows that 


K dT dew _ 7 dT 

-RvTfew dz 7 -1 Tdz 


7 

(7-i)r 



_ 3 _ 

Cpd, 


or assuming 7 « 7 d 

/ v ^-^dew 9 

R.Tl w dz R d T 

which then gives 


r d , 


= g T dew 
dz el v T 


The value of the right-hand side in the atmosphere ranges from 
about 1.7 to 1.9°Ckm~ 1 . Therefore we can assume that on the 
average T dew decreases with height at the rate of 


T dew = 1.8 °Ckm _1 . 


By dropping the subscript 0 in equation (7.40) it follows that the 
height of LCL is given by the approximate relationship 


ZlCL — z 


T — T 1 

x dew 

8 


(7-45) 


Equations (7.35) or (7.36) and (7.45) give the temperature and 
height (in km) of the level where clouds begin to form from the 
temperature and dew point temperature at some reference level z. 
Reaching saturation via adiabatic ascent and the estimation of T LCL 
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Figure 7.4 

Illustration of the 
estimation procedure of 
the lifting condensation 
level and saturation 
temperature. 



and the height of LCL can be visualized in Figure 7.4. The setup in 
this figure is the same as in Figure 7.3. Since the parcel is initially 
unsaturated, as it ascents, it follows the dry adiabat passing from 
its initial state (p,T) (recall that 9 m ss 9 d ). As it rises (follow the 
arrow) it cools and approaches the w a „ = lOgkg -1 line, which it 
intersects at a pressure level (height) Plcl(-Slcl)- At that point, the 
temperature is T LCL . As the figure illustrates, T LCL < T dew . 


7.2.5 Saturated ascent 

Once saturation is reached, further ascent results in a further 
increase in relative humidity, which means that the vapor pressure 
becomes greater than the equilibrium vapor pressure. At this point, 
the system will return to equilibrium after the “extra” water vapor 
condenses and water droplets (or ice depending on the temperature) 
form on condensation nuclei. From this point on two extreme pos- 
sibilities may be considered. (1) The condensation products remain 
in the parcel. This is a reversible process because if we reverse 
the ascent the products will evaporate. It is also adiabatic as we 
assume that no heat is exchanged between the parcel and the envi- 
ronment. Since it is reversible and adiabatic it is also isentropic. (2) 
All condensation products fall out and the parcel consists always 
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of dry air plus saturated water vapor. This makes the parcel an 
open system. Obviously, this process is neither reversible nor adi- 
abatic and thus it is not isentropic. We will call this process a 
pseudo-adiabatic process. In reality these two extreme cases may 
never exist as some condensation products are likely to remain sus- 
pended in the parcel. However, because of the conditions imposed 
the extreme cases are easier to treat. Nevertheless, the study of the 
extreme cases offers useful insights and approximations for the real 
cases in between. 


Reversible saturated adiabatic process 

Since the process is assumed to be reversible, the heterogeneous 
system of dry air, water vapor, and liquid water (or ice if it is too 
cold) must be at equilibrium at all times. In this case the total 
entropy of the system, S, is conserved: 

S = Sd + S v + S„ = constant. 


If we divide all the terms in the above equation by the mass of dry 
air we arrive at the following expression 

S m v to w 

— — Sd H s v H s w . 

m d TOa rrid 

Recalling that the parcel is saturated {w = w 3W ) and the definition 
of the total water mixing ratio (w t = (m v + m w )/ro d), the above 
equation can be written as 

g 

— = s d + w sv s v + (w t - ■u; sw )s w . (7.46) 

m d 

For a reversible process 


SQ dU pdV 
dS = — = — + ~y~ 

In addition, 


1 

T 




+ ^dV. 

(7.47) 


dS = 




From equations (7.47) and (7.48) it follows that 


and 



1 

T 




(7.48) 


(7.49) 


(7.50) 
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For liquid water and water vapor in equilibrium at temperature T 
and pressure p equation (7.50) gives 

rr S v -S w U v — U w 

1 = h v 

V - V V -V 

r v r w r v r w 

or 

T(S V - S w ) = U v - J7 W + p(V v - V w ). 

If we use equation (6.3) the above equation becomes 


T(S V 


S w ) = H V -H„ = L V 


or 


T(s v 5 W ) — l v . 


Substituting s v from the above equation into equation (7.46) results 
in 


S IvUlsw 

— = Sa + ui t s w H — — 

m d T 


(7.51) 


Now recall that 


s d = c pd hr T - R d In p d + constant [p = p d + e sw (T)] 
and, from equation (5.9), 

S' w = c w In T + constant. 

Accordingly, equation (7.51) becomes 

S l vo 

— = (c„ a + WtCw) In T — R d In p d + vV ^ w + constant. (7.52) 
m. d T 

Since here S and m d are conserved, it follows that 

(Cp d + ic t c w ) In T — R d In p d + ^ v ^ sw = constant'. (7.53) 

As previously, here we can assume that w t , c pd , c* ^ f(T). However, 
u> sw through equations (7.5) and (6.17) is a function of T. If we 
assume that l v is a function of T (and that may be a good idea 
here as the parcel may rise to very cold temperatures), then we can 
write equation (7.53) in differential form as follows: 

(Cp d + w t c w )d In T - R d d In p d + d = 0. (7-54) 

Equations (7.53) and (7.54) describe reversible saturated adiabatic 
processes. Note again that here T and p correspond to saturation. 
If we now define 6' as 

1000 \ ^Acpd + ^tCw) 


9' = T 


Pd 


(7.55) 
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it follows that 


(c pd + w t c w )d In O' = (Cp d + w t c„)d In T — i? d d In p d . (7.56) 


Combining (7.56) and (7.54) gives 


(c pd + w t c„)d In O' = —d 


l-v'W gw 

T 


or 


d In O' = -d 


LyUl SW 

(cpd. T ay c w; )T _ 


Integrating the above equation results in 


O' exp 


l v W SV j 

(Cpd + ui t Cw)T 


= constant. 


(7.57) 


The above equation defines the family of curves that describe 
reversible saturated adiabatic processes. We can evaluate equation 
(7.57) at a state where all the vapor has condensed (u> sw = 0). This 
defines a new temperature, the equivalent potential temperature 0 e 
given by 


0 e = O' exp 


l v w sw 

(Cpd + w t c w )T 


= constant. 


(7.58) 


Although 0 e is defined for saturated air it can also be defined for 
unsaturated air at temperature T, pressure p , and mixing ratio w 
by considering it to be lifted to the LCL. At LCL we have 

( 1000\ ^/Cpd+^tCw) 

) 

Pd ) 


(equation (7.55)) and thus 


— 2lcl 



-Rd/Ccpd+^tCw) 

exp 


*v(T lci > 
(c p d + ^tCw)^LCL 


(7.59) 


Note that in the above two expressions p d = Plcl — e sw (T LCL ). 
Another way to express 0 e for unsaturated air at T, p, w is given by 
Emanuel (1994): 


6L = T 


^ 1000 \ ^d/Ppd+^tCw) 


, — wR w / (c pd +WtCvr) 


Pd ) 

x exp [l v w/ (Cpd + w t c w )T] 


(7.60) 


where r is the relative humidity. Here p d = p — e(T). When the 
initial mixing ratio is zero, both (7.59) and (7.60) give (as expected) 
0 e = 0 , the potential temperature. 
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Recall that the moist potential temperature 0 m (7.26) is 

( 1000 \ Rd / c v d< ' 1 -° 26 ‘ w } 

T) 

Comparison between 0 m (ss 6) and 9' indicates that 9‘ is only 
approximately equal to 6. Given the definition of 8 ' , it follows from 
(7.60) that the equivalent potential temperature is approximately 
equal to the potential temperature a parcel would have if it were 
lifted to very low pressure so that all its water vapor condensed 
(w = 0). 


Pseudo-adiabatic processes 

The description of reversible saturated adiabatic processes depends 
on the value of w t . The problem with this dependence is that 
although w sw can be determined at a given state, ( p,T),w t can- 
not. The liquid water mixing ratio can take on any value and as 
such at a given ( p , T) point there may exist an infinite number of 
reversible saturated adiabats. If we were to describe such a pro- 
cess in the (p, T) domain this would present a great inconvenience. 
This problem can be dealt with by defining a new process called 
a pseudo-adiabatic process. In such a process we do away with liq- 
uid water by assuming that it is removed as soon as it is produced. 
Obviously, our system is now an open system and this process is not 
reversible. However, we can think of the whole process as a two-step 
process. During the first step we have a reversible saturated adia- 
batic process, with condensation of a mass of water (dm v ). During 
the second step the water produced is leaving the system. For the 
first step the entropy may be defined as previously but with the 
liquid water term omitted. Accordingly, we set in equation (7.52) 
w t = w aw and we obtain 


5 


m d 


(c p d+w sv ,c w ) In T R d In p d T 


l v w sw 

T 


+ const ant. 


(7.61) 


Since in this stage S and m d are conserved we can write the above 
equation in differential form 

d[{c p d + w sw c w ) In T] - R d d In p d + d = 0. (7.62) 


After water forms it is removed immediately, thereby decreasing the 
entropy of the system. However, this process has no effect on the 
values of T and p. As such equation (7.62) describes the variation 
of p and T in a pseudo-adiabatic process. Equation (7.62) is very 
similar to equation (7.54) but here w a w depends on temperature 
whereas w t in equation (7.54) does not. This makes equation (7.62) 
very difficult to solve analytically. However, numerical techniques 
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can be employed to define, similarly to the case of the equivalent 
potential temperature, a pseudo-equivalent potential temperature, 
0 e p . As is shown in Bolton (1980), within 0.3 °C, 9 ep is given by 


w{l + 0.81u>) ( -2.54^ 

V Alcl ) _ 

(7.63) 

where (T, p, w) is any state of the parcel saturated or not (w here 
is dimensionless and T, T LCL are in K). 

The pseudo-equivalent potential temperature can be interpreted 
as the actual temperature achieved by the parcel under the fol- 
lowing thermodynamic processes: (a) pseudo-adiabatic expansion 
to zero (or to a very low) pressure level at which we can assume 
that all water vapor has condensed and fallen out, and (b) sub- 
sequent dry adiabatic descent to 1000 mb. Note that according to 
the definition of 9 e (the condensation products do not leave the 
parcel) no similar meaning can be attached to 9 e . Nevertheless, 
the differences in the temperature of a parcel in a reversible sat- 
urated adiabatic and in a pseudo-adiabatic process are not great 
and are often neglected as processes in the atmosphere lie some- 
where in between true reversible saturated adiabatic processes and 
pseudo-adiabatic processes. This can be easily verified if we assume 
that u> sw « w t < c p d (i.e. c pd + w s „c v « c pd ) and p d > e ws (i.e. 
p w p d )- In this case equations (7.54) and (7.62) become identical: 
c pA d In T — R d d In p + d (w avr l v /T) = 0. 

Graphically a saturated ascent can be described with the following 
diagram (Figure 7.5). At LCL the parcel is saturated. Now, as it 
continues to rise, it follows a path (depicted by the broken line) 
dictated by the saturated adiabat passing through (plcl; T lcl ) (in 
practice by the pseudo-adiabat, i.e. from equation (7.62) or (7.63)). 
At some very low pressure level (say, 200 mb), we may assume that 
all vapor has condensed and the condensation products have fallen 
out. At that limit, the pseudo-adiabat approaches a dry adiabat 
(solid line). This dry adiabat corresponds to 9 ep (in other words if 
followed down to 1000 mb, it gives the pseudo-equivalent potential 
temperature) . 


= T 


1000 


— U . ZOW ) 


exp 


7.2.6 A few more temperatures 

The dry adiabatic descent from zero pressure level to the 1000 mb 
level defined the pseudo-equivalent potential temperature, 0 ep . At 
the original (p,T) level it defines the pseudo-equivalent tempera- 
ture, T ep (Figure 7.6). Since the descent is dry adiabatic we have 
that 
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Figure 7.5 

Illustration of the 
estimation procedure of 
the pseudo-equivalent 
potential temperature. 




Since T = 9 (p/1000) i:d it follows that 

T e p = T°-f. (7.64) 

During the process by which T ep is defined, much more water vapor 
condenses than it would have if it condensed during an adiabatic 
isobaric process at an initial level below LCL. Therefore (recall 
the definition of T ei ), if all that vapor were to condense at the 
initial level it would result in a temperature greater than T ei . Thus, 
T ep > T ei , which makes T ep the highest temperature. 

Many more temperatures can be defined but I will spare you this! 
Just keep in mind that at a given state, if the air is unsaturated, a 
dry adiabat can be followed to 1000 mb to define a corresponding 
potential temperature. As such from the state (T ei ,p) the isobaric 
equivalent potential temperature, 9 ei , can be defined. Because T ep > 
T ei it follows that 6 ep > 9 ei . It also follows that 9 ep > 6. 

Similarly, if the air is saturated, a pseudo-adiabat can be fol- 
lowed to either the original pressure level or to the 1000 mb level 
to define yet more temperatures. For example, and recalling that 
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Figure 7.6 

Illustration of the 
estimation procedure of 

T ep ,T W p, and $ W p. 



T w is attained at saturation in an adiabatic process, following 
the pseudo-adiabat corresponding to point (7 l GL . p LCL ) to (T,p) 
level defines the pseudo-wet-bulb temperature, T wp , and continuing 
down to 1000 mb defines the pseudo-wet-bulb potential tempera- 
ture 9 „ p (see Figure 7.6). Obviously, since both 9 ep and 9 wp are 
defined by the same curve they must be related and conserved 
in a pseudo-adiabatic process. 0 wp can be estimated from equa- 
tion (7.63) by setting T = T LCL = 9 wp and w = w' = w a (p = 
1000 mb, T = 0 wp ), i.e. 


9 


ep 


9 wp exp 


w'(l + 0.81-u/) 




(7.65) 


Similarly, following the pseudo-adiabat passing through point 
(T w ,p) to 1000 mb defines the isobaric wet-bulb potential 
temperature 0 W . Here we should note that T wp and T w (or 0 wp 
and 6 W ) are not exactly the same. T w and 9 W correspond to a 
true saturated adiabatic process whereas T wp and 0 wp correspond 
to a pseudo-adiabatic process. Since water remains around in a 
true saturated adiabatic process it turns out that T wp < T w and 

^wp ^ U w • 

It is interesting to note here that if at the initial point ( p , T ) (see 
Figure 7.6), the amount of vapor content ( w ) is higher, then the dry 
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adiabatic ascent will last a shorter time and therefore the LCL will 
be found at higher pressure and temperature. The pseudo-adiabat 
passing through this point will be of a greater temperature than 
before. This means that higher 9„ p (or 9 ep ) values would correspond 
to higher moisture (in accordance with 7.63). 

Here is then the final breakdown of all the temperatures and 
potential temperatures and their relationships: 

T LC l < T dew < T wp < T w < T < T virt < T ei < T ep 

(7.66) 

0wp < 0w < 0{~ 0m) < 0virt < 0e i < 0 ep . 


7.2.7 Saturated adiabatic lapse rate 


Since condensation of water vapor releases heat during a reversible 
saturated adiabatic ascent, the cooling of the ascending air slows 
down. As such the corresponding lapse rate should be smaller than 
that of the dry or moist unsaturated adiabatic ascent. 

We could derive an expression for this saturated lapse rate, T s , by 
starting with equation (7.54) 


dT dpd 

C P rri did 

T p d 


Us, 

~Y~ 


= o 


where c p = c pd + u> t c w . Multiplying all terms of the above equation 
by T/c p gives 

dT - ^-d{p - e sw ) + — d(l v w svf ) - dT = 0, 

CpPd Cp CpJ 

where we have considered that p = p d + e sw . The second term, using 
the C-C equation in the form de svl /dT = e sv ,l v /R v T 2 and the fact 
that i? d e sw /i? v p d = w s „, can be written as 


didT w sv/ l v 

dp H iprdT. 

CpPd Cpd 

Substituting this into the previous equation, dividing all terms 
by dz, assuming that dp « dpd, and using the hydrostatic 
approximation, leads to 


dT 

dz 


9_ 

c p 


1 d 


c p dz 


(l v w svl ) . 


(7.67) 


Note that since in the above equation c p ^ c pd , the term —g/c p is 
not exactly the dry adiabatic lapse rate. By performing the differ- 
entiation on the right-hand side we can write equation (7.67) as 
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or 


or 


dT _ g w a „ dl v U dw a „ 

dz c p c p dz Cp dz 

dT g w sw dl v dT l v dw a „ 

dz Cp c p dT dz c p dz 

A w aw dl v \ dT _ g l v dw aw 

\ Cp dT J dz c p c v dz 


Using equation (6.12) the above equation becomes 




dT 

dz 


g l v dw sw 

c p c p dz 


(7.68) 


Typical values for c w ,c pv ,c p , and w a „ are 4220 J kg' 1 K _1 , 
1850 J kg' 1 K' 1 , 1050 4 kg' 1 K' 1 , and 0.01, respectively. Such values 
make the term rc sw (c w — c pv )/c p of the order of 0.02. As such, it can 
be neglected to arrive at the following approximation 

dT g l v dw sv , . . 

~dz~~~^ v ~ 7p~d^' ( ' > 


Now all we have to do is to find an expression for dw aw /dz. Log- 
differentiating equation (7.5) leads to 


1 dw a „ _ _ 1 dp 1 de sw 

w sw dz p dz e sw dz 


1 dp 1 de a „ dT 
p dz e sw dT dz 

where p is the total pressure of the parcel or the surroundings. If 
we again use the hydrostatic approximation and the C-C equation 
the above equation becomes 


1 dw avr _ 4 dT 

w a „ dz RT R V T 2 dz 

Combining equations (7.70) and (7.69) leads to 

r _ _dT^ _ 1 + l v w sv /RT 

dz c p 1 + l 2 w av , / CpFRT 2 


(7.70) 


(7.71) 


where c p = c p a + WtC w and R is the gas constant for the surround- 
ings. In other thermodynamics books you may find some slight 
differences in the formula giving T s . Most of the differences are due 
to certain assumptions and approximations. For typical values of 
saturation mixing ratio and temperature in the lower troposphere, 
T a ranges around 5°Ckm~ 1 , which is half the dry adiabatic rate. 
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7.3 Other processes of interest 
7.3.1 Adiabatic isobaric mixing 

Let us consider two air masses of temperatures 7\ and T 2 and 
humidities qi and q 2 , which are mixing at the same pressure. This 
corresponds to a horizontal mixing. If we assume the mixing to be 
an adiabatic process, then it is also isenthalpic. If we assume that 
no condensation takes place, then we can write for the change in 
enthalpy of the whole system 


AH = niiAhi + m 2 Ah 2 = 0 (7.72) 


where 


Ah x = c pl AT = c pl (T - Tl) 

Ah 2 = c p2 AT = c p2 (T — T 2 ) 

and T is the final temperature of the mixture. Considering equation 
(7.11), equation (7.72) becomes 

mic pd (l + 0.87qi)(T — Tl) + m 2 c p d(l + 0.87q 2 )(T — T 2 ) = 0. 


Solving for T we get 


(toiTi + m 2 T 2 ) + 0.87(miqiTi + ?n 2 q 2 T 2 ) 
to + 0.87 to v 


(7.73) 


In the above expression m = m i + m 2 is the total mass and 
m. v = TOi qi + m 2 q 2 is the total water vapor mass. Since there is no 
condensation, m v remains constant. If the final specific humidity is 
q, then the water vapor mass is qm. It follows that 


q = 


TOiqi + m 2 q 2 

TO. 


(7.74) 


i.e. the final q is the weighted average of qi and q 2 . Using equation 
(7.74) we can rewrite equation (7.73) as 

(tolTl + to 2 T 2 ) + 0.87(TO 1 q 1 T 1 + TO 2 q 2 T 2 ) 

T= m(l + 0.87,) ' ( " 5) 


If we neglect all vapor terms the above equation becomes 


T 


TolTl + to 2 T 2 
m 


(7.76) 


indicating that, approximately, the temperature of the mixture is 
also the weighted average of the initial temperatures. One can 



7.3 OTHER PROCESSES OF INTEREST 


129 


show that the potential temperature and the vapor pressure of the 
mixture will also, approximately, be the weighted averages of the 
initial values (see problem 7.1): 


m 1 9 1 + m 2 6 2 

u ~ 

(7.77) 

m 


m 1 ei + m 2 e 2 

(7.78) 

e ~ . 

m 


Equations (7.76) and (7.78) indicate that in a vapor pressure- 
temperature diagram (Figure 7.7) if the two masses correspond 
to points A 1 and A 2 then the mixture will correspond to the point 
A which will lie on the straight line connecting A 1 and A 2 . From 
the similarity of triangles AA' 1 A 1 and A 2 A' 2 A it follows that point 
A will be found at a distance such that A±A/AA 2 = AA' 1 /A 2 A' 2 = 
(T) —T)/(T — T 2 ) = m 2 /m 1 (here equation (7.76) was used). The 
thin solid lines in Figure 7.3 represent the equal relative humidity 
curves corresponding to points A ll A 2l and A. Because of their cur- 
vature it can easily be seen that the final relative humidity will be 
greater than the weighted average. This can easily be visualized if 
we consider the simple case where the air masses have the same rel- 
ative humidity, i.e. points A 1 and A 2 lie on the same equal relative 
humidity curve (Figure 7.8). Then point A will lie to the left of 
the curve indicating that the final relative humidity is greater than 
the weighted average, (mr + nnr ) /2 m = r. We conclude that iso- 
baric adiabatic mixing increases the relative humidity. This is a very 
interesting result because it implies that two air masses that are ini- 
tially unsaturated (but close to saturation) may result in a mixture 
that is supersaturated. In this case condensation will occur. This 
mechanism can produce mixing fogs and other important processes 
in the atmosphere. 


Figure 7.7 

When two initially 
unsaturated air masses 
(Ai at ri and A 2 at r 2 ) 
mix isobarically, the final 
product (A) may be 
saturated. 
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Figure 7.8 

When two saturated 
(r = 1) air masses (Ai 
and A 2 ) mix isobarically 
the final product ( A ) is 
supersaturated. 


r= 1 



T 


7.3.2 Vertical mixing 


Let us now look at vertical mixing occurring because of turbulent 
and/or convective processes. Such a process is much more com- 
plicated than the previous one because now the vertical variation 
of T and p come into the picture. However, we can simplify the 
treatment of vertical mixing by (1) considering two isolated air 
masses m i and m 2 at pressure levels pi and p 2 ( p\ > p 2 ) and tem- 
peratures X) and T 2 , which are taken to a pressure level p where 
they mix isobarically, and (2) redistributing the mixture in the 
layer A p = p 1 — p 2 . The excursions to level p can be considered as 
adiabatic ascent and descent (during which the specific humidities 
qi and q 2 remain constant) resulting in new temperatures for the 
air masses: 



Once at level p the air masses mix isobarically. From what we 
discussed in the previous section it follows that in this case 


T 

q 


m\T[ + m 2 X \ 
m 

TOiqi + m 2 q 2 
m 


(7.79) 


and 


8 


m 1 6 1 + m 2 0 2 
m 


where 9 1 and 0 2 are the initial potential temperature values, 
which remain constant during the adiabatic ascent and descent. 
The subsequent redistribution of the mixture in the layer A p 
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will preserve 9 assuming that it consists of adiabatic ascents and 
descents. It follows that when the layer is well mixed, 9 will not 
vary with height. We can generalize equation (7.79) by considering 
mixing of n air masses and employing the general expression for 
the weighted average. For example, 


- f. m 9dm 
9= ^ , 

TO 

r_qdm 

H 5 

m 

f m wdm 
w = — . 

TO 

Using the hydrostatic approximation we can express the above 
equations as 

_^ = Ip x P dz = JrT xdp 
fo pdz _pi - p 2 

where x is 9 , or q, or w. 

The temperature variation with pressure can then be expressed 
according to 

T = 9 (-^—Y ■ (7.80) 

VlOOO/ v ; 

If in the initial distribution (before mixing) there is a level 
( p,T < T), where r = r, then this level is defined as the mixing 
condensation level. 


7.3.3 Freezing inside a cloud 

As air rises and cools and clouds begin to form, it is often the case 
that the developing cloud would reach the freezing level. From that 
point on, freezing of cloud droplets (or rain drops) would occur 
and water vapor and ice will co-exist. Does this freezing affect the 
further development of the cloud? We can get an answer to this 
question by making the following simple assumptions: 

1) Initially the air is saturated with respect to water at 
temperature T. 

2) As water freezes latent heat is released into the environment. 

3) Once water freezes, water vapor is supersaturated with respect 
to ice. Depotition is then initiated releasing more latent heat. 
This continues until the vapor pressure becomes equal to the 
saturation vapor pressure over ice. 

4) The whole process occurs at a constant pressure level and it is 
adiabatic. Since it is adiabatic and isobaric, it is also isenthalpic. 
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Accordingly, we can divide this process into three “treatable” steps: 


Step 1: Water freezes at a constant temperature T. In this case we 
have a change of phase at constant pressure and the change 
in enthalpy is (see Chapter 6) 

A // 1 ^Qp=constant -bf / f w 

where ?n w is the amount of water that freezes. The minus sign 
is due to the fact that heat is released. Dividing the above 
equation by the mass of dry air yields 

A M 

i = —l{W w (T) (7.81) 

m d 

Step 2: Vapor condenses on ice at a constant temperature T until the 
water vapor pressure become the saturation vapor pressure 
over ice at a temperature T' > T. Like in Step 1, the change 
in enthalpy during this step will be 

AH 2 = —l s (m„ - mi) 

or 

A H 

2 = -l B [w v (T) - Wi (r)] (7.82) 

TYld 

Step 3: The whole system is warmed from T to T'. The change in 
enthalpy during this term will be (recall 4.18) 


A H 3 = CpAT 


or 

A H 3 

= c pd AT 7.83 

m d 

where it is assumed that C p /m d « c pd . Since the whole 
process is isenthalpic, it follows that 

-l { w w (T) - l s [w w (T) - Wi(T')} + c pd A T = 0 


or 


AT = 


IfW^jT) + l s [w w (T) - Wj (T')] 

Cpd 


or using (7.5) 


AT = 


^fCsw(T) + d s [e aw (T) — e si (T')] 


pCpd 
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In the above equation and in the range of typical conditions in 
our atmosphere, the right-hand side is positive. As such the effect 
of freezing is to increase the temperature inside the developing 
cloud. Typically, this temperature increase is of the order of 3 K, a 
significant warming during cloud development. 


Examples 

(7. 1 ) The initial state of a parcel of air is p = 1 000 mb, T = 30 °C, w = 
1 4 g kg- 1 . Find (a) r, (b) T virt> (c) T dew , (d) T w , (e) T ei , (0 T LCL , 
(g) ^■LCL* (h) 9 e p, (i) / ep , (j) $ w p> 0<) 9 v i rt • 


(a) 


(b) 

(c) 

(d) 


(e) 

(f) 

(g) 

(h) 

(i) 

(j) 

(k) 


From equation (6.17) we estimate that e sw (T) = 
43.1 mb. From equation (7.3) we then find that 
e(T) ss 22.0 mb. Thus, r ss 0.51, w s „ « 28gkg _1 . 

T virt = (1 + 0.61q)T = [l + 0.61 (^) 


T= 305.55K. 


From equation (7.16) we obtain T dew = 291.75 K. 
From equation (7.22) we have 


T„ = T + —(w - w sw ) 

Cpd 

where l v (T = 303 K) = 2.43 x 10 6 Jkg- 1 . Now recall 
that ie sw in the above formula is the saturation mix- 
ing ratio at T w . Making use of the approximations in 
equations (7.5) and (6.17) we have 


T =T 

-L -or -L 


l v e 
C/xi P 


e(T) - 6.11 exp 19.83 - 


5417 


Numerical solution of the above equation gives T w « 
295.4K. 

From equation (7.23) for c w = 4218 J kg -1 we find that 
T ei w 335 K. 

From equation (7.36) it follows that T LC l ~ 289 K. 
From equation (7.45) and assuming that the initial 
level is z = 0 we find that z L CL = 1.4 km. 

Using equation (7.63) with T = 303 K, p = 1000 mb, 
w = 0.014 g kg - 1 5 Tlcl = 289 K, we obtain that 9 ep « 
345 K. 

From equation (7.64) it follows (since T = 9) that 
T e p = 9 ep = 345 K. 

Considering that w' = w s (1000 mb, T = 9„ p ) = 
(e/p)6.11 exp[19.83 — (5417/0 wp )], numerical solution 
of equation (7.65) yields 9 wp w 295 K. 

From equation (7.28) and since p = 1000 mb it follows 
that 9 virt = T vilt = 305.55 K. 
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Note that all estimated temperatures and potential temper- 
atures are consistent with equation (7.66). 


(7.2) The temperature of saturated air at 1000 mb is 10 °C as it begins 
to drop through radiation loss. As a result, some water vapor 
condenses and radiation fog forms. During the formation of fog, 
12 x 10 3 J kg -1 , of heat is lost to the surroundings. Find the final 
temperature and the decrease in vapor pressure. 

The process is non-adiabatic and thus non-isenthalpic. 
From equation (7.20) we have 


dh = c p dT + l v dw 


where c p = c p a + wye*, and we have assumed that dH/m a ~ 
dh. During the process the air remains in a state of satura- 
tion, thus w corresponds to saturation. Also, the process 
is isobaric, so dh = Sq and we can write the above 
equation as 

Sq = c p dT + l v dw s „. (7.84) 

Using the approximate relationship tu sw ss ee sw /p, we find 
that dw sw ~ (e/p)de sw or, using the C-C equation, 

dw ~ = lwh dT - < 7 - 85 > 

By substituting (7.82) into (7.81) we obtain 


Sq 



el 2 e \ 

c(r V C 'SW \ 

pR v T V 


dT. 


In the above equation, c p is to a good approximation 
independent of temperature. Also, for small temperature 
changes, l v is nearly independent of temperature. The sat- 
uration vapor pressure, however, is dependent on T (via 
equation (6.17)). This makes the integration of the above 
equation very difficult. One way to overcome this difficulty is 
to assume that e sw /T 2 is nearly constant for small tempera- 
ture differences and then treat the differentials as differences. 
In this case, for T = 10 °C, e sw = 12.16 mb, Z V (10°C) = 
2.4774 x 10 6 Jkg-\(5g = -12 x 10 3 Jkg" 1 , and c p « c pd = 
1005 J kg -1 K _1 , we find that AT = —5.3 °C. Thus, the final 
temperature is 4.7 °C. Note that this result can be improved 
if we do not assume that c p ss c p d . However, this will require 
the value of w t which is not known. We could of course 
assume that w t « ui sw (10°C) and consider c p s=s c pd + u> sw c w . 
In this case, we find that the difference in the result is 
small (~ 0.066 °C). In order to find the decrease in vapor 
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pressure we combine equation (7.81) with dw a „ « (e/p)de sv , 
to arrive at 

Sq = c„dT + — de 3W . 

P 

Then incorporating the C-C equation gives 


Here again if we treat the differentials as differences and use 
the values mentioned above we find that Ae sw = —4.326 mb. 

(7.3) Two equal air masses, both at 1000 mb, mix adiabatically. 
Their initial temperatures and mixing ratios are Ti = 293.8 K, 
Wi = 16.3 g kg -1 , T 2 = 266.4 K, and w 2 = 1.3 g kg -1 . First 
show that the mixture is supersaturated. Then find the final 
temperature, the final mixing ratio, the amount of liquid water 
content per unit mass of air, and the amount of liquid water per 
unit volume of air. 

As the two air masses mix isobarically, the initial tempera- 
ture of the mixture will be (equation (7.75)) T = 280.2 K. 
Since qi = 16.0386 g kg -1 and q 2 = 1.2983 g kg” 1 we find 
that, for the mixture, q = 8.6685 g kg” 1 (equation (7.74)), 
w = q/(l — q) = 8.7443gkg _1 , and e = wp/(e + w) = 
13.86 mb (using the approximate equations (7.76) and (7.78) 
results in T = 280. 7K and e = 14.1 mb). 

For T = 280.2 K equation (6.17) yields e sw = 10.05 mb. This 
result indicates that initially the mixture is supersaturated 
(e > e sw ). Thus, immediately after mixing the excess vapor 
will condense so that the mixture reaches the equilibrium 
state (T',e'). This step, graphically shown in Figure 7.9, 
will result in a decrease in vapor pressure and an increase 
in temperature (remember the mixing is adiabatic, so the 



Figure 7.9 

Diagram for example 7.3. 
When two initially 
saturated air masses mix 
isobarically, the final 
product is supersaturated 
(T, e). In order for the 
mixture to reach the 
equilibrium state (T',e'), 
the excess vapor will 
condense. The amount of 
vapor that condenses is 
proportional to the 
difference e — e'. 


e 



(h, e,) 


r= 1 


T 
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heat gain due to condensation remains in the mixture) . The 
increase in temperature is denoted by the line connecting 
points (T, e) and (T', e) and the decrease in water vapor pres- 
sure by the line connecting points (T',e') and ( T',e ). Since 
the process from ( T , e) to (T', e') is isobaric and adiabatic it 
is isenthalpic. As such, equation (7.21) applies: 

(c p d + w t c„)(T - T') = (w' - w)l v 


or 

(c pd + w t c w )(T - T') = (e' - e)l v ^ 
or 


(c pd + w t c w )(T — T') 


6.11 exp 


^19.83 




where, assuming that the mixture is a closed system, 
w t = w = 8.7443 g kg -1 . Numerical solution of the above 
equation, with c w = 4218Jkg _1 K _1 and l v = 2.4774 x 
10 6 Jkg _1 , gives T' « 282.8 K. This corresponds to e' = 
e sw (T') = 12 mb and hence to w' = 7.55gkg _1 . 

The amount of liquid water produced per unit mass of air 
is given by the difference q — q' « w — w' « |(e — e') « 
1.16 g kg -1 . This amount must be equal to the decrease of 
water vapor per unit mass, —dm v , which is proportional to 
the decrease in water vapor pressure, — cie sw . From the ideal 
gas law we know that the mass of water vapor per unit vol- 
ume at temperature T is e sw /i? v T and its variation with 
temperature is 


de s „/dT = de sv ,/R v T — (e sw / R v T 2 )dT. 


Using the C-C equation we can verify that the second term 
on the right-hand side of the above equation is much smaller 
than the first term. As such, the amount of water vapor 
lost is approximately equal to de av ,/R v T. It follows that the 
amount of liquid water per volume is equal to de a „/R v T , or 
equal to 

e — e! , 

fl-jr » 1-44 gm 3 . 


(7.4) The following problem (Salby, 1996) is a classic example, demon- 
strating the effect of the downslope wind called Chinook wind 
(Chinook in American Indian means “snow-eater”). Moist air 
from the Pacific is forced to rise over the continental divide. On 
the west side, the surface lies at 800 mb where the temperature 
is 293 K and the mixing ratio is 15gkg _1 . If the summit lies at 
600 mb and assuming that any condensation products fall out, 
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find the surface temperature and relative humidity at a place on 
the leeward side at 830 mb. 

For T = 293 K equation (6.17) gives e sw = 23.4 mb. From 
equation (7.4) we then find that the saturation mixing ratio 
at state (' T,p,w ): (293 K, 800 mb, 15gkg _1 ) is 18.74 g kg -1 . 
Since w < w s , the parcel is initially unsaturated. As the 
air begins to rise over the mountain it cools and its relative 
humidity increases. From equation (7.36) we find that, for 
r = 15/18.74 = 0.8, T LCL = 288. 7K. 

The ascent to the level where the air first becomes saturated 
is described by Poisson’s equation 



where k = fc d (l — 0.26 re) = 0.2849. It follows that Plcl ~ 
760 mb. At this point the air is saturated. As it continues 
to rise condensation takes place and all condensation prod- 
ucts fall out. As such, the ascent above p LCL is described 
by a pseudo-adiabatic process. Therefore, 9 ep from 760 to 
600 mb is conserved. From equation (7.63) we find, for T = 
Tlcl = 288.7 K, p = p LCL = 760 mb, and w = w SW (T LCL ) = 
15gkg _1 , that 9 ep « 359 K. Evaluating the same equation 
at p' = 600 mb gives 


359 = T 


/ ^QQQ ^ 0.285(1— 0.28u/) 

V v 

x exp 


'3376 


w’(l + 0.81«/) ( - 2.54 

V Ao 


(7.86) 


where T', p', and w' is the state of the air at the summit. In 
the above equation we know p' but T' and w' are unknowns. 
However, w' can be expressed as a function of T' alone by 
combining equations (7.5) and (6.17): 


w' = w S w(T') 


ee sw (T') 


6. lie 

exp 

V' 


^19.83 


5417 \ 
~T r ) ' 


Thus, in effect, equation (7.83) has only one unknown, T' . 
Solving this equation numerically gives T' ss 279 K. It fol- 
lows that at the summit w' « 9.6 g kg -1 . From that point 
on the air descends adiabatically with the mixing ratio 
remaining constant (9.6 g kg - ). From Poisson’s equation we 
then find that at 830 mb, T ~ 306 K. At this temperature 
ie sw ~ 38.5 g kg -1 and at p = 830 mb the relative humid- 
ity is r = 0.25. Thus the air on the leeward side is warmer 
and drier than that on the west side. Note that we assume 
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that the parcel from 760 to 600 mb remains saturated with 
respect to liquid water. This is a safe assumption because 
T lcl = 288.7 K and the vertical distance between 760 and 
600 mb in a standard atmosphere is about 1.5 km. So, we 
would expect that the temperature at the summit would be 
above 0°C. 


(7.5) What is the maximum rate of precipitation when saturated air 
rises? 

Consider a saturated parcel of air rising. For simplicity, we 
assume that this parcel is a unit mass. As the parcel rises 
water vapor condenses. If we assume that all condensed 
products (rain or snow) fall out we can estimate the upper 
limit for the rate of precipitation. In order to achieve this we 
consider equation (7.54) that describes reversible saturated 
adiabatic processes, which we simplify by assuming that l v = 
constant and c p d + u^c* ~ c p 

c pd d In T - R d d In p d + l v d ( ) = 0 

or 

c pd d In T - R d d In p d + dT = 0 

or by multiplying all terms by T 

c pd dT — TR d d In p d + l v dw svl — v ^ sw = o 
or 

dT — — d In p d + l v dw sw = 0 
P 

or using the hydrostatic approximation 

^c pd - dT + gdz + l v dw sv , = 0. 

For typical conditions in the atmosphere l v w s „/T is very 
small compared to c p d- We can thus neglect it to arrive at 

c pd dT + gdz + l v dw sv , = 0 



or 


or 


l v dw aw 

dz 


dT 

~ Cpd d ^ +9 


l v dw a „ _ dT g 
c pd dz dz c pd 
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Now, recall that P d = — and that — ff- = P s . Thus, the 
above equation becomes 

K dw sw 

, = 

Cpd dz 

or 

dw sw Cpti ( / d r„) 

dz l v 

or 

dw sw dt __ c pd (r d - r s ) 
dt dz l v 

where is the rate of ascent or speed of the parcel, U . It 
follows that 

dw svr Cpd ( I i] f 's ) jj 

dt l v 

The above equation gives the mass of condensed water per 
unit mass of air and unit time. If we multiply it by the mass 
of air, m, we can obtain an expression for the mass of water 
condensed and precipitated per unit time 


or 


or 


P= c ? d i r j — ^ Um 


p = Cpd rs) u P v 


P = Cpd (/d ^ UpAAz 


where A is the area and A z the thickness of the rising 
volume. Thus, per unit area and unit time 


p = c pd (Id r *) UpAz ' 

lv 


In the above equation, P is given (in MKS units) in 
kgm -2 s -1 . In practice, the amount of condensed water and 
precipitation per unit area is expressed by the depth of the 
collected water in millimeters (mm). Considering that the 
density is about 10 3 kgm -3 , then lkgrn -2 is equivalent to 
1 mm of precipitated water. If, in addition, we consider that 
lh = 3600s, then we can express the above equation in 
mmh -1 , which is the unit for the rate of precipitation used 
in practice 


P = 3600 


c P d(^d ~ r s ) 

L 


U pAz. 


I 
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Problems 

(7.1) Show that when two air masses mix adiabatically and isobar- 
ically without condensation taking place, the final potential 
temperature and final vapor pressure are the weighted 
averages of the initial values (make necessary assumptions) . 

(7.2) Show that 

Plcl _ P 
^SW (Tlcl) ^SW (Tdew) 

Then use equation (6.17) to derive a relationship between 
Tlcl, 7dew,P, & n( l Plcl- 

(7.3) An air mass has a temperature of 20 °C at 970 mb pres- 
sure with a mixing ratio of 5 g kg - . After some time this 
mass has acquired a temperature of 10 °C and a pressure 
of 750 mb. Assume that no condensation or mixing between 
the air mass and the environment takes place and calcu- 
late the initial and final values of vapor pressure, relative 
humidity, and dew point temperature. (7.74 mb, 5.98 mb, 
0.33, 0.49, 276.5 K, 272.8 K) 

(7.4) Show that evaporation of dm grams of water into ?n d grams 
of dry air (under constant temperature) requires an amount 
of heat absorbed given by 

6Q = m d l v dw. 

(7.5) A Chinook wind blowing at 800 mb has a temperature of 
38 °C and a mixing ratio of 4gkg~\ Is this the same as the 
air at 1000 mb on the windward side of the mountains having 
a temperature of 294.5 K and a mixing ratio of 10gkg -1 or 
is it the same as the air at 800 mb with 278 K and 5gkg _1 ? 
(It is the same as the air at 1000 mb) 

(7.6) Show that if an isothermal layer of temperature T 0 con- 
tained between isobars pi and p 2 mixes vertically, the final 
temperature distribution is given by 

T(p) = [T 0 / (1 - k)} (p{~ k - p\~ k ) p k /( Pl - p 2 ). 

(7.7) What will be the lowest possible temperature to which very 
dry air at T = 40 °C can be cooled by evaporation? Consider 
p = 1000 mb and l v constant and equal to 2.4 x 10 6 Jkg -1 . 
(~288K) 

(7.8) When it is cold outside our breath creates mixing clouds. Is 
it possible to create a mixing cloud when the temperature of 
our breath is less than the air temperature? Assume that the 
temperature of breath is 30 °C. (Hint: Plot e sw as a function 
of T using equation (6.17) and make assumptions about this 
curve for high temperatures. Then use the theory discussed 
in section 7.3.1.) 
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(7.9) On a summer day in Milwaukee the maximum temperature 
is predicted to reach 35 °C with dew point around 28 °C. 
What will be the lowest temperature to which the air can 
be cooled by evaporation in this case? Does the answer make 
sense? Elaborate on the difference between this answer and 
the answer in problem 7.7. Consider l v constant and equal 
to 2.4 x 10 6 Jkg- 1 . (~303K) 

(7.10) Moist air moves over land where the temperature is 20 °C 
and the relative humidity is 70%. Through contact with the 
ground the air cools. At what temperature will fog form? 
Consider l v constant and equal to 2.45 x 10 6 J kg -1 . (287.5 K) 

(7.11) The main cabin in an airplane is maintained at a pressure of 
900 mb and at a temperature of 25 °C. If it suddenly decom- 
presses adiabatically to a pressure of 500 mb, what should 
the relative humidity be to avoid cloud formation in the 
cabin? (6.29%) 

(7.12) Outside air has a temperature of — 15 °C and a relative 
humidity of 0.6. If the air indoors is 25 °C, what is the rela- 
tive humidity inside the room (assuming that the air inside 
is only heated, not humidified)? If the room has a volume 
of 100 m 3 how much water vapor must be supplied so that 
the relative humidity rises to 50% and what will then be the 
total mixing ratio? If the temperature change due to evap- 
orating water is neglected, what amount of heat must be 
added for this to happen? Assume l v constant and equal to 
2.5 x 10 6 Jkg” 1 . (3.6%, 1.08kg, 8.97 g kg -1 , 2.7 x 10® J) 

(7.13) Which of the following quantities are conserved in (a) 
reversible adiabatic unsaturated transformation, (b) rever- 
sible adiabatic saturated transformation, (c) adiabatic 
unsaturated isobaric transformation, (d) adiabatic satu- 
rated isobaric transformation? 

w, q, w t 7 ^|t 7 ^7 ^SW 7 ®sw j T , Plcl , Tlcl , T w , 0 , $ v irt? h- 

(7.14) A refrigerator of volume 2 m 3 is initially filled with air at 
temperature 303 K and with a relative humidity of 50%. 
At what temperature condensation will form on the walls? 
Assuming that the desired final temperature is 275 K, how 
much water vapor must condense to achieve the final tem- 
perature? What is the total amount of heat that will be 
given away to the surroundings during the whole process? 
Assume l v = 2.5 x 10 6 Jkg _1 , c w = 4218 J kg _1 K _1 , p d = 
1.293 kgm" 3 . (~ 19°C,19.8g, 117368 J) 

(7.15) If the temperature inside a room is 25 °C and the tem- 
perature outside is — 10 °C calculate the maximum rela- 
tive humidity that can be accommodated inside the room 
without windows fogging when (a) the interior part of the 
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windows is not thermally isolated (i.e. heat is transferred 
freely) from the exterior part, and (b) it is thermally iso- 
lated. Assume l v = 2.5 x 10 6 Jkg -1 , c w = 4218 Jkg _1 K _1 . 
(0.32,1.0) 

(7.16) If the wet-bulb temperature of a parcel at 1000 mb is 15.6 °C 
and its mixing ratio is bgkg^ 1 , find the relative humidity 
of the parcel after it is lifted dry adiabatically to 900 mb. 
Assume l v = 2.47 x 10 6 Jkg” 1 . (36.14%) 

(7.17) Raindrops at 10 °C evaporate into air of temperature 20 °C. 
If the saturated mixing ratio of air at 10 °C is 8 g kg -1 , what 
is the mixing ratio of the environmental air? Make assump- 
tions and consider l v = 2.47 x 10 6 Jkg -1 . (3.93 g kg -1 ) 

(7.18) If water vapor comprises 1% of the volume of air (i.e. 
it accounts for 1% of the molecules) what is the virtual 
temperature correction? (T virt = 1.0038 T) 

(7.19) If all the water vapor in a column of air extending from the 
surface to the troposphere were to condense and fall to 
the surface, it would occupy the lower d w of the height of 
the column. In this case, we can write that 




w 



p v dz 


where p w and p v are the densities of liquid water and water 
vapor, respectively. Show that the upper limit for d„ is given 

by 

, TqCso 


where r is the lapse rate, which is assumed to be constant, 
T 0 is the surface temperature, and e s0 is the saturation vapor 
pressure with respect to water at T 0 . Get an estimate of d„ 
by assuming realistic values for T 0 and IT What do you 
think of the estimated d w ? Is it realistic? [d w is often called 
precipitable water.] 

(7.20) Show that the mixing ratio of a parcel of air is given exactly 


1 + (1 — r)^f ' 


(7.21) In Meteorology Book I, Aristotle says: “Both frost and dew 
are found when the sky is clear and there is no wind.” What 
do you think? Is he right or wrong? Explain why. 



CHAPTER EIGHT 


Vertical stability in the 
atmosphere 


Over large scales the atmosphere is very nearly hydrostatic, which 
means that the pressure gradient force balances the force of gravity. 
Because the net force in the vertical is zero, the atmosphere over 
those scales exhibits slow vertical motions of constant speed, i.e. 
there is no upward acceleration of air. However, on smaller scales 
hydrostatic equilibrium may be invalidated. In this case, acceler- 
ated motion leads to convection which embraces many phenomena 
observed in the atmosphere, from the structure of planetary bound- 
ary layers to the dynamics of hurricanes. In this chapter we will 
investigate the conditions that allow such accelerating motions. 
In particular, we will examine the fate of a parcel at equilibrium 
with its environment when it is subjected to a small perturbation. 
We will make the following assumptions: (1) the environment is 
in hydrostatic equilibrium; (2) the parcel does not mix with its 
surroundings; (3) the parcel’s movement does not disturb the envi- 
ronment; (4) the process is adiabatic; and (5) at a given level the 
pressure of the environment and the pressure of the parcel are equal. 
In problem 8.1 you will be asked to elaborate on the validity of the 
above assumptions. 


8.1 The equation of motion for a parcel 

Because the environment is in hydrostatic equilibrium the following 
equation holds: 


For the parcel, the above equation is not applicable because as the 
parcel is displaced (upwards or downwards) it has some acceleration 
(< d 2 z/dt 2 ). Assuming that gravity and the pressure gradient force 
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are acting on the parcel, Newton’s second law dictates that per unit 
volume (primes identify the parcel) 


,d 2 z 

P W = ~ P9 


dp 

dz 


or 


z = -g-a 


,dp^ 

dz 


Because of assumption (5) above, dp' /dz = dp/dz. Thus, the above 
equation can be written as 


z = —g — a 



or 


or 



(8.1) 


The right-hand side gives the force per unit mass acting on the 
parcel due to the combination of gravity and pressure gradient and 
is called the buoyancy of the parcel. Using the definition of virtual 
temperature and assumption (5), we can write the ideal gas law for 
both environment and parcel as: 


P pR d?vi r t 


and 


P = 


p'RdK h 


( 8 . 2 ) 


The reason for adopting these expressions rather than p = pRT and 
p = p'R'T' is that this way we deal with T virt and T( irt rather than 
R,R',T, and T'. This makes the analysis more straightforward. 
Combining equations (8.1) and (8.2) yields 


z = g 


T virt - T vir 

T virt 


(8.3) 


What we are interested in here is to investigate the effect on the 
motion, described by equation (8.3), of small displacements (z <C 1) 
from an original equilibrium level. If, for simplicity, we take this 
level to be the z = 0 level where the temperature is T virtj0 and 
express both T virt and TL rt in terms of Taylor’s series we get 


T virt = T vi 


dT vilt t 1 d 2 T vilt 2 


dz 


~ Z+ 2~ dz* 


z' + ■■ 


and 


T' = T ■ 

virt vr 


dT! ri 


virt.O 


dz 


1 d 2 V 

_ _ ^ virt r 2 

2 dz 2 
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If we now neglect terms of order higher than the first and define the 
environmental virtual temperature lapse rate as —dT vilt /dz = irt 
and the parcel’s virtual temperature lapse rate as —dT^ irt /dz = _T' irt 
we arrive at 


z = 


g{r v irt 


^virt.O 



(8.4) 


In the above equation we can manipulate the term 1/ (T viltfi — r vilt z) 
to arrive at 


1 


T v irt.O 


- r 

' vi r 


1 1 


Tvirt.O 


1 - 


^virt 2 
^virt,0 


1 


T v irt,0 



Then equation (8.4) becomes 


because <C 1. 

^virt.O 


z = 


g 


T v i r t,0 




r 1 . 

* virt 


)z 


or by eliminating terms involving z 1 


z = 


(^v irt - -Cirt) 


1 virt.O 


or 


z + 


g 


Tv irt,0 


(r'. 

V vii 


Tviit) Z = 0. 


(8.5) 


8.2 Stability analysis and conditions 

The solution of the differential equation (8.5) depends on the 
constants. Three possibilities exist. 

(i) r' irt - r virt > o 

In this case equation (8.5) takes the form 2 + \ 2 z = 0 and has 
the solution 

z(t) = A sin A t + B cos A t 

where the oscillatory components are characterized by 



This is the so-called Brunt-Vaisala frequency. Since we 
assumed that the initial level is z = 0, it follows that in this 
case B = 0 and as such z(t) = A sin A t. This indicates that 
the parcel will oscillate in time about its original position with 
a period r = This represents a stable case as the parcel 

does not leave the original level. 
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( 2 ) r' irt - r virt < o 

In this case equation (8.5) takes the form z — A 2 z = 0 and has 
the solution 

z(t) = Ae xt + Be~ xt 

where now 

A = \l ' (-^virt — -^virt) > 0- 

V virt,0 

Since at t = 0, 3(0) = 0 it follows that A+B = 0. This indicates 
that A = — B ^ 0 (the possibility of A = B = 0 leads to 
the trivial solution z(t) = 0 which we ignore). Since A / 0 
it then follows that as t — > oo, a parcel’s displacement grows 
exponentially. Note that since for t — > oo, dz/dt = A\e xt , the 
parcel’s motion is an accelerating motion. This is the unstable 
situation where the parcel leaves the original level and never 
returns. 

(3) r' irt - r virt = o 

In this case equation (8.5) becomes 3 = 0 and has the linear 
solution 

z(t) = At T B 

indicating that the displacement grows linearly with time. The 
parcel’s motion is now one of a constant speed ( dz/dt = A). 
This is the neutral case where the parcel leaves the original 
level and never returns. Thus, the only difference between neu- 
tral and unstable situations is whether or not the motion is an 
accelerating one. 

The physical meaning of the above is that as long as the parcel’s 
lapse rate is greater than that of the environment and the parcel is 
forced upwards (downwards), it will become colder (warmer) than 
the environment and it will sink (rise) back to the original level. 
If the parcel’s lapse rate is less than that of the environment and 
the parcel is forced upwards (downwards) it will become warmer 
(colder) than the environment and it will continue to rise (sink) 
thus moving away from the original level. The difference between 
the two lapse rates is the net force acting on the parcel. If the net 
force is zero the motion is one of a constant speed. Otherwise the 
net force acts in the direction of the initial impulse if /T rt — T virt < 0 
or in the opposite direction if JT rt — T virt > 0. 

Let us now examine in detail what happens when a parcel is dis- 
placed in a layer, which is characterized by a virtual temperature 
lapse rate T virt . From equation (7.9) we have that for an unsatu- 
rated parcel TL rt = (1 + 0.61w/)T', where w' is the initial mixing 
ratio, which because of assumption (2) remains constant. Thus, 

dT' dT' 

= (l + 0 . 61 u/)^— 
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Km = (i + o.6i«/)r m 

or using equation (7.29) 


r; irt = (i + o.6i«/)(i - o.87 w')r d 


or 

Kin » (i - o.26w')r d . 

The term 0.26«/ in the above relation is rather small and therefore 
we may neglect it. Thus, for unsaturated parcels 

K* » ^d. (8-6) 


Note that for the environment, w is not constant but a function of 
z. As such, for the environment we have 


or 


dT v irt 
dz 


(l + 0.61u;)^+0.61T^ 
az az 


r • 

J vir 


(l + 0.61w)r-0.61T 


die 

dz 


(8.7) 


The second term on the right-hand side of equation (8.7) is not 
negligible and therefore we cannot arrive at a relation similar to 
equation (8.6) for T virt and T (the environmental lapse rate). This is 
why we have to use the virtual temperature in our stability analysis 
rather than the actual temperature. If T virt were approximately 
equal to T, then we could use T rather than T virt in our stability 
analysis. Because this is not the case, it is more appropriate to use 
virtual temperatures. 

We can then state the condition for stability for an unsaturated 
parcel in a layer with T virt as follows: 

if r v irt > r d layer is unstable, 

if r vilt = r d layer is neutral, (8.8) 

if r v irt < r d layer is stable. 


For saturated parcels the mixing ratio decreases with height. 
Proceeding as above we have 


or 


T’ m = (1 + 0.61w')V 


dT'. 


dz 


dT' 

(1 + 0.61u;') — +0.61T' 
dz 


dw' 

dz 


or 


r; irt = (i + o.6iu;')r s -o.6iT'^. 



148 


8 VERTICAL STABILITY IN THE ATMOSPHERE 


Figure 8. 1 

Relative position of 
fundamental lines and 
regions of stability. 



In this case the second term of the right-hand side is much smaller 
than the first term. Thus, we can approximate the above equation as 


Accordingly, the condition for stability for a saturated parcel in a 
layer with T virt can be stated as 

if T v irt > r a layer is unstable, 

if .T virt = r a layer is neutral, (8-9) 

if -T virt < r a layer is stable. 

Since T d = 9.8°Ckm” 1 and r a ss 5°Ckm _1 , conditions (8.8) and 
(8.9) can be combined in one as follows: 

if -T virt > r d layer is absolutely unstable, 

if r a < r virt < r d layer is conditionally unstable, (8.10) 

if r vilt < r a layer is absolutely stable. 


The word “absolutely” indicates that the stability criterion holds 
for any type (saturated or unsaturated) of parcel. The term “condi- 
tionally unstable” means that the layer is stable for displacement of 
unsaturated parcels and unstable for saturated parcels (Figure 8.1). 
Now recall that for unsaturated environments equation (7.28) 
applies: 


^virt. — T vi 


1000 V 


- ' 


Logarithmic differentiation of the above equation gives 


1 d9 v irt _ 1 dr virt _ Aq dp 
T virt dz p dz 


6>virt dz 
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or 


or 


or 


1 dd vi r t _ 

9 V irt dz 

1 dd v irt 


1 p __ I _ 

'fi 1 virt I 

-L virt P 






virt 


1 


_ = ^ _9_ 

-'vitl dz virt d virt V ^pd 

1 d0 virt 1 


flvirt ^2 


T, 


(dd — d^virt) . 


( 8 . 11 ) 


virt 


Combining (8.8) with (8.11) yields an alternative way to express 
stability conditions for unsaturated parcels: 


if 


if 


if 


d9 v 


dz 

dO v irt 

dz 

d@ v irt 

dz 


> 0 layer is stable, 

= 0 layer is neutral, 
< 0 layer is unstable. 


(8.12) 


For a saturated parcel similar conditions can be applied if we substi- 
tute $ virt by 9 e (or in practice 9 ep ) which is invariant along saturated 


adiabats (or pseudoadiabats) 



if 

d9 e 

d? >0 

layer is stable, 


if 

~~7~- — — 0 

dz 

layer is neutral, 

(8.13) 

if 

^ <° 
dz 

layer is unstable. 



The above set of criteria applies to the stability as a parcel rises 
in a motionless layer. If this is not true (for example, when entire 
layers are lifted or lowered), then the stability of the layer may be 
affected. Here is why it happens. 

Let us consider a stable layer in hydrostatic equilibrium whose bot- 
tom has a pressure pi and whose top has a pressure p 2 - Let us 
further assume that the difference in pressure, A p, between top 
and bottom remains constant during the lifting or sinking. If the 
whole process is an unsaturated adiabatic process, then equation 
(8.11) applies. Since in this process 0 virt is conserved then the dif- 
ference in 0 virt between top and bottom is conserved. It follows 


that 


1 dd v i r t 

9 V irt dz 


1 

9 virt 


d9v 


dp 


-(- P9 ) 


Tp (dd ~ d'Virt) 

virt 
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or 


or 


1 d0 v irt Rd / ^ \ , , 

— = (1 d - i virt) = constant 

0 vilt dp pg 

Rd{Rd — F virt ) = constant x g x p. 


During lifting pressure decreases, and according to the above equa- 
tion r v irt — > r d . Since the layer is initially stable this tendency will 
cause the layer’s lapse rate to approach the region of instability 
(recall Figure 8.1). During sinking the opposite is true. As such, 
lifting of a layer tends always to destabilize the layer and sinking 
of a layer tends always to stabilize the layer. 

If the lifting process is strong enough to cause the layer to saturate, 
then the picture is quite different. In this case the stability changes 
depend on the way saturation is reached, which in turn depends on 
the vertical distribution of moisture and the vertical structure of 
temperature among other factors. 

For example, let us consider an initially stable unsaturated layer 
at level 1 in which the temperature profile is BT and in which 
dd e /dp < 0 (i.e. the bottom of the layer has more moisture than 
the top of the layer). As Figure 8.2 illustrates when the layer is 
lifted, the bottom of the layer (B) and the top of the layer (T) 


Figure 8.2 

Graphical representation 
of convective instability. 
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will cool at the dry adiabatic rate (solid lines). However, the bot- 
tom will reach saturation sooner than the top. Thus, at level 2 the 
bottom will start cooling at the moist adiabatic rate (broken line) 
whereas the top will continue cooling at the dry adiabatic rate. 
Because of that, at position 3, independently of whether or not the 
top has reached saturation, the temperature profile will be rotated 
to the left in relation to the initial profile at level 1. This coun- 
terclockwise rotation causes an initially stable situation to become 
less stable (see Figure 8.1). Thus, when dd e /dp < 0 the initially 
stable unsaturated layer is potentially or convectively unstable. 

If in the same layer dd e /dz > 0 (i.e. the top has more moisture than 
the bottom) , then the top will reach saturation first and it will begin 
to cool at the moist rate while the bottom will continue cooling at 
the dry rate. In this case the initial temperature profile BT turns 
clockwise (Figure 8.3), which means that the initially stable layer 
will tend to become even more stable. Thus, when dd e /dz > 0 
the layer is potentially or convectively stable. If in the same layer 
dO e /dz = 0 (i.e. all levels of the layer have the same moisture), then 
both top and bottom will reach saturation along the same moist 
adiabat. In this case the initially stable layer is called potentiality 
or convectively neutral. 


Figure 8.3 

Graphical representation 
of convective stability. 
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If the whole layer becomes saturated, the stability criteria for the 
layer translate to the following 

dd e 

if — — > 0 the saturated layer will be stable with 
respect to saturated parcel processes, 

dO 

if — ^ = 0 the saturated layer will be neutral with 
respect to saturated parcel processes, 

dO e 

if — — < 0 the saturated layer will be unstable with 

respect to saturated parcel processes. (8-14) 

where dO e /dz refers to the initial profile of 9 e (or in practice 9 ep ) 
in the layer, which initially is not saturated. In the atmosphere, 
a potentially unstable layer tends to form cumulus-type clouds 
and often convective precipitation. On the other hand, a poten- 
tially stable layer produces stratiform clouds with light (if any) 
precipitation. 


8.3 Other factors affecting stability 

In our initial assumptions we accepted that the parcel does not mix 
with the environment, that the process is adiabatic, and that the 
environment is not disturbed (i.e. there are no compensating ver- 
tical motions in the environment as a parcel or a layer rises in it). 
When these assumptions are violated they affect the temperature 
lapse rate of the parcel and so the stability criteria derived previ- 
ously do not apply. A detailed treatment of these subjects is beyond 
the scope of this book. Interested readers can consult Iribarne and 
Godson (1973). 


Examples 

(8.1) The table below provides pressure, temperature, and dew point 
temperature measurements at various pressure levels. (I) Inves- 
tigate the stability conditions of parcels rising in each layer. (2) If 
each layer were lifted en masse until it became saturated would 
it become convectively stable, neutral, or unstable? (Assume that 
ly is independent of T and equal to 2.45 X 1 0 6 J kg 1 .) 

First, we need to establish whether or not at a given level 
the parcels are saturated. The temperature is at all levels 
greater than the dew point temperature except at 850 mb. 
Accordingly, for the stability of the parcels, conditions 
(8.12) apply for all layers except for the layer 850-800 mb. 
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p 

(mb) 

T 

(°C) 

T"dew 

(°C) 

1000 

30 

23 

950 

27 

21 

900 

23 

20 

850 

20 

20 

800 

18 

10 

750 

15 

5 

700 

10 

2 

650 

5 

0 


For the layer 850-800 mb, conditions (8.13) apply. For the 
convective stability of each layer conditions (8.14) apply. 
Thus, in order to answer the questions we need to calcu- 
late 0 virt and 0 ep at each level. In order to estimate 0 virt 
and 0 ep (recall equations (7.28) and (7.63)) we need to 
find at each level w and T LCL . The mixing ratio can be 
estimated by combining equations (7.3), (7.7), (7.16), and 
(6.17). T lcl can be estimated from equation (7.36). After 
plenty of calculations we arrive at the following table; here e s 
is in mb. 


p T T dew 


1000 303 296 

Layer 1 

950 300 294 

Layer 2 

900 296 293 

Layer 3 

850 293 293 

Layer 4 

800 291 283 

Layer 5 

750 288 278 

Layer 6 

700 283 275 

Layer 7 

650 278 273 


w 

T virt 

$virt 

dep 

r 

e s 

Tlcl 

0.0182 

306.4 

306.4 

357.4 

0.66 

43.1 

294.3 

0.0169 

303.0 

307.5 

355.1 

0.70 

36.0 

292.6 

0.0166 

299.0 

308.1 

355.0 

0.83 

28.2 

292.3 

0.0176 

296.1 

310.2 

360.2 

1.00 

23.4 

293.0 

0.0099 

292.8 

312.1 

340.1 

0.61 

20.6 

281.3 

0.0149 

290.6 

315.5 

361.7 

0.51 

17.0 

276.0 

0.0063 

284.1 

314.6 

333.3 

0.58 

12.2 

273.4 

0.0058 

279.0 

315.6 

332.9 

0.70 

8.6 

272.0 


(1) According to the results in the above table and condi- 
tions (8.12) and (8.13), we find that for parcels rising 
in each of the layers, layer 1 is stable (d0 viTt /dz > 0), 
layer 2 is stable (d0 vilt /dz > 0), layer 3 is stable 
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(dO viIt /dz > 0), layer 4 is unstable ( dd ep /dz < 0), 
layer 5 is stable ( dd vilt /dz > 0), layer 6 is unstable 
(dO vilt /dz < 0), and layer 7 is stable ( dO viIt /dz > 0). 

(2) Similarly, we find that layer 1 is convectively unstable 
( dd ep /dz < 0), layer 2 is convectively neutral ( d0 ep / 
dz « 0), layer 3 is convectively stable ( dO ep /dz > 0), 
layer 4 is convectively unstable, ( d9 ep /dz < 0), layer 
5 is convectively stable (dd ep /dz > 0), layer 6 is 
convectively unstable ( d9 ep /dz < 0), and layer 7 is 
approximately convectively neutral ( d0 ep /dz « 0). 

(8.2) In an atmospheric layer the virtual temperature is constant and 
equal to 10 °C. If an initial fluctuation causes the parcel to rise 
adiabatically, calculate the energy per unit mass that has to be 
given to the parcel in order for it to rise I km above the bottom 
of the layer before it begins to sink. 

When a parcel rises in the atmosphere, a certain amount of 
work is performed by or against the buoyancy force depend- 
ing on whether the motion is along or against the direction 
of the buoyancy force. If the buoyancy force is directed 
downwards (negative buoyancy) a certain amount of work is 
done against buoyancy and if the buoyancy force is directed 
upwards (positive buoyancy) a certain amount of work is 
done by the buoyancy force. In the above problem since the 
layer is isothermal the layer is stable. Thus, the buoyancy on 
the parcel is negative and a parcel forced to rise will have to 
return to the initial level. But before it begins its return, it 
will reach a maximum height which depends on how strong 
the initial impulse was. Recall from Chapter 4 that 


6W = Fdz = madz = mzdz. 


Thus, the work done, W, when a parcel is forced to rise from 
a level i to a level / is 


[ f 

W= mzdz. 


(8.15) 


Using equation (8.3), the above equation becomes 


W = I mg 


TUAz) - T vilt (z) 
T viIt (z) 


dz 


(8.16) 


or 


r f v (z) f f 

W = gm / TT^Adz - gm / dz 
Ji T vilt {z ) Ji 
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or, assuming that the ascent is dry adiabatic, 

w f* ( T virt,o - r d z) 
w = ( J m / — Tf — j—, dz - gm(f - i). 

Ji virt v^y 

Since T vilt (z) = T virtj0 and assuming that the bottom of the 
layer is at 2 = 0 (i.e. i = 0 m and / = 1000 m), we find that 
the work done per unit mass is equal to —169.7 Jkg -1 , i.e. 
work is done against the buoyancy force. Since the ascent is 
assumed to be adiabatic, SQ = 0. Then from the first law 
we obtain 

du = 169.7 J kg -1 . 


This is the amount of energy per unit mass that must be 
given to the parcel for the above process to take place. 

Note that using the hydrostatic approximation, equation 
(8.16) can be expressed as 


or using 


or 



^virt ~~ T v ir 

T virt 


dp 


P pdddT-v irt 

W = —R d m J f (IN rt — T virt ) ^ 

W = - R d m ( (Tl rt - T virt )d In p. 
J i 


(8.17) 


In a (T — Inp) diagram the above equation is proportional 
to the area enclosed between the vertical profile of the 
environment and the vertical profile of the parcel’s virtual 
temperature. Equation (8.16) or (8.17) gives the maximum 
work done by the buoyancy and when divided by the total 
mass defines (remember we are always dealing with adia- 
batic processes, i.e. 5Q — 0) what is called the convective 
available potential energy (CAPE) and the convective inhi- 
bition energy (CINE). More on CAPE and CINE and their 
applications will be discussed in the next chapter. 


Problems 

(8.1) Elaborate on the five assumptions made in derivation of the 
stability conditions of a parcel that is subject to a small 
displacement from equilibrium. 
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(8.2) If the virtual temperature profile in a layer is T virt {z) = 
T V i r t,oa/(a + z), where a > 0 and T virtj0 > add, derive the 
vertical profile of the virtual potential temperature. What 
is the condition for stability for unsaturated parcels in this 
layer? (Unstable for 2 < z c , stable for z > z c , and neutral 
for z = z c , where z c = y/(aT viltfi / r d ) - a) 

(8.3) If the profile of virtual potential temperature in a layer is 

^virt — 

2 = eai - — ’ 

<Adrt,0 fl 2 

what is the condition for stability in the layer for unsaturated 
parcels? (Stable when z > ai ln(ai/a 2 ), unstable when z < 
ai ln(oi/o 2 )) 

(8.4) An atmospheric layer is isothermal. A dry parcel is sub- 
jected to an upward displacement and begins to oscillate 
about its original level. Plot the oscillation period as a 
function of the temperature of the layer. What do you 
observe? 

(8.5) If we define the geopotential <j> according to dc/> = gdz with 
((<(0) = 0 show that under hydrostatic conditions 


Az = (— In— 'j T virt , 

V g PzJ 

where Az is the thickness of a layer bounded by p\ and 
p 2 ( Pi > P 2 ) and T V i r t is the mean virtual temperature 
of the layer. Using the above relationship show that if 
a layer is lifted en masse while the mean virtual tem- 
perature remains constant, then the upper level pressure 
change, dp 2l is related to the lower level change, dpi, via the 
equation 


dpi _ dp 2 
Pi Pi 


(8.6) Show that for a constant virtual temperature lapse rate 
atmosphere (T virt (z) = T virt - T virt z) 


P = Po 



T vi 


-| s/fld r virt 


T v i r t, 0 


(8.7) Show that for an atmosphere where T vilt (z) = T vir t,o 


P = P 0 e 9 2 /fldT v irt,0 


(8.8) From the following data determine the stability and convec- 
tive stability of each layer. 
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p ( mb ) 

T virt (° C ) 

Tdew(°C) 

1000 

30.0 

22.0 

950 

25.0 

21.0 

900 

18.5 

18.0 

850 

16.5 

15.0 

800 

20.0 

10.0 

750 

10.0 

5.0 

700 

- 5.0 

- 10.0 

650 

- 10.0 

- 15.0 

600 

- 20.0 

- 30.0 


(8.9) In an unstable layer of air extending from the ground to the 
pressure level of 900 mb, the virtual temperature decreases 
with height at a rate of 25°Ckm -1 . A parcel of air at the 
surface is given an initial velocity of 1ms -1 . If the virtual 
temperature at the surface is 7°C, and assuming that the 
parcel is and remains dry, find its position and velocity after 
lmin. (81 m, 2.1 ms -1 ) 

(8.10) Sea-breezes and lake-breezes are established when the air 
over the land becomes warmer than the air over the water. 
The warmer air over land rises thus destroying the pressure 
stratification and resulting in horizontal pressure gradients 
between land and water. The height of the breeze is the 
altitude where the horizontal pressure gradient vanishes. 
Assuming that the vertical virtual temperature profiles over 
land and water are isothermal and the same, show that this 
height is given by 

R d ln(p w /pi) 

av 2V irtiW T virtl ' 

where the subscripts 1 and w denote the conditions over land 
and water, respectively. 

(8.11) Consider problems 8.6, 8.7, and 8.10. How would the results 
change if these problems were stated in terms of the tem- 
perature profiles rather than virtual temperature profiles? 




CHAPTER NINE 


Thermodynamic diagrams 


Chapters 1 to 6 laid down the fundamental physical and mathemat- 
ical concepts pertaining to thermodynamics. While we always kept 
the discussion close to our atmosphere, it was not until Chapter 7 
that, in a more applied mood, we presented how these concepts 
are applied to yield quantities useful to atmospheric processes. 
However, even if we understand all the mathematics, we still need 
an efficient way to present and visualize thermodynamic processes 
in the atmosphere. Thermodynamic diagrams can do that. Up to 
this point we made an effort to visualize the thermodynamic pro- 
cess using a (p, V) or a (p, T) diagram. However, such diagrams, 
while simple, may not be very convenient when it comes to their 
utilization. 

Since the purpose of a diagram is efficiently and clearly to display 
processes and estimate thermodynamic quantities, the following are 
very desirable in a thermodynamic diagram: (1) for every cyclic 
process the area should be proportional to work done or energy 
(area-equivalent transformations), (2) lines should be straight (easy 
to use), and (3) the angle between adiabats and isotherms should 
be as large as possible (easy to distinguish). The (p, V) diagram 
satisfies the first condition ( pda = dw), but the angle between 
isotherms and adiabats is not very large (Figure 4.5(a)). Because 
of this, while it is used for illustration purposes, this diagram is not 
used in practice. 


9.1 Conditions for area-equivalent transformations 

When we are constructing a new diagram, in effect, we go from the 
x = a, y = —p domain to a new domain characterized by two new 
coordinates, say u and w. These two coordinate systems are shown 
in Figure 9.1. If we want the first condition to be satisfied then an 
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Figure 9.1 

Area-equivalent 

transformations. 


y=-p 


w 





x = 


a 


u 


area dA = \ dx x dy | in the x, y domain is mapped onto an area 

dA' = | du x dw | in the u, w domain, such that dA = (dA' (x here 
is the cross product between two vectors and ( > 0 is a constant). 
Since any pair of points x and y are transformed to u and w we 
have that x = x(u.w) and y = j/(u,w). We can thus write 


~r* dx dx 
dx= — du +7— aw 
au aw 


Then 


— + dy — > dy ~r> 
dy= -A- du +7— dw 
au aw 


dA = 


/ dx dy 

dx dy\ 

V du dw 

<9w du J 


du x dw 


where we have used the facts that 


(9.1) 


du x du = dw x dw= 0 


and 


du x dw = — dw x du . 
In equation (9.1) | du x dw | = dA' . Thus, 

dA = (dA' 


(9.2) 


where 


dx 

dx 

9u 

dw 

dy 

dy 

du 

dw 


is the Jacobian of the coordinate transformation. It follows that 
the condition for area-equivalence is that the Jacobian should be a 
constant (if J = 1, then this is an equal-area transformation). The 
above approach tells us whether or not the two new coordinates 
provide an area-equivalent transformation, but both u and w must 


9.2 EXAMPLES OF THERMODYNAMIC DIAGRAMS 


161 


be specified. If only one is specified, then the following approach 
provides the other one. Since dA = (dA' it follows that 


j) ydx = <j> —pda = ( wdu 
or 

( pda + (wdu) = 0. 

The above condition is true if pda + (wdu = dz is an exact 
differential. Thus, z = /(a, u) and 


dz 


( dz 


dz(a,u) = ( fa ) da+ ( — ) du. 


\du 


It follows that the sufficient conditions for an area-equivalent 
transformation are 

' dz N 


p = 


da 


and 


/ f)? 


From the above two equations we have 

d 2 z 


and 


Thus, if 


c 


dp 

du 


/ dw 

l9a 


dadu 

d 2 z 

dadu 


( dw\ (dp 


\ da J u 


du 


(9.3) 


then the areas will be equivalent (again here if £ = 1 the areas will 
be equal). From equation (9.3) if we specify u we can determine 
what w should be in order for the transformation to be an area- 
equivalent transformation. 


9.2 Examples of thermodynamic diagrams 
9.2.1 The tephigram 

If we suppose that u = T then 

f dp\ = ( dp\ 

\&r) a - 
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Using the ideal gas law, the above leads to 

dp\ _ R 
du) a' 

From equation (9.3) we then have 

c 


/dw\ _ R 

V da. 


or 


c 


/ dw 

l da 


da = — da 
a 


or after integrating 


£w = R In a + /(T) 


(9.4) 


where f(T) is the constant of integration which we can choose. By 
manipulating Poisson’s equation we can write 


T 

~e 


RT \ k 
1000a ) 


or 

In a = — (In 0 — In T) + In T + In R. — In 1000 
k 

or 

R In a = c p In 8 + f(T). 

If we combine the above equation with equation (9.4) and choose 
f(T) = — f'{T ) and £ = c p , we arrive at an area-equivalent diagram 
with coordinates 


w = In 9 (9-5) 

u = T. 

Obviously, in this diagram the dry adiabats are straight lines and 
perpendicular to the isotherms (Figure 9.2). The equation for the 
isobars can be derived from Poisson’s equation for p = constant: 

In 9 = In T + constant. 

In a coordinate system (In 9, T) the above equation describes loga- 
rithmic curves. As expected, the work done on a unit mass of air in 
a reversible cycle (—q) is proportional to the cycle: —q = — §Tds = 
—c p Td In 9 = —c p A teph where A teph is the area of the cycle in the 
tephigram* (positive if described counterclockwise). If the air is dry, 

Cp Cpd. • 


* The word “tephigram” originated from the coordinates of this diagram, which 
are T and entropy. Entropy was originally denoted as cf> (phi), but is now denoted 
by S. 
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Figure 9.2 

Isobars, isotherms, and 
dry adiabats on a 
tephigram. Also shown is 
the w s = lOgkg -1 
saturation mixing ratio 
line (solid bold line) and 
the pseudo-adiabat 
6 e = 40 °C (broken line). 



9.2.2 The emagram 

Let us assume again that u = T. By taking the logarithm of the 
equation of state we get 

In a = — In p + In R + In T. 

Combining this equation with equation (9.4) yields 

w = -C _1 i?lnp+C _1 [i? In R + R In T + f(T)}. 

If we now choose f(T ) = —R In R — Rln T and ( = R, we arrive 
at the area-equivalent diagram with coordinates 

w = — In p (9-6) 

u = T. 

In this diagram (which is called an emagram - Energy per unit 
MAss diaGRAM) the isobars are parallel straight lines and per- 
pendicular to the isotherms. The dry adiabats are again derived 
from Poisson’s equation 

In 9 = In T + k d In 1000 — fc d In p 


— lnp=— — In T + constant. 
k d 


or 
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Figure 9.3 

Same as Figure 9.2 but for 
an emagram. 



Temperature (°C) 


Since the coordinates are In p and T, the above equation is not a 
straight line. In the emagram the dry adiabats are slightly curved 
lines (Figure 9.3). Their angle with the isotherms is about 45°. 
Here, again, the work done on a unit mass of air in a reversible 
cycle (— § pda ) is proportional to the area of the cycle: 


pda = j) \—d{pa) + adp] 
= <j> [— RdT + adp) 
= — y RdT + <j> 


adp 


= ® adp 


RT— 

P 

Td { — In p) 


= -RA e 


where H em is the area in the emagram enclosed by the cycle, which 
is positive if described counterclockwise. If the air is dry, R = R d . 
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9.2.3 The skew emagram (skew (T-ln p) diagram) 

Let us now assume that u = — In p. Then 


or 


or 


Integrating for p 


c 


/ <9w 


dp 

d In p 


c 


/ <9w 

V^a 


= ~P 


In p 




1 da - —pda. 


da 


In p 

constant yields 


C w = —pa + F(ln p) 


or 

£w = —RT + F(ln p). 

We may choose F(ln p) = — £ In p where £ is a constant. Then 
Cw = —RT — £ In p. 

If we choose £ = I?, we arrive at a new diagram with coordinates 
w = — T — £,/ R. h i p 
u = — In p. 


Since the sign of an area in a diagram involves only the direction 
in which a cycle is carried out, the new coordinate system can be 
written as 

w = T + 4 In p (9.7) 

R 

u = — In p. 


This represents a correction (in fact a rotation) to the emagram 
which provides larger angles between any adiabats and isotherms 
(Figure 9.4). An isotherm in this diagram has the equation 


w = constant + 


A 

R 


In p 


or 


w = constant — — u. 

R 


In the (u, w) domain this equation represents a straight line with a 
slope equal to £/R. Thus, we may choose a £ for which the isotherms 
are perpendicular to the dry adiabats. Note that the equation for 
dry adiabats is not linear. From Poisson’s equation we have that 
In 9 = In T + k d In 1000 — k A In p or In p = (l/fc d ) In T+ constant. 
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Figure 9.4 

Same as figure 9.2 but for 
a skew (T— In p ) diagram. 



While here In p is a coordinate, In T is not and so the dry adiabats 
are not exactly straight lines. Here again it is easy to show that 
the work done in a cyclic process is proportional to the area of the 
cycle. The skew emagram and tephigram are the most commonly 
used diagrams. Other diagrams have been proposed (see problems 
9. 1-9. 3), but their practical use is limited. A detailed skew T-ln p 
diagram is supplied at the end of the book (Figure A.l). 

For a complete description of a thermodynamic process in a dia- 
gram, in addition to isotherms, isobars, and adiabats, we also need 
to define the constant mixing ratio lines and the saturated adiabats 
(see also Chapter 7). From the definition of mixing ratio we have 


p = e aw (T) + 


ee sw (T) 


or 


P 


ee aw (T) 


(9.8) 


For a constant u> sw the above equation defines a family of curves, 
which in the (T, In p) domain are approximately (but not quite) 
straight lines: 


P 


6. lie 


w s . 


exp 19.83 — 


5417 \ 


or 


, , 6. lie 5417 

In p ss In h 19.83 


T 
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or 


In p 


A + 


B 

T 


where A = In (6.11e/u> sw ) + 19.83 and B = —5417. The above 
equation for typical values of w sw and T is nearly a straight line. 
Since w = w aw (T dew ,p), the dew point of unsaturated air defines 
the mixing ratio w at (T,p). It follows that the mixing ratio at 
(T, p) defines the dew point temperature. 

The saturated adiabats are determined numerically from equation 
(7.63). They can be labeled in terms of the estimated d ep or in terms 
of the related (through equation (7.64)) 9„ p . 


9.3 Graphical representation of thermodynamic 
variables in a T — In p diagram 

Figure 9.5 shows the graphical procedure to obtain all the ther- 
modynamic variables derived analytically in Chapter 7, for a skew 
T-ln p diagram. We start with the state (T, p) where the air is 
assumed to be unsaturated and to have a mixing ratio w. The 
intersection between the w = constant line and the temperature 
axis gives the dew point temperature. By following the dry adiabat 
passing through point ( T,p ) (line ai) down to 1000 mb we obtain 
the potential temperature 9. By following it up until we intersect 
the w = constant line we arrive at the lifting condensation level 
(PlclX where the parcel first becomes saturated and has a temper- 
ature T lcl . Through point (T LC l,Plcl) passes a pseudo-adiabat s 


Figure 9.5 

Graphical procedure to 
obtain thermodynamic 
variables, from a skew 
(T-ln p) diagram, given 
an initial state ( T,p ). 


200 mb 


DRY ADIABAT a. 


PSEUDO-ADIABAT s 


p 

1 LCL 

w = CONSTANT — 

A 


V- DRY ADIABAT a 2 

1 000 mb 

/T 

Tlcl T^^ 

T W p i 

l T V\\ 
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defined by 9 ep . If we follow this pseudo-adiabat down to the original 
level p , we obtain the pseudo- wet-bulb temperature, T wp , and con- 
tinuing down to the 1000 mb level we obtain the pseudo-wet-bulb 
potential temperature, 0 wp (this also indicates the one-to-one rela- 
tionship between 9„ p and 9 ep hence the labeling of pseudo-adiabats 
in terms of either 9 ep or 9 wp ). If we follow the pseudo-adiabat s up 
to a pressure level where we can assume that all vapor has con- 
densed and has fallen out of the parcel (typically this level is taken 
to be the 200 mb level), and then follow the dry adiabat passing 
through that point (line a 2 ) down to the original level we obtain 
the pseudo-equivalent temperature, T ep . Continuing down to the 
1000 mb level yields the pseudo-equivalent potential temperature, 
9 ep . Note that T ei ,T w ,0 w and 9 e are not directly estimated in dia- 
grams. As we discussed in Chapter 7, even though they are very 
close to each other, in actuality T ei ^ T ep , T w ^ T wp , 9 W ^ 0 wp , and 
9 ei / 0 ep (recall equation (7.66)). 


9.3.1 Using diagrams in forecasting 

Figure 9.6 shows in a skew T-ln p diagram a sounding taken at 
Davenport, IA, on 21 June 1997 at 0000 UTC. The solid line is the 
vertical profile of environmental temperature and the heavy broken 


Figure 9.6 

Sounding at Davenport, 

IA, on 21 June 1997 at 
0000 UTC. As explained 
in the text, the 
information provided by 
this diagram would not 
support the development 
of deep convection that 
day. Nevertheless, deep 
convection did occur later 
in the day and very severe 
weather resulted in parts 
of the Midwest. See text 
for details. 
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line to its left is the vertical profile of dew point temperature. The 
parcels at the surface are unsaturated. The lifting condensation 
level of these parcels is about 860 mb where T LC L = 20 °C. Above 
the LCL the light broken line is the pseudo-adiabat passing through 
the point (T LCL , Plcl) and below the LCL the light broken line is the 
dry adiabat passing through point (T,p) at the surface. The pseudo- 
adiabat intersects the temperature profile at point A. Up to this 
point the parcel is colder than the environment. As such the layer 
from surface to point A is stable. This means that the buoyancy 
in this layer is negative and that work must be done against the 
buoyancy if parcels originating at the surface are to reach that 
level. The work needed is proportional to the darker shaded area 
(negative area). This work or the magnitude of this negative area 
is often called the convective inhibition energy (CINE) given by 
equation (8.17) for i = p sur f ace and / = p A (here and in practice we 
make no distinction between T and T virt ). 

/ PA 

{T'-T)dlnp. 

surface 

In our example, p sur f ace ~ 980 mb and p A s=s 740 mb. Assuming that 
the average T' — T in this layer is about —0.8 °C and that R s=s i? d , 
we find that CINE ss 64 Jkg -1 . If we now recall equation (8.15) we 
have that 


f PA dv 1,1, 

i’'" 

where v is the vertical velocity. Assuming that v A ~ 0 it follows 
that if a parcel that originated at the surface were to reach point 
A, the minimum initial impulse has to be 

u min = \/2 CINE « 11.4 ms- 1 . 

Therefore, unless an initial impulse of 11. 4ms” 1 is applied to the 
parcels they will not make it to point A. This point is called the 
level of free convection (LFC) . Above this level the parcel is warmer 
than the environment all the way up to 150 mb. The buoyancy from 
LFC to 150 mb is positive and now work is done by buoyancy. This 
work is proportional to the lighter shaded area (positive area). In 
this layer the parcels that made it to LFC are now free to rise 
and accelerate. Note that owing to the temperature inversion in 
the stratosphere, above 150 mb the temperature normally increases 
with height. As such, there is a point ( B in our example) where 
the temperatures of the parcel and the environment become equal 
again. Above this level, which is referred to as the level of neutral 
buoyancy (LNB), work is done against buoyancy and the parcels 
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begin to decelerate. The magnitude of the positive area is called 
the convective available potential energy (CAPE) and is given by 


CAPE = 


-R 


(T' -T ) din p. 


/ LFC 


Assuming that between LFC and LNB the average difference 
between the temperature of the parcels and the environment is 
about 7°C and that R ft R d we can estimate that CAPE ft 
3200 J kg -1 . As we did above, we can show that in this case the 
maximum velocity a parcel will attain is 

i> max = V 2 CAPE sd 80 ms' 1 . 


Both CINE and CAPE are very useful as they provide information 
on whether or not convection will occur (CINE) and on how severe 
a storm might become (CAPE) . 

In our example CAPE is rather high. Thus, if the parcels make it 
to LFC, they will then accelerate greatly and deep convection will 
ensue. The problem in predicting deep convection in this particular 
example, however, is that the value of CINE is considerable. Such 
values may inhibit development of convection. In practice, situa- 
tions of high CAPE values should be monitored to see whether the 
negative area is removed by warming of the air in the lower levels, 
or whether upward acceleration is aided by a low level jet or fron- 
togenesis. In our case all of the above happened thereby allowing 
deep convection and severe weather to develop later in the day in 
the Midwest. 

Together with CAPE and CINE we may calculate the precipitable 
water from soundings. If we recall the definition of precipitable 
water in problem 7.19 and the definition of specific humidity in 
equation (7.2) we have 

1 

d„ = — / q pdz 

Pw Jo 

or using the hydrostatic approximation 

^ /*Psurface 

d„ = / q dp. 

gp w Jo 

The above equation indicates that precipitable water can be esti- 
mated by integrating the specific humidity profile. The specific 
humidity profile as a function of pressure can be obtained from 
the vertical profiles of T and T d using the equations presented in 
Chapter 7. 


Example 

(9.1) An air parcel at 1 000 mb has a temperature of 20 °C and a 
mixing ratio of lOgkgA The parcel is lifted to 700 mb by passing 
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over a mountain. During the ascent 80% of the water vapor 
that condenses falls out of the parcel. Find the temperature and 
potential temperature of the air when it returns to 1 000 mb on 
the other side of the mountain. 

The parcel is located at the point with coordinates T = 
20 °C and p = 1000 mb (w = lOgkg” 1 ). At this point, 
the saturation mixing ratio is ~ 16gkg _1 . Therefore, the 
parcel is unsaturated and as it begins to rise it will fol- 
low the dry adiabat at that point (which because we are 
at 1000 mb is 20 °C). The point at which this dry adia- 
bat intersects the w = lOgkg” 1 line is at p = 910 mb and 
T = T LC l ~ 12.4 °C. This is the lifting condensation level 
and from that point on the parcel will follow the correspond- 
ing pseudo-adiabat (0 ep = 49.5 °C, which in the skew T-ln p 
diagram is labeled in terms of the related 0 wp = 16.5 °C). 
At 700 mb the saturation mixing ratio is ~ 6.66 g kg -1 . 
Therefore, 10 — 6.66 = 3.34 g of vapor per kg of dry air 
have condensed during the saturated ascent from 910 mb to 
700 mb. Since 80% of this amount is removed by precipi- 
tation, it follows that 0.668 g kg -1 of liquid water remains 
in the parcel as it begins its descent. This means that ini- 
tially the descent will be dictated by the same saturated 
adiabat until all the water evaporates, i.e. until the mixing 
ratio increases from 6.66 to 6.66 + 0.668 = 7.328 g kg” 1 . This 
will happen at the point where the pseudo-adiabat inter- 
sects the w = 7.328 line. This is the point with coordinates 
p ss 750 mb and T ss 4°C. From this point on, the descent 
is dictated by the corresponding dry adiabat ( 9 « 28 °C). 
Thus, on the other side of the mountain at p = 1000 mb, 
T = d = 28 °C. 


Problems 

(9.1) Show that the Refsdal diagram with coordinates u = In T, 
w = —Tin p is an area-equivalent diagram and that neither 
its isobars nor its dry adiabats are straight lines. 

(9.2) Show that the Strive (or pseudo-adiabatic) diagram with 
coordinates u = T, w = — p k& is not an area-equivalent dia- 
gram. Show that isobars, isotherms, and dry adiabats are all 
straight lines. 

(9.3) Using the Jacobian method show that the diagram 
with coordinates u = T, w = c p In 6 is an equal-area 
diagram. 

(9.4) Given the sounding shown in the table, 
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p(mb) 

T{° C) 

950 

22.5 

900 

18.0 

850 

15.0 

800 

16.0 

750 

12.0 

700 

7.0 

650 

4.0 

600 

-1.5 

500 

-10.0 

400 

-20.0 


determine using a diagram (1) the mixing ratio, (2) 
the relative humidity, (3) the potential temperature, (4) 
the pseudo-wet-bulb temperature, (5) the pseudo-wet-bulb 
potential temperature. Having determined the above, com- 
pute the wet-bulb temperature, the equivalent tempera- 
ture, the pseudo-equivalent potential temperature, and the 
pseudo-equivalent temperature. Assume that l v = 2.45 x 
10 6 Jkg -1 ,c w = 4218 J kg -1 , and that the dew point at 
950 mb is 15.7°C. (12.4gkg”\ 0.68, 27.5°C, 18°C, 20°C, 
19 °C, 51 °C, 66 °C, 60 °C) 

(9.5) Approximate CAPE if in the previous problem we extend 
the sounding as follows: 


p(mb) 

T(°C) 

300 

-35.0 

200 

-55.0 

150 

-60.0 


Would you expect deep convection to occur? 

(9.6) Consider the following scenario (Figure 9.7). On a given day, 
on a skew T— In p diagram the temperature profile from the 
surface (1000 mb) to 200 mb has a constant lapse rate of 
7 °C km -1 . Above 200 mb the temperature is constant. If the 
parcels at the surface are saturated, what will the maximum 
vertical velocity be on that day assuming that the saturated 
adiabatic ascent can be approximated by the constant lapse 
rate of 6.5 °C km” 1 ? The temperature at the surface is 30 °C. 
Assume R = R d (~ 54ms” 1 ). 

(9.7) A parcel of air with an initial temperature of 15 °C and dew 
point 0 °C is lifted adiabatically from 1000 mb. Determine 
the lifting condensation level and temperature at that level. 
If the parcel is lifted a further 250 mb what will its final 
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Figure 9.7 

Figure for problem 9.6 


p(mb) 


PLNB 


200 


1000 



temperature be and how much liquid water is produced? 
(~ 800 mb, w —3 °C, ~ -22 °C, ~ 2.7gkg" 1 ) 

(9.8) In the previous problem if all the condensed water falls as 
precipitation, can you roughly estimate the height of water 
collected on a rain collector of area 1 m 2 ? Here you have to 
make your own assumptions. How does this result compare 
with the result obtained using the formula in problem 7.19? 
Can you explain any differences? 

(9.9) An air parcel at 900 mb has a temperature of 10 °C and a 
mixing ratio of 5gkg -1 . The parcel is then lifted to 700 mb 
by passing a mountain. If in this process 70% of the con- 
densed water is removed by precipitation, what will its 
temperature, potential temperature, and pseudo- wet-bulb 
potential temperature be when it returns to its initial level 
of 900 mb on the other side of the mountain? (~ 14 °C, 
~ 23.5 °C, ~ 11 °C) 

(9.10) If a temperature profile in °C could be described as T(z) = 
28 — 8 z for z > 1km and T = 20 °C for z < 1km, and 
the mixing ratio at the surface is 16gkg _1 , estimate using a 
diagram the level of free convection. (~ 2.4km) 



CHAPTER TEN 


Beyond this book 


The previous chapters have presented the basics in atmospheric 
thermodynamics. As we know, in atmospheric sciences the ulti- 
mate goal is to predict as accurately as possible the changes in 
weather and climate. Thermodynamic processes are crucial in pre- 
dicting changes in weather patterns. For example, during cloud and 
precipitation formation vast amounts of heat are exchanged with 
the environment that affect the atmosphere at many different spa- 
tial scales. In this final chapter we will present the basic concepts 
behind predicting weather changes. This chapter is not meant to 
treat the issue thoroughly, but only to offer a glimpse of what comes 
next. 


10. 1 Basic predictive equations in the atmosphere 

In the Newtonian framework the state of the system is described 
exactly by the position and velocity of all its constituents. In the 
thermodynamical framework the state is defined by the tempera- 
ture, pressure, and density of all its constituents. In a dynamical 
system such as the climate system both frameworks apply. Accord- 
ingly, a starting point in describing such a system will be to seek 
a set of equations that combine both the mechanical motion and 
thermodynamical evolution of the system. 

The fundamental equations that govern the motion and evolution 
of the atmosphere (and for that matter of the oceans and sea ice) 
are derived from the three basic conservation laws: the conserva- 
tion of momentum, the conservation of mass, and the conservation 
of energy. For the atmosphere the equation of state relates temper- 
ature, density, and pressure. In summary, from the conservation of 
momentum we derive the following set of predictive equations of 
motion (Washington and Parkinson (1986)) 
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with 


du ( tan (f> 
dt-( f+u — 


v = — 


1 1 dp 


a cos (j> p dX 


F x 


dv ( , tan <b 

is + V +u ^r 


1 

u — w— + F^ 

pa o<p 


du> 1 dp 

~rr = — w — <7 + ^2 

dt p az 


d d u d v d d 

dt dt^~ a cos (j) dX a d<j> + dz 


( 10 . 1 ) 

( 10 . 2 ) 

(10.3) 


u = a cos ( 


^d\ 

dt 


v = a 


d(j> 

dt 


dz 

u> = — 
dt 

where u, v, and ui are the horizontal and vertical components of the 
motion, <j) is the latitude, A is the longitude, a is the radius of 
the earth, / is the Coriolis force, and F is the friction force. In 
the above equations the forces that drive the motion are the local 
pressure gradients, gravity, Coriolis, and friction. 

Equation (10.3) can be approximated by the hydrostatic equation 
(by assuming that du>/dt = 0 and F z = 0) 


1 dp 

9 = — IT 
p dz 


(10.4) 


which relates density to pressure. 

From the law of conservation of mass we can derive the continuity 
equation 


dp 

dt 


1 


a cos (j> 


-( PU ) 


dcj) 


(pv cos <f>) 


dz 


(pu>) 


(10.5) 


which provides a predictive equation for density. The last two equa- 
tions in this basic formulation come from thermodynamics. The first 
law of thermodynamics (which expresses the law of conservation of 
energy) 


, dT 1 dp dQ 

p dt p dt dt' 


(10.6) 


where dQ/dt is the net heat gain, provides a predictive equation 
for temperature. The above equations (10.1, 10.2, 10.3 (or 10.4), 
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10.5, and 10.6) make a system of five equations with six unknowns 
(u,v,u,p, p,T). The equation of state 

p = pRT (10.7) 

provides the additional equation that connects pressure, density, 
and temperature, thus resulting in a system of six equations (called 
the primitive equations system), with six unknowns (assuming of 
course that dQ/dt,F x , and F ri , are constants and known). Note, 
however, that this system is not a closed one because dQ/dt, F x , and 
Fjj must be determined from the other variables. These terms are 
very important for climate simulations but for short-term weather 
prediction they are often ignored. 


10.2 Moisture 


The above system of predictive equations was derived without the 
inclusion of moisture. Even though a model can be derived without 
moisture in it, including moisture can dramatically improve mod- 
eling and prediction. In a manner analogous to the continuity of 
mass, the changes of moisture must be balanced by the moisture’s 
sources and sinks. An equation of continuity for water vapor mixing 
ratio can be written as 


dw 

dt 


-M + E 
P 


(10.8) 


where M is the time rate of change of water vapor per unit volume 
due to condensation or freezing, E is the time rate of change of 
water vapor per unit mass due to evaporation from the surface and 
to horizontal and vertical diffusion of moisture occurring on scales 
below the resolution of the model, and w is the mixing ratio. Often 
the above equation is combined with equation (10.5) to obtain 


d(pu) 

dt 


1 

a cos (j> 


d\ 


(puju) 


d(f> 


( pujv cos <f>) 


dz 


(pu>w) 


M+pE. 


(10.9) 

When water vapor changes to water (or ice) and vice versa, heat 
is added to or removed from the atmospheric system. In this case 
dQ = d(lw s /T). 

The above general and simple equations constitute only the start- 
ing point for studying and investigating the fascinating dynamics 
of atmospheric motion. The equations have to be modified and 
adjusted in order to include the effect of many other factors (for 
example, radiative processes) and to be applied over a specific range 
of scales (mesoscale models, for example). Such studies are crucial 
in the effort to understand the climate and ecology of our planet. 
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APPENDIX 


Table A.l. Units 


Physical quantity 

Unit 

MKS system 

egs system 

Length 

m 

cm 

Time 

s 

s 

Velocity 

ms -1 

cms' 1 

Acceleration 

ms” 2 

cms' 2 

Mass 

kg 

g 

Density 

kgm~ 3 

gem 3 

Force 

kg ms -2 (newton, N) 

gems -2 (dyne, dyn) 

Pressure 

kg m 1 s 2 (pascal, Pa) g cm 1 s 2 (microbar, fjb&r) 

Energy 

kgm 2 s 2 (joule, J) 

g cm 2 s 2 (erg) 

Specific energy 

m 2 s“ 2 (Jkg- 1 ) 

cm 2 s- 2 (ergg- 1 ) 


Relations: 

1 N = 10 5 dyn 

1 Pa = 10 /ibar 

1 bar = 10 6 fi bar = 10 5 Pa 

1 atmosphere = 1013 mb = 1013 x 10 2 Pa = 1013 liPa 
1 J = 10 7 erg 
1 cal = 4.185 J 
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Table A. 2. Selected physical constants 


General constants 

Avogadro’s number, N a 

6.022 x 10 23 mol -1 

Boltzmann’s constant, k 

1.38 x 10" 23 JR- 1 

Gravitational acceleration at sea level, g 

9.807ms" 2 

Universal gas constant, R * 

8.314 JK" 1 mol" 1 

Constants for dry air 

Molecular weight, M d 

28.97 g mol" 1 

Density, p d 

1.293 kg m" 3 (at STP)* 

Specific gas constant, R d 

287 J kg" 1 K" 1 

Isobaric specific heat (at 273K),c P d 

1005 J kg" 1 K" 1 

Isochoric specific heat (at 273K),Cyd 

718 J kg" 1 K" 1 

Ratio of specific heats, 7 d 

1.4 

Constants for water 

Molecular weight, M v 

18.015 g mol" 1 

Density (water) 

1000 kg m" 3 (at STP) 

Density (ice) 

917 kg m" 3 (at STP) 
461.51 J kg" 1 K" 1 

Specific gas constant, R v 

e = A/d/JWv = R d / 7? v 

0.622 

Isobaric specific heat (at 273 K) 

1850 J kg" 1 K" 1 

(vapor, c pv ) 

Isochoric specific heat (at 273 K) 

1390 J kg" 1 K" 1 

(vapor, Cvv) 

Specific heat capacity (at 273 K) 

4218 J kg" 1 K" 1 

(water, c pw « Cvw ~ c w ) 

Specific heat capacity (at 273 K) 

2106 J kg" 1 K" 1 

(ice, c pi = Cvi = Ci) 

Latent heats 

see Table A3 


STP: Standard temperature and pressure (273 K, 1013 mb) 
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Table A. 3. Specific latent heats 


T (°C) 

h (lO^kg- 1 ) 

if (HfJkg- 1 ) 

U (HfJkg- 1 ) 

-100 

2.8240 



-90 

2.8280 



-80 

2.8320 



-70 

2.8340 



-60 

2.8370 



-50 

2.8383 

0.2035 

2.6348 

-40 

2.8387 

0.2357 

2.6030 

-30 

2.8387 

0.2638 

2.5749 

-20 

2.8383 

0.2889 

2.5494 

-10 

2.8366 

0.3119 

2.5247 

0 

2.8345 

0.3337 

2.5008 

5 



2.4891 

10 



2.4774 

15 



2.4656 

20 



2.4535 

25 



2.4418 

30 



2.4300 

35 



2.4183 

40 



2.4062 

45 



2.3945 

50 



2.3893 
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Figure A. I 

The skew T— Inp 
diagram. Saturated 
adiabats are labeled in 
terms of 6„ p . 



INDEX 


absolute humidity, 101 
absolute temperature, 15, 60 
Amagat- Andrews diagram, 

82 

Avogadro 

hypothesis, 17 
number, 16 

Boltzmann’s constant, 10 
Boyle’s law, 16 
Brunt-Vaisala frequency, 145 
buoyancy, 144, 169, 170 

Carnot cycle, 55—58 
efficiency, 58 
lessons learned from, 58—62 
chemical potential, 79 
Chinook wind, 136 
Clausius-Clapeyron equation, 
85-87, 107, 114, 117, 
126-127, 134 

and changes in the melting 
and boiling points of 
water, 91 

approximations of, 87-91 
Clausius postulate, 59 
cloud formation 
on ascent, 114 
on descent, 114 
coefficient of thermal expansion 
isobaric, 13 
isocharic, 14 
continuity equation, 176 
convective available potential 
energy (CAPE), 155, 170 
convective inhibition energy 
(CINE), 155, 169 
convective stability, 151 
critical point, 82 

Dalton’s law, 20-21, 99 
degrees of freedom 

in Hamiltonian dynamics, 9 
in thermodynamics, 9 
destabilization of lifted layers, 
149-152 


dew point temperature, 

105-106, 107, 112 
variation (of dew point) with 
height, 106 
disorder, 68-69 
dry adiabatic 

processes, 38-40 
lapse rate, 41-42 
dry air 

composition, 22 
molecular weight, 22 

energy 

conservation of, 5 
definition of, 26 
internal, 5 
kinetic, 5 
potential, 5 
relation between internal 

energy and temperature, 
34 

enthalpy 

definition of, 33 
of fusion, 84 
of sublimation, 84 
of vaporization, 84 

of temperature dependence 
of, 87-88 

entropy 

and atmospheric processes, 
71-72 

and disorder, 68-69 
and irreversible processes, 
60-62 

and potential temperature, 

70 

and reversible processes, 61 
definition of, 62 
for ideal gases, 62 
for liquids and solids, 62 
of an isolated system, 64 
of a mixture, 72-74 
equal distribution of energy 
theorem, 9 


equations of motion in the 
atmosphere, 176 
equivalence between heat and 
work, 31 

equivalent (isobaric) potential 
temperature, 124 
equivalent potential 

temperature, 121 
equivalent (isobaric) 

temperature, 111 
equilibrium 

between vapor and liquid, 

80 

metastable, 3 
stable, 3 
unstable, 3 
vapor pressure 
for ice, 81 
for water, 81 
temperature dependence of 
phase transformation, 
84 

equilibrium state, 3 
exact differential 
definition, 7 
conditions for, 8 

freezing in a cloud, 131-133 
fundamental relations, 66-67 


fusion, 81 
latent heat of, 

84 

gas constant 

for dry air, 22 

for moist air, 

102 

specific, 18 
universal, 19 
Gibbs formula, 

81 

Gibbs function, 

66 

Guy-Lussac 



1 st law, 13 
2 nd law, 14 


heat capacities (see thermal 
capacities) 

Helmholtz function, 66 
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INDEX 


hydrostatic approximation, 

42, 126, 143 
humidity 

absolute, 101 
relative, 101 
specific, 101 

ideal gas 

definition, 10 
law, 12, 17-19 

Joule’s law, 34, 96 

Kelvin’s postulate, 58, 59 
kinetic theory of heat, 9 

latent heat 

definition of, 84 
of fusion, 84 
of sublimation, 84 
of vaporization, 84 
lifting condensation level (LCL) 
definition of, 115 
estimation, 116-117 

Maxwell’s relations, 76 
mixing and irreversibility, 74 
mixing ratio, 101 
saturation, 101 
water, 109 

mixture of gases, 20, 21 
moist adiabatic 
lapse rate, 113 
processes, 109, 112, 113 
moist air 

definition of, 99 
mean molecular weight, 102 
moisture in the equations of 
motion, 177 
molar heat capacity, 33 
mole definition of, 16 
molecular weight of dry air, 

22 

Poisson’s equations, 38-40 
potential temperature, 43-44 
of unsaturated moist air, 

112 

precipitation rate, 138 
precipitable water, 142, 170 
process 

adiabatic expansion of 

unsaturated moist air, 

112 

adiabatic isobaric, 109 
adiabatic isobaric mixing, 

128 

isenthalpic, 109 
isentropic, 119 
isobaric cooling, 104 
pseudo-adiabatic, 122-123 


saturated adiabatic, 

119-122 

saturated ascent, 118 
vertical mixing, 130 
pseudo-equivalent temperature, 
123 

equivalent potential 

temperature, 124 
wet-bulb temperature, 125 
wet-bulb potential 

temperature, 125 

raindrops, growth by 
condensation, 93 
relative humidity, 101 

saturated adiabatic lapse rate, 
126 

saturation temperature (at 
LCL), 115 
specific heat capacity 

at constant pressure, 33 
at constant volume, 33 
specific humidity, 101 
specific gas constant, 18 
stability criteria, 

for saturated parcels, 148, 
149 

for unsaturated parcels, 

147, 149 

of parcels when the layer is 
lifted, 149-152 
stability, convective, 151 
state 

equation of, 2, 12 
functions, 2 
mechanical, 1 
thermodynamical, 2 
variables, 2 
sublimation, 81 
latent heat of, 84 
supercooled water, 81, 94-95 
supersaturation, 129-130 
system 

closed, 1 

homogeneous, 2, 79 
heterogeneous, 2, 79 
isolated, 1 
open, 1, 2 

table of constants, 185 
table of units, 181 
temperature 

dew point, 105-109 
equivalent (isobaric), 111 
equivalent (isobaric) 
potential, 124 
equivalent potential, 121 
Kelvin scale, 15 
potential, 43-44, 112 
of unsaturated moist air, 
112 


pseudo-equivalent, 123—124 
pseudo-equivalent potential, 

124 

pseudo-wet-bulb, 125 
pseudo-wet-bulb potential, 

125 

relationship between dew and 
wet-bulb, 112 
saturation (at LCL), 115 
thermodynamic definition of, 
67 

virtual, 102-103 
virtual potential, 113 
wet-bulb, 109-111 
wet-bulb potential, 125 
temperature lapse rate 
dry adiabatic, 41-43 
moist adiabatic, 113 
saturated adiabatic, 

126-127 

thermal (heat) capacity 
at constant pressure, 33 
for dry air, 36 
for moist air, 103-104 
at constant volume, 33 
of diatomic gases, 36 
of monatomic gases, 35 
thermodynamics 
definition of, 1 
statistical nature of, 67-70 
1 st law of, 32, 37 
2 nd law of, 62, 64 
combining the two laws, 66 
thermodynamic diagrams 
and forecasting, 168-170 
conditions for area-equivalent 
transformation, 

159-161 

emagram, 163 
graphical representation in, 
167 

skew emagram (T-lnp 
diagram), 165 
tephigram, 161, 162 
transformations 


adiabatic, 

4, 38, 65 

cyclic, 4, 

38 

irreversible 

3, 8 

isentropic, 

65, 71 

isobaric, 

4, 37, 65 

isochoric, 

4, 37, 65 

isothermal, 

4, 37, 65 

polytropic, 

40-41 

reversible, 

3, 8 

triple point, 

82 


units, 181 

universal gas constant R* , 
19 


vaporization, 81 
latent heat of, 84 
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variables 

dependent, 2 
extensive, 2 
independent, 2 
intensive, 2 

virtual potential temperature, 
113 


virtual temperature, 102-103 
water 

density of, 82 
mixing ratio, 109 
precipitable, 142, 170 

supercooled, 81, 94-95 

transformations of, 80-82 


wet-bulb temperature, 

109-111 

wet-bulb potential temperature, 
125 
work 

and kinetic energy, 29 
definition of, 27 



